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DOCUMENTATION  OF  WAVE-HEIGHT  AND  TIDAL  ANALYSIS  PROGRAMS 
FOR  AUTOMATED  DATA  ACQUISITION  AND  CONTROL  SYSTEMS 


PART  I :  INTRODUCTION 
Background 

1.  Over  the  past  decade,  automated  processing  technology  has  evolved 
from  large,  expensive  computers  to  minicomputers  and  microprocessors.  With 
this  evolution,  the  physical  size  and  cost  of  automated  processing  systems 
have  greatly  decreased  with  a  minimal  decrease  in  system  capabilities.  Cost 
reductions  for  such  systems  have  resulted  in  economic  justification  of  the  use 
of  minicomputers  and  microprocessors  to  a  specific  task  or  group  of  specific 
operations,  whereas  large  computers  can  be  justified  economically  only  for 
multiple  operations  and  tasks.  The  automation  of  physical,  hydraulic  modeling 
techniques  has  lagged  automation  efforts  in  many  other  fields  mainly  because 
of  cost  justification  and  the  requirements  of  highly  specialized  instrumenta¬ 
tion  (e.g.,  sensors).  However,  needs  for  such  automation  have  existed  for 
many  years.  These  needs  are  the  results  of  requirements  for  (a)  real-time 
model  control  decisions,  (b)  quasi-real-time  data  analyses,  and  (c)  more 
accurate  and  reliable  model  data  for  engineering  and  environmental  interest 
studies . 

2.  In  the  Wave  Dynamics  Division  (WDD)  of  the  Hydraulics  Laboratory, 

U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  recent  increases  in 
the  number  and  complexity  of  physical  model  studies  conducted  each  year  and 
the  use  of  such  models  to  solve  hydraulic  problems  of  extremely  large  harbors 
(Whalin  et  al.  1974,  Outlaw  et  al.  1977)  and  estuaries  (requiring  large  model 
areas)  have  vividly  demonstrated  a  need  for  automation  of  modeling  procedures 
including  operation,  data  collection,  and  data  analysis.  The  physical  size  of 
models  involving  the  study  of  long-wave  phenomena,  the  vast  complex  of  basins 
and  channels  requiring  detailed  study,  and  the  large  number  of  tests  necessary 
to  evaluate  the  effect  of  improvement  plans  on  harbor  circulation  and  oscilla¬ 
tion  eliminate  manual  procedures  for  data  collection  and  analysis.  Over  the 
last  several  years  the  WDD  has  been  very  successful  in  automating  (Durham  and 
Downing  1978)  major  aspects  of  its  physical  models  for  wave  (Durham  and  Greer 
1975)  and  tidal  inlet  studies  (Durham,  Greer,  and  Whalin  1976).  Major 


automation  efforts  were  devoted  to:  (a)  model  control,  (b)  model  data  acqui¬ 
sition,  and  (c)  model  data  reduction  and  analyses.  The  two  automated  systems 
of  WDD  have  been  given  the  title  "Automated  Data  Acquisition  and  Control  Sys¬ 
tem,"  whose  acronym  is  ADACS.  Although  each  physical  model  has  some  unique 
requirements,  many  automation  requirements  are  similar  and  are  shared  by  all 
physical  models.  Therefore  ADACS  for  the  wave  and  tidal  inlet  models  ol  WDD 
were  developed  with  as  much  generality  and  flexibility  designed  into  each 
system  as  was  possible  considering  the  specific  requirements  of  various  model 
Such  approach  to  the  design  of  each  ADACS  provides  each  system  with  the  capa¬ 
bility  of  being  used  by  future  physical  models  and  sharing  transducers  for 
measuring  various  model  parameters  and  for  making  major  modifications  to  each 
system.  The  general  configuration  of  each  ADACS  consists  of  the  following 
four  subsystems:  (a)  a  data  recording  and  control  generating  subsystem  which 
is  basically  a  minicomputer  with  required  peripheral  devices  and  appropriate 
interfacing  to  other  subsystems;  (b)  model  sensors  and  interfacing  equipment 
for  measuring  model  parameters  such  as  water  velocity,  temperature  of  water 
and/or  air,  changes  in  water-surface  elevations  (waves  and/or  tides),  etc.; 
(c)  model  controls  and  interfacing  equipment  for  wave  and  tide  generators, 
sensor  calibrations,  etc.;  and  (d)  a  data  reduction,  analysis,  and  display 
subsystem  which,  in  most  cases,  is  the  minicomputer  in  the  first  subsystem 
with  mass  storage  and  display  devices.  The  design,  development,  and  general 
application  of  ADACS  to  wave  and  tidal  inlet  models  are  presented  in  detail 
by  Whalin  et  al.  (1974),  Outlaw  et  al.  (1977),  Durham  and  Downing  (1978), 
Durham  and  Greer  (1975),  and  Durham,  Greer,  and  Whalin  (1976),  which  are 
suggested  to  readers  requiring  detailed  information  on  ADACS.  However,  for 
general  information  as  to  system  configuration  and  subsystem  makeup,  a 
schematic  of  ADACS  for  wave  and  tide  models  is  presented  in  Figure  1  for 
reference  during  reading  of  this  report. 


Data  Analyses 

3.  Digital  computer  programs  have  been  written  to  process  hydraulic 
model  wave  and  tide  data  which  are  collected  by  the  ADACS  and  are  time-based 
voltages  representing  water-surface  elevations.  In  addition,  velocity 
measurements  (time-based  voltages)  are  collected  and  processed  by  ADACS. 

These  programs  accept  model  and  testing  parameters  as  well  as  model  wave  data 


from  a  magnetic  tape  generated  by  ADACS .  ADACS  computer  programs  for  the 
acquisition  of  wave  and  tide  data  and  control  of  various  model  operations  are 
documented  elsewhere. 

4.  Previous  model  data  reduction  techniques  were  time-consuming,  sub¬ 
jective  analyses  performed  by  individuals.  These  manual  techniques  were  con¬ 
fined  to  analyzing  analog  wave  and  tide  records.  For  wave  data,  significant 
wave  heights  of  progressive  waves  and  heights  and  periods  of  standing  waves 
were  obtained.  For  tide  data,  mean  tide  level,  tidal  ranges,  and  phases  were 
obtained.  With  the  development  and  installation  of  ADACS  on  hydraulic  models 

i  an  opportunity  and  need  existed  to  automate  such  analysis  techniques.  Such 

automated  procedures  provide  uniformity  and  accuracy  of  results,  allow  so¬ 
phisticated  mathematical  and  statistical  procedures  to  be  employed,  expedite 
the  reduction  of  large  volumes  of  data,  and  provide  model  results  in  a  report 
quality  format. 

5.  Data  analyses  can  be  performed  by  either  the  ADACS  or  the  Honeywell 
DPS-1  of  the  Automatic  Data  Processing  Center  at  WES.  Schematically,  the 
various  procedures  for  wave  and  tidal  data  analyses  are  as  follows: 

a.  Program  initialization. 

(1)  Input  test  parameters  and  option  flags. 

(2)  Read  and  decode  data  files. 

(3)  Demultiplex  data  files  and  scale  data. 

b.  Wave  record  analyses. 

(1)  H  ±  oH  and  T  ±  a T  .* 

(2)  H 

rms 

(3)  where  x  is  a  specified  percent  of  the  highest  wave 

heights,  normally  x  =  0.333  ;  thus  is  significant 

wave  height. 

(4)  Option  to  plot  wave  heights  versus  time. 

c.  Least-squares  harmonic  analysis  for  tide  records. 

(1)  Amplitude  and  phase  for  specified  periods. 

(2)  Relative  phases  and  amplification  factors  between  gages. 

(3)  Analysis  of  residual  variance. 

(4)  Graphic  output  of  above  results. 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  A). 
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6.  The  major  purpose  of  this  report  is  to  document  the  ADACS  computer 
programs  that  were  developed  for  wave  and  tidal  record  analyses.  The  wave 
record  analysis  program,  which  is  routinely  used  to  analyze  all  model  wave 
data,  and  tide  record  analysis  programs,  used  in  analyzing  model  tide  data, 
are  documented  in  detail  in  this  report.  Required  inputs  for  each  program  and 
various  tabular,  graphic,  and  magnetic  tape  output  options  are  fully  described 
and  illustrated.  These  computer  programs  are  designed  to  be  used  (a)  by  ADACS 
configured  as  described  in  this  report  (Figure  1),  and  (b)  for  analyses  of 
wave  and  tide  data  collected  from  all  physical  wave  models  in  WT" 


PART  II:  WAVE-HEIGHT  PROGRAM 


Program  Functions 

7.  This  program  can  process  up  to  50  channels  of  wave  data  which  are 
recorded  on  9-track  magnetic  tape  by  ADACS.  Wave-height  estimates  and  domi¬ 
nant  wave  period  are  computed  for  each  channel  of  data  with  options  to  list 
in  several  tabular  forms  these  results  and  various  other  data  statistics. 
These  wave-height  estimates  can  be  recorded  in  permanent  data  files  on  a  mag¬ 
netic  tape  for  later  reference  or  generation  of  report  tables,  etc.  Also, 
plots  of  the  scaled  data  with  wave  peaks  and  troughs  labeled  can  be  generated 
for  each  data  channel. 

Program  Specifics 


Crest-trough  search  techniques 

8.  Wave  data  taken  by  ADACS  consist  of  a  number  of  samples  of  water- 
surface  elevations,  f(NAt)  ,  taken  at  equally  spaced  time  intervals,  At  , 
for  N  samples  with  n  =  1,...,N  .  Thus  the  sampling  duration  or  record 
length  is  NAt/T  wave  periods  and  T/At  samples  per  wave  period  in  the 
record.  Estimations  of  these  parameters  assume  that  the  waves  for  each  wave 
gage  or  data  channel  are  monochromatic. 

9.  A  simple  procedure  for  analyzing  this  wave  record  is  to  bracket  the 
elevation  data  in  each  of  the  NAt/T  data  segments  or  windows  and  search 
these  windows  for  maximum  and  minimum  water-surface  elevations  (crests  and 
troughs).  The  most  difficult  part  of  this  procedure  is  choosing  the  windows 
correctly. 

10.  Starting  at  the  beginning  of  the  record  it  is  possible  to  bracket 
each  of  NAt/T  wave  periods  successively  to  find  the  maximum  elevations  dur¬ 
ing  each  period.  Before  the  crests  and  troughs  are  used  to  define  a  wave 
height,  they  are  checked  for  acceptability  according  to  the  following  two 


conditions : 


a.  Each  crest  follows  a  trough. 

b.  Computed  time  intervals  between  deviations  of  crest  to  crest 
and  trough  to  trough  are  less  than  a  specified  variation  in  the 
period  of  the  generated  wave. 
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A  general  flow  chart  of  the  wave  peak-trough  search  routine  is  given  in  Fig¬ 
ure  2.  Using  these  selected  wave  heights  (crest-trough),  various  statistics 
(mean,  root  mean  square,  and  mean  of  largest  one-third  wave  heights)  are  cal¬ 
culated.  The  positions  relative  to  time  of  the  crests  and  troughs,  which  are 
used  to  define  the  wave  heights,  are  also  used  to  calculate  a  dominant  wave 
period  of  the  record. 

Scaling  of  data 

11.  All  values  output  by  this  program  are  scaled  to  prototype  feet  by 
using  either  linear,  quadratic,  or  spline  interpretation  schemes  based  on  a 
spline  fitting  of  the  calibration  data.  However,  as  mentioned  previously,  the 
data  as  recorded  on  magnetic  tape  are  voltage  readings  in  either  integer  or 
scaled  fraction  form  which,  when  scaled  by  the  chosen  technique,  represent 
water-surface  elevations.  The  analysis  procedure  mentioned  previously  is 
applied  to  the  unsealed  data  (voltages).  Since  the  analysis  procedure  is 
searching  only  for  maximums  and  minimums,  the  procedure  will  work  on  any  type 
of  data  whether  the  data  are  voltages  or  water-surface  elevations.  Extreme 
voltages  also  are  the  extreme  water-surface  elevations  when  the  voltages  are 
scaled . 

12.  The  analysis  is  done  on  the  unsealed  data  for  the  following  two 
reasons : 

a.  Integer  voltages  require  one-half  of  the  memory  required  for 
real  scaled  water-surface  elevations. 

b.  ADACS  can  execute  integer  arithmetic  faster  than  real  arithme¬ 
tic.  Thus,  by  performing  the  analysis  on  the  unsealed  data, 
computer  memory  and  computer  run  time  are  minimized. 

All  values  to  be  printed  or  saved  on  magnetic  tape  are  scaled  to  water-surface 
elevations.  The  scaling  technique  and  scaling  coefficients  are  input  from 
magnetic  tape  files  generated  during  data  acquisition. 

13.  Data  are  scaled  by  calibration  coefficients  located  in  the  calibra¬ 
tion  coefficient  record  (Appendix  C)  on  magnetic  tape.  Flags  located  in  the 
transducer  flag  record  on  magnetic  tape  contain  information  yielding  the  type 
of  fitting  procedure  to  use  with  the  calibration  coefficients.  As  mentioned 
previously,  the  three  fitting  procedures  are  linear,  quadratic,  and  spline. 

The  following  equation  is  used  for  both  linear  and  quadratic  fitting 
procedures 


SE  =  ( (SF3  x  USE  +  SF2)  x  USE  +  SF1]  x  MSF  +  MWL 


.  f * 
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Figure  2.  General  flow  chart  of  wave  peak- trough  search  routines 


SE  =  scaled  water  elevation 
USE  =  unsealed  water  elevation 
SF1  =  scale  factor  constant 

SF2  =  scale  factor  linear  term 

SF3  =  scale  factor  quadratic  term 

MSF  =  model  scaling  factor 

MWL  =  reference  mean  water  level 

This  one  equation  can  be  used  for  both  linear  and  quadratic  fits  since  SF3 
will  be  zero  when  the  linear  fit  is  used.  SF1,  SF2,  and  SF3  come  from  the 
calibration  coefficients  record  on  magnetic  tape.  For  the  cubic  spline  fit¬ 
ting  (Cheek,  Radhakrishnan,  and  Tracy  1971)  procedure,  the  necessary  inputs 
include  the  actual  observations  from  which  the  spline  coefficients  were  de¬ 
rived.  This  information  is  stored  in  the  calibration  voltage  record  (Appen¬ 
dix  C).  Thus  it  is  necessary  to  provide  some  information  regarding  the 
calibration  of  the  sensor. 

14.  In  order  to  calibrate  each  sensor,  the  output  voltage  from  the 
sensor  is  monitored  and  recorded  as  the  sensor  is  moved  vertically  a  known 
distance  into  or  out  of  the  still  water.  A  precision,  linear-position 
potentiometer  is  located  on  the  sensor  stand  and  is  coupled  directly  to  the 
sensor  by  a  gear-train  driven  by  an  electric  motor.  By  moving  vertically  the 
coupled  sensor  and  potentiometer  wiper  with  the  electric  motor  and  by  moni¬ 
toring  the  output  voltage  from  the  potentiometer,  the  sensor  can  be  moved 
vertically  a  precise  distance.  By  systematically  moving  the  sensor  through 
11  quasi-equally  spaced  locations  over  the  range  being  used,  21  voltage  sam¬ 
ples  are  obtained  from  which  an  averaged  voltage  value  for  each  of  the 
11  locations  is  calculated  (Figure  3).  This  averaging  technique  minimizes 
the  effects  of  any  slack  in  gear  drives  and  any  hysteresis  of  the  sensors  re¬ 
sponse.  The  21  voltage  readings  taken  from  both  the  sensor  and  the  potenti¬ 
ometer  wiper  are  recorded  on  magnetic  tape  in  the  calibration  voltage  record. 
The  11  averaged  readings  are  fitted  to  a  least  squares  linear  fit.  If  the 
maximum  deviation  from  the  fit  is  greater  than  a  previously  specified  toler¬ 
ance,  a  least  squares  quadratic  fit  is  performed.  If  the  maximum  deviation 
from  this  fit  is  greater  than  a  previously  specified  tolerance,  the  cubic 
spline  fit  is  computed.  If  the  spline  fit  fails  a  specified  tolerance  for  the 
slopes  of  the  end  points,  the  previously  calculated  quadratic  fit  is  used  and 
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Figure  3.  Waverod  calibration 

a  message  to  this  effect  is  supplied  to  the  operator. 

15.  In  computing  the  information  that  the  spline  fit  routine  must  use, 
the  21  voltage  readings  from  both  the  sensor  and  the  potentiometer  wiper  must 
be  averaged  for  the  11  locations  recorded  during  the  sensor's  calibration.  In 
turn,  the  voltages  recorded  from  the  potentiometer  wiper  must  be  scaled  to  a 
displacement  by  the  proper  potentiometer  calibration  coefficient  recorded  in 
the  reference  potentiometer  calibration  coefficient  record  (see  Appendix  C) . 
Thus,  for  each  displacement  that  is  scaled  from  the  potentiometer  wiper  volt¬ 
ages,  a  corresponding  voltage  from  the  sensor  is  available.  This  is  the  in¬ 
formation  the  spline  fitting  procedure  must  have. 

In-line  assembly  language 

16.  The  EAI  Pacer  100  FORTRAN  compiler  allows  the  programmer  to  use 
some  EAI  Pacer  100  assembly  language  instructions  in  line  with  FORTRAN  coding. 
The  only  three  assembly  language  instructions  used  in  this  program  are  the 
Jump,  Link,  and  Octal  instructions. 

17.  The  Octal  instruction  is  used  in  this  program  simply  to  set  aside 
one  memory  cell  and  has  the  following  form: 

123  OCT  0 


where  123  is  a  statement  number. 
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18.  The  Link  instruction  is  very  much  like  the  FORTRAN  CALL  instruction 
and  has  the  following  form: 

L  .  123 

where  123  is  the  statement  number  upon  which  the  link  instruction  will  act 
and  is  always  an  Octal  instruction  in  this  program.  If  the  instruction 

L  .  123 

were  executed,  program  execution  would  transfer  to  the  statement  following 
statement  number  123  and  a  return  pointer  would  be  stored  in  statement  number 
123. 

19.  The  Jump  instruction  is  used  in  the  following  two  forms  in  this 
program: 

J  .123 
J  .123,2 

where  123  is  the  statement  number  to  be  acted  upon.  In  the  first  form 
(J  .123),  the  Jump  instruction  is  used  exactly  like  the  FORTRAN  "GO  TO" 
statement.  That  is 

J  . 123  and 
GO  TO  123 

both  cause  program  execution  to  transfer  to  statement  number  123.  The  latter 
form  (J  .123,2)  is  used  similarly  to  the  FORTRAN  "RETURN"  statement.  After 
a  Link  (L  .123)  instruction  has  been  executed,  a  return  pointer  is  stored  in 
statement  number  123.  A  Jump  instruction  of  the  following  form 

J  .123,2 

causes  program  execution  to  transfer  to  the  return  pointer  stored  in  state¬ 
ment  number  123. 


20.  In  the  following  sample  program,  all  forms  of  the  three  assembly 
language  instruction  (Link,  Jump,  and  Octal)  are  demonstrated.  Program 
execution  follows  sequentially  the  statement  numbers  associated  with  each 
statement . 

1  L  .2 

5  J  .6 

2  OCT  0 

3  1  =  1 

4  J  .2,2 

6  END 


21.  Several  program  initializations  are  input  by  teletype.  The  func¬ 
tions  and  pertinent  related  information  for  each  of  these  inputs  are  discussed 
in  this  section. 

a.  Tape  Unit=?  -  The  magnetic  tape  unit  on  which  the  data  tape 

is  located  is  input  as  either  14  or  15  for  an  EAI  Pacer  100 

with  two  tape  drives. 

b.  List  Gage  Numbers?  -  A  "YES"  response  is  input  to  cause  each 

gage's  number  to  be  listed  to  the  teletype  as  its  data  are 

processed.  This  is  done  so  that  the  operator  can  monitor  pro¬ 

gram  progression  and  alter  sense  switch  settings  for  different 
gages  if  desired.  A  "NO"  response  suppresses  gage  number 
listing. 

£•  Title=?  -  If  a  title  for  the  tabular  output  and  plots  is  de¬ 
sired,  it  is  input.  A  maximum  of  40  characters  may  be  input. 

d.  Data  Scaled  Fractions?  -  Data  as  recorded  on  magnetic  tape  can 
be  in  either  integer  or  scaled  fraction  form.  Integers  will  be 
ranged  from  ±10,000  representing  millivolts.  Scaled  fractions 
also  will  be  ranged  from  ±10,000  millivolts  but  their  bit  con¬ 
figuration  is  much  different  than  that  for  integers.  Voltages 
as  read  by  the  EAI  Pacer  100  always  will  be  in  scaled  fraction 
form.  However,  since  integers  are  more  commonly  used  and  are 
more  readily  understood,  the  analysis  procedure  expects  data 
to  be  in  integer  form.  If  this  input  is  "YES,"  the  analysis 
program  converts  the  data  from  scaled  fraction  form  to  integer 
form.  If  the  input  is  "NO,"  the  analysis  program  expects  the 
data  to  have  already  been  converted  to  integer  form. 


£■  Average  Model  Depth=?  -  The  average  model  depth  is  input  as 

model  feet.  This  is  used  in  computing  model  scale  factors  for 
a  distorted  model . 

f-  File  Name=?  -  The  data  file  name  to  analyze  is  input  as  ASCII 
characters . 

£•  //PER  to  Analyze  and  //REC  to  Skip=?  -  The  number  of  wave  periods 
to  analyze  and  the  number  of  magnetic  tape  data  records  to  skip 
before  the  analysis  begins  are  input. 


Sense  switch  settings 

22.  Sense  switches  are  single  throw  switches  located  on  the  control 
panel  of  the  Pacer  100.  They  are  manually  set  and  can  be  read  by  the  CPU. 
The  FORTRAN  statement,  IF  (SENSW(l))  GO  TO  10,  causes  program  execution  to 
branch  to  statement  number  10  if  sense  switch  A  (1)  is  set  or  is  .TRUE.  It 
follows  that  sense  switches  are  treated  as  logical  variables. 

23.  The  advantage  to  using  sense  switches  is  that  they  may  be  set  or 
reset  while  a  program  is  executing  without  impeding  program  execution.  The 
functions  of  the  sense  switches  in  this  program  are  as  follows: 


a.  Sense  Switch  A  [SENSW(l)]  -  If  sense  switch  A  is  set,  the 

"limited"  printer  output,  described  in  paragraphs  29-35  of  the 
Program  Outputs  section,  is  generated. 


b.  Sense  Switch  B  [SENSW(2)1  -  If  sense  switch  B  is  set,  the 

operator  may  change  the  number  of  data  cycles  to  analysis  and 
also  may  skip  data  records  on  magnetic  tape  before  the  analysis 
begins . 


c.  Sense  Switch  C  [SENSW(3)]  -  If  sense  switch  C  is  set,  a  plot 
of  the  scaled  prototype  data  is  generated. 

d.  Sense  Switch  D  [SENSW(4)]  -  If  sense  switch  D  is  set,  the  orig¬ 
inal  data  are  copied  to  a  disc  file  and  demultiplexed  from  the 
disc  file.  This  is  done  to  speed  program  execution. 

e.  Sense  Switch  E  [SENSW(5))  -  If  sense  switch  E  is  reset,  this 
program  analyzes  consecutive  data  files-  until  two  consecutive 
end  of  file  marks  are  encountered  or  until  the  end  of  tape 
mark  is  encountered.  If  sense  switch  E  is  set,  program  execu¬ 
tion  is  halted  upon  completion  of  the  current  data  file  being 
analyzed . 


f.  Sense  Switch  F  [SENSW(6)]  -  If  sense  switch  F  is  set,  magnetic 
tape  output,  described  in  paragraphs  36-38  of  the  Programs  Out¬ 
put  section,  is  generated. 

g.  Sense  Switch  G  ( SENSW ( 7 ) )  -  If  sense  switches  G  and  A  are  set, 
the  summary  output,  described  in  paragraphs  29-35  of  the  Pro¬ 
gram  Outputs  section,  is  generated. 

h.  Sense  Switch  H  [SENSW(8)]  -  If  sense  switch  H  is  set,  analysis 
of  the  next  data  gage  is  omitted.  By  exercising  the  option, 


described  in  paragraph  21  of  the  Program  Inputs  section  of 
listing  gage  numbers,  program  progression  can  be  monitored  and 
any  gage’s  analysis  can  be  omitted  by  setting  this  sense  switch. 


Magnetic  tape 

24.  Testing  parameters,  gage  identification,  scaling  factors,  test 
identification,  and  test  data  are  all  recorded  or  either  derived  from  infor¬ 
mation  recorded  on  magnetic  tape.  A  general  tape  format  is  described  in  this 
section  with  more  specific  information  available  in  Appendix  C.  The  data 
tests  are  recorded  in  files  on  magnetic  tape.  There  are  two  general  file 
formats  containing  a  variable  number  of  records. 

25.  The  two  file  formats  can  be  called  a  calibration  file  and  a  non¬ 
calibration  file.  In  a  calibration  file,  several  records  are  devoted  to  gage 
identification  and  data  scaling  coefficients.  Since  this  information  remains 
constant  for  a  series  of  data  tests,  there  is  no  need  to  record  it  in  every 
data  file.  A  calibration  file  and  a  noncalibration  file  are  similar  except 
that  some  records  included  in  the  calibration  file  are  excluded  from  the 
noncalibration  file.  Thus  this  exclusive  information  recorded  in  a  calibra¬ 
tion  file  is  applicable  to  all  succeeding  data  files  until  updated  by  another 
calibration  file. 

26.  Each  file  is  separated  by  an  end-of-file  (EOF)  mark  and  the  file 
formats  are  as  follows: 

a.  File  Identification  Record  -  This  record  is  five  16-bit  words 
long.  The  first  three  words  contain  a  six-character  ASCII  code 
denoting  the  file  name.  The  fifth  word  contains  the  length  of 
the  next  record  (usually  20  words). 

b.  First  File  Header  Record  -  Testing  and  model  identification 
information  is  contained  in  this  record.  The  15th  element  of 
this  record  determines  whether  this  is  a  calibration  file  or  a 
noncalibration  file.  The  20th  word  of  this  record  contains  the 
length  of  the  next  record  (usually  50  words). 

c.  Second  File  Header  Record  -  This  record  contains  testing  param¬ 
eters  and  model  control  information.  The  model  control  infor¬ 
mation  is  not  used  in  reducing  the  data. 

d.  Transducer  Flag  Record  -  This  record  contains  two  flags  for 
each  sensor  sampled  during  data  acquisition.  These  flags  con¬ 
tain  data  scaling  information,  gage  type,  gage  number,  and  in¬ 
formation  revealing  calibration  coefficient  record  and  calibra¬ 
tion  voltage  record  lengths.  This  record  exists  only  in  a 
calibration  file. 

e.  Calibration  Coefficient  Record  -  This  record  length  is  deter¬ 
mined  by  the  number  of  gages  having  calibration  coefficients 


associated  with  them.  A  flag  in  the  Transducer  Flag  Record 
determines  whether  a  gage  has  coefficients  in  the  Calibration 
Coefficient  Record.  There  are  11  coefficients  for  each  gage 
having  coefficients  recorded  in  this  record.  If  no  gages  have 
entries  in  this  record,  the  record  is  not  recorded  on  tape. 

This  record  exists  only  in  a  calibration  file. 

f.  Gage  Name  Record  -  This  record  contains  two  elements  (four  ASCII 
characters)  for  each  gage.  These  ASCII  characters  are  used  for 
gage  labeling  and  identification.  This  record  exists  only  in  a 
calibration  file. 

g.  Weighted  Mean  Depth  Record  -  This  record  contains  an  entry  for 
each  wave  gage  representing  the  average  model  depth  in  feet 
between  the  gage  and  the  wave  generator.  These  values  are  used 
in  computing  Keulegan's  friction  coefficients  (Keulegan  1950). 
This  record  exists  only  in  a  calibration  file. 

h.  Gage  Distance  Record  -  This  record  contains  an  entry  for  each 
wave  gage  representing  the  distance  between  each  gage  and  the 
wave  generator  in  prototype  feet.  This  record  exists  only  in 
a  calibration  file. 

i.  Gage  Coordinate  Record  -  This  record  contains  three  floating 
point  coordinates  for  each  gage.  An  X  and  Y  coordinate  deter¬ 
mining  location  and  a  Z  coordinate  representing  water  depth  for 
each  gage  are  recorded.  This  record  exists  only  in  a  calibra¬ 
tion  file. 

Reference  Potentiometer  Calibration  Coefficient  Record  -  This 
record  contains  the  coefficients  necessary  to  convert  the 
millivoltage  readings  collected  from  the  reference  potentiom¬ 
eters  (discussed  in  PART  III)  to  inches.  There  is  one  co¬ 
efficient  (slope)  for  each  of  the  gages  calibrated.  This 
record  exists  only  in  a  calibration  file. 

k.  Calibration  Voltage  Record  -  The  21  voltage  samples  collected 
from  both  the  reference  potentiometer  and  the  wave  sensor  dur¬ 
ing  calibration  (paragraph  14)  are  recorded  in  this  record. 

These  actual  readings  are  necessary  if  the  data  scaling  re¬ 
quires  that  a  cubic  spline  technique  be  used.  This  record 
exists  only  in  a  calibration  file. 

l.  Data  Records  -  The  remaining  records  are  data  records.  The 
length  of  these  records  is  equal  to  the  number  of  data  gages 
times  the  number  of  data  samples  per  wave  cycle  divided  by  the 
numl  er  of  data  records  per  wave  cycle.  The  number  of  data 
records  is  equal  to  the  number  of  wave  cycles  sampled  times 
the  number  of  records  per  wave  cycle. 


Program  Outputs 


All  teletype  outputs  generated  by  this  program  are  to  aid  in 


program  operation.  The  teletype  inputs  described  in  paragraph  21  are  input 
upon  request  from  teletype  outputs;  that  is,  each  teletype  input  is  requested 
by  a  teletype  output  in  a  question-answer  fashion.  The  only  other  teletype 
outputs  that  may  be  generated  by  this  program  are  gage  numbers.  If  requested, 
as  each  gage  is  processed  the  gage's  number  is  output  to  the  teletype.  This 
would  be  done  so  that  the  operator  could  monitor  program  progression  and  alter 
sense  switch  settings  for  different  gages  if  desired. 

28.  Examples  of  teletype  output  as  well  as  teletype  input  are  illus¬ 
trated  in  Plates  1  and  2.  Notice  that  all  teletype  outputs  except  gage  number 
listings  in  Plate  2  are  followed  by  question  marks.  All  other  listings  repre¬ 
sent  teletype  inputs.  Notice  also  that  the  analysis  in  Plate  2  will  be  for 

20  wave  periods  with  the  first  two  data  records  on  magnetic  tape  skipped  or 
not  analyzed. 

Printer 

29.  There  are  three  general  forms  of  printer  output  for  this  program. 
The  output  can  vary  from  very  descriptive  analysis  to  little  more  than  a  sum¬ 
mary  of  the  analysis  results.  The  form  of  output  desired  is  selected  by 
sense  switch  options.  For  this  documentation  the  three  forms  of  output  will 
be  referred  to  as  full  (Plate  3),  limited  (Plate  4),  and  summary  (Plate  5) 
outputs . 

30.  The  first  lines  of  output  generated  by  the  program  (Plate  6)  are 
test  identification  information.  This  is  information  that  comes  from  the 
first  and  second  header  records  on  magnetic  tape.  For  the  full  and  limited 
output  options,  this  information  is  output  on  a  separate  page.  For  the 
summary  output,  the  test  identification  is  on  the  same  page  as  the  analysis 
output . 

31.  For  the  full  output  analysis,  a  page  is  devoted  to  each  gage  with 
gage  number,  gage  arrangement,  prototype  water  depth,  and  gage  coordinates 
printed  out.  The  crest-trough  search  technique  itemizes  the  crests  and 
troughs  found,  in  prototype  minutes,  from  the  start  of  the  test.  The  wave- 
height  value  and  time  also  are  listed.  Statistics  in  both  model  inches  and 
prototype  feet  of  the  data  are  then  printed  out  which  include  average  wave 
height  with  standard  deviation,  rms ,  average  of  highest  one-third  wave  heights 
with  standard  deviation,  average  of  highest  one-third  wave  heights  with 
Keulegan's  (1950)  friction  coefficients  applied,  and  computed  wave  period  with 
standard  deviation. 
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32.  There  are  two  error  messages  that  the  search  routine  may  output. 

The  first  is  shown  in  Plate  7  as  "***  only  23  periods  of  Data  Found 

This  is  not  a  critical  situation  since  the  program  was  only  searching  for 
24  wave  periods.  In  fact,  it  is  quite  common  to  find  one  or  two  less  periods 
than  searching  for  due  to  the  nature  of  the  gage  arrangement  in  the  model, 
the  data  acquisition  technique,  and  the  search  technique.  However,  the  other 
error  situation  may  be  more  critical.  Plate  8  illustrates  a  search  error.  A 
search  error  simply  means  that  a  crest-trough  pair  as  found  by  the  search 
technique  did  not  meet  the  acceptability  requirements  described  in  PART  III. 
The  number  of  periods  analyzed  before  the  error  occurred  is  printed  out  to 
show  how  much  good  data  were  considered  before  the  error  occurred.  As  a  note, 
a  search  error  usually  occurs  on  gages  with  extremely  low  amplitudes  (e.g., 
<0.05  prototype  feet),  or  if  data  were  acquired  before  a  model  had  properly 
set  up. 

33.  A  limited  output  option  is  identical  in  every  respect  with  the 
full  output  option  except  that  the  crests,  troughs,  and  wave  heights  are  not 
itemized . 

34.  The  summary  output  is  designed  to  provide  a  minimum  amount  of 
printout  but  still  provide  the  most  commonly  needed  statistics  from  the  anal¬ 
ysis.  This  form  of  output  fits  on  one  page  of  computer  paper.  The  test  iden¬ 
tification  information  in  Plate  6  is  located  at  the  top  of  the  output  page  for 
summary  output  (Plate  5).  Each  gage  has  one  line  of  output  containing  gage 
name,  average  wave  height  in  prototype  feet  with  standard  deviation,  average 
of  highest  one-third  wave  heights  in  prototype  feet  with  standard  deviation, 
calculated  model  and  prototype  periods  in  seconds  with  standard  deviations, 
average  of  highest  one-third  wave  heights  with  Keulegan's  friction  factors 
applied,  error  code,  and  number  of  data  cycles  or  periods  found  by  the  search 
technique . 

35.  The  columns  of  information  at  the  extreme  right  of  Plate  5  refer  to 
the  error  code  and  the  number  of  data  cycles  found.  First  the  gage  number  is 
listed  again,  then  an  error  code  of  0,  1,  or  2  is  listed.  An  error  code  of  0 
means  that  no  error  occurred.  An  error  code  of  1  means  that  the  search  rou¬ 
tine  found  fewer  cycles  of  data  than  it  was  searching  for.  This  is  the  same 
error  as  shown  in  Plate  7  by  "***  only  23  Periods  of  Data  Found  ***."  An 
error  code  of  2  means  that  a  "search  error"  was  encountered.  The  last  column 
refers  to  the  number  of  data  cycles  found  by  the  search  technique.  For  an 


error  code  of  0,  the  number  of  cycles  found  will  be  the  same  number  searched 
for.  For  an  error  code  of  1  the  number  of  cycles  should  be  no  less  than  three 
or  four  cycles  less  than  the  number  searched  for.  For  an  error  code  of  2  the 
number  of  cycles  found  will  be  unpredictable.  However,  the  statistics  listed 
will  correspond  only  to  those  wave  cycles  found. 

Magnetic  tape 

36.  As  an  option,  an  output  magnetic  tape  may  be  generated  by  this 
program.  Since  this  program  was  written  to  be  run  on  an  EAI  Pacer  100  mini¬ 
computer  with  only  two  magnetic  tape  units,  the  program  will  automatically 
consider  the  output  tape  to  be  the  tape  unit  which  is  not  declared  to  be  the 
input  tape  unit. 

37.  The  output  tape  is  simply  a  direct  copy  of  some  (not  all)  of  the 
input  tape  records  with  the  addition  of  another  record  containing  analysis 
results.  This  output  tape  is  desirable  for  later  data  comparison  in  the  form 
of  computer-generated  analysis  result  tables  or  plots. 

38.  There  is  no  differentiation  between  a  calibration  file  and  a  non¬ 
calibration  file  for  the  output  tape.  The  file  format  is  aa  follows: 

a.  File  Identification  Record  -  Same  as  input  tape. 

b.  First  File  Header  Record  -  Same  as  input  tape. 

c.  Second  File  Header  Record  -  Same  as  input  tape. 

d.  Transducer  Flag  Record  -  Same  as  input  tape. 

e.  Gage  Name  Record  -  Same  as  input  tape. 

f.  Weighted  Mean  Depth  Record  -  Same  as  input  tape. 

g.  Gage  Distance  Record  -  Same  as  input  tape. 

h.  Gage  Coordinate  Record  -  Same  as  input  tape. 

i.  Analysis  Results  Record  -  This  record  contains  the  average 
highest  one-third  of  the  wave  heights  for  each  gage's  data. 
There  is  one  entry  in  this  record  for  each  gage. 

Graphics 

39.  As  an  option,  this  program  is  capable  of  generating  a  plot  of  the 
scaled  prototype  data.  The  peaks  and  troughs  as  found  by  the  search  tech¬ 
nique  described  in  paragraph  10  are  labeled  on  the  plot  by  asterisks.  Only 
the  first  3000  data  points  may  be  plotted;  this  is  due  to  buffer  limitations 
within  the  program.  Plate  9  is  an  example  of  a  data  plot. 
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PART  III:  TIDAL  ANALYSIS  PROGRAM 


40.  The  tidal  analysis  program  can  process  three  types  of  data:  water- 
surface  tidal  elevations,  velocities,  and  temperatures.  The  latter  two  types 
of  data  (velocities  and  temperatures)  are  both  handled  in  the  same  manner. 

They  are  simply  scaled,  using  scaling  factors  which  are  input  from  magnetic 
tape,  and  plotted  with  respect  to  model  time. 

41.  Water-surface  elevation  data  are  scaled,  using  scaling  factors 
computed  from  informatiom  on  magnetic  tape,  and  also  plotted.  The  data  are 
fitted  to  a  series  of  cosine  functions  (harmonic  analyses)  by  a  Legendre 
method  of  least  squares  approximation  (Hamming  1950) ,  which  is  fully  described 
in  the  following  sections  of  PART  III.  Results  of  the  least  squares  analysis 
is  printed  in  tabular  form.  At  this  point,  the  analysis  can  take  one  of  two 
possible  courses.  The  model  data  can  be  compared  with  the  least  squares  fit 
or  it  can  be  compared  with  a  reference  data  set  (usually  prototype  data).  An 
overlay  plot  of  the  original  and  comparison  data  can  be  generated. 

42.  The  residual  of  the  comparison  and  original  data  may  be  plotted 
with  the  option  of  computing  a  power  spectral  density  on  the  residuals.  The 
power  spectral  density  estimate  is  computed  by  means  of  FFT  (Fast  Fourier 
Transform)  techniques  (Electronics  Associates  Inc.  1976).  The  procedure  is 
to  first  compute  the  direct  FFT  of  residual  data  and  then  compute  the  abso¬ 
lute  value  squared.  Thus,  if  the  result  of  the  FFT  is 


FFT(n)  =  f  f t  j (n)  +  fft2(n) 

then 

PSD(n)  =  [ f ft x (n) ] 2  +  [fft2(n)]2 

where 

FFT(n)  =  the  Fast  Fourier  Transform  results 
fftj(n)  =  real  part  of  FFT  results 
fft2(n)  =  imaginary  part  of  FFT  results 
PSD(n)  =  results  of  power  spectral  density 
The  result  of  the  power  spectral  density  can  be  plotted  as  an  option. 
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Program  Specif i cs 
Least  squares  harmonic  analysis 

43.  The  least  squares  fitting  procedure  mentioned  earlier  considers 
the  data  to  be  an  equally  spaced  series,  f(nAt)  ,  with  n  =  1 , . . . ,N  where 
N  is  the  total  number  of  data  samples  in  the  series  and  At  is  the  time 
interval  between  consecutive  samples.  The  series  is  assumed  to  be  periodic 
of  the  form 

J 

f(nAt)  =  f  (nAt)  +  e(nAt)  =  a  +  I  [a.  cos  (ui.nAt) 

0  i=l  1  1 


+  b.  sin  (uj.nAt)]  +  e(nAt) 


(1 


where 


f"(nAt)  =  series  less  noise  (true  signal) 
e(nAt)  =  noise  associated  with  transfer  function  of  instruments,  etc. 


aQ  =  mean  of  f(nAt) 


J  =  total  number  of  components  and  is  <  N/2 


w.  =  2n/T. 
t  i 


.th 


T^  =  wave  period  of  i  component 
44.  The  fitting  procedure  minimizes  the  error  produced  by  e(nAt)  by 
using  a  Legendre  method  of  least  squares  such  that  the  following  is  true: 


N  2  N  ,  2 

E  =  1  e  (nAt)  =  I  [f(nAt)  -  f  (nAt)]  ->  minimal 


(2 


n=l 


n=l 


where  E  =  error  that  the  least  squares  procedure  minimizes.  To  minimize  the 
error,  set 


9E  n  ,  9E  _  „  ,  .  _  .  T 

§T  "  0  and  9b-  ~  °  f0r  J  -  !»..., J 


Since 


IS-  =  gfr  {  I  f  '(nAt)  ] 2  -  2f  (nAt )}  x 


N 


1  [a.  cos  (ui.nAt) 

i  =  l  1  1 


+  b  .  sin  (u> .  nAt)  ] }  + 


I  [a.  cos  (ui.nAt)  +  b.  sin  (lu.nAt)] 


for  j  =  1 , . . . , J  and 


§77  = 


ft) 


then 


9E  N  r  J 

~ —  =  I  r-2f'(nAt)  cos  (ui.nAt)  +  2  cos  (uj.nAt)  I  [a.  cos  (uj.nAt) 

o3  -iL  1  i  .  -  j  i 

j  n=l  J  i=l  J 


+  sin  (u^nAt)]  for  j  =  1, 


Similarly, 


aF  n  r  J 

=  I  -2f'(nAt)  sin  (uj.nAt)  2  sin  (uj.nAt)  I  [a.  cos  (uj.nAt) 
db.  j  j  J  .  ,  l  l 

j  n=l  J  J  i=l 


+  b.  sin  (w.nAt)]  for  j  =  1, 
ii 


By  setting  9E^ 
9a. 


0  and  3E 

9b. 

J 


=  0  for  j  =  1,...,J  the  equations 


reduce  to 


N 
I 
n=  1 


cos  (uKnAt)  I  ( a ^  cos  (urnAt)  +  b^  sin  (u^nAt)] 


i=l 


N 


=  I  [f'(nAt)  cos  (uj.nAt)]  for  j 


n=l 


and 


N 

I 

n=  1 


J 

sin  (ui.nAt)  I  [a.  cos  (uj.nAt)  +  b.  sin  (uKnAt)] 

J  i=l 


=  I  (f'(nAt)  sin  (uj.nAt)]  for  j 


H 


which  is  a  set  of  2J  simultaneous  equations  for  a  and  h. 

45.  In  matrix  notation,  these  equations  can  be  written  in  the  follow 
ing  form: 

ICJ  [ A j  =  [F] 


C1  C1 
C2  C1 


ClC2---CiCj  C!S1 

C2  C2-'-C2  Cj  C2  S1 


C1  S2-‘-Cl  Sj 

C2  V--C2  Sj 


CJ  C1 


S1  C1 


S2  C1 


C.  C  c„  c.  C.  s, 

J  2  3  j  j  1 


S1  S'  -S,  ^ 


S2C2---S2Cj  S2S1 


C.  S...C.S.  X  a.  =  c  f 
J  2  J  J  J  j 


si  V--s,  sj 


S2  s2...s2  s. 


LSjC‘  Vr-'Vj  V.  5 :  s.' "  * s :  sj  J  [»,  *_ 

where  C  is  a  symmetric  matrix  with 


C  =  I  cos  (w  nAt)  j  =  1,2,. .  .  ,J 
J  n=l  J 

N 

S  =  I  sin  (w  nAt)  j  =  1,2,. ...J 
J  n=l  J 

N 

Cif  =  2  (f'(nAt)  cos  (w.nAt) ]  j  =  1,2,...,J 

J  n=l  J 


S  f  =  I  [  f  '(nAt)  sin  (uj.nAt)]  i  =  1,2 . J 


The  method  of  solution  for  a 


and  bj  for  j  =  1 ,  .  .  .  ,  J  is: 


[A]  =  [cf1  [f] 


equal  to 


46.  It  can  be  shown  that  I  a.  cos  (w.nAt)  +  b.  sin  (w.nAT)  is 

i  i  i  i 


I  Amp.  x  cos  lui.(nAt)  -  PHS .  ] 
i=l  1  1  1 


where 


2  2 

Amp^  =  a^  +  b^  ,  commonly  called  amplitude 


PHS.  =  ARCTAN 
i 


,  commonly  called  phase  angle 


The  latter  equation  is  easier  to  work  with  and  is  used  in  this  program  where 
results  of  the  lease  squares  harmonic  analysis  are  expressed  as  Amp.  and 
PHS.  . 

l 

Scaling  data 

47.  As  mentioned  previously,  both  velocity  and  temperature  data  are 
scaled  from  scaling  factors  read  directly  from  magnetic  tape.  Tidal  data  are 
handled  much  differently.  There  is  only  one  transducer  for  all  tidal  gages, 
with  a  multiplexer  selecting  which  gage  the  transducer  is  acting  upon;  thus 
the  same  scaling  factors  can  be  applied  to  all  tidal  gages. 

48.  With  every  data  scan  during  data  acquisition,  reference  tidal  gages 
also  are  scanned  (usually  three).  These  reference  gages  are  located  in  a 
still-water  pool  that  does  not  oscillate  with  the  model.  The  sensors  of  ea_h 
reference  gage  are  extended  a  known  displacement  into  the  water.  This  "known" 
displacement  is  recorded  in  the  gage  coordinate  record  in  the  "depth"  location 
for  each  reference  gage  (Appendix  F).  By  relating  this  known  displacement  to 
the  reading  obtained  from  each  reference  gage  during  a  scan,  it  is  possible  to 
compute  new  calibration  coefficients  for  the  tidal  transducer  with  each  data 
scan.  The  purpose  for  computing  new  update  coefficients  with  each  scan  is  to 
counteract  any  electronic  drift  within  the  sensor. 

49.  These  tidal  calibration  coefficients  are  computed  by  performing  a 
linear  least  squares  fit  on  the  reference  channels. 

Reference  elevations 


50.  Each  tidal  gage  has  a  zero  reference  reading  recorded  on  magnetic 


tape.  This  zero  reference  reading  is  taken  at  a  known  elevation  on  the  model 
with  the  model  water  surface  still.  This  zero  reference  reading  is  consid¬ 
ered  to  be  the  known  elevation  of  the  model.  Thus  each  virtual  reading  is 
computed  in  the  following  manner: 

VR  =  (AR  -  ZRR)SFS  +  KE  +  SFC 


where 

VR  =  virtual  reading 
AR  =  actual  reading 
ZRR  =  zero  reference  reading 
SFS  =  scale  factor  slope 

KE  =  prototype  water  elevation  at  time  zero  reference  readings  made 
SFC  =  scale  factor  constant 

51.  Reference  gages  are  read  prior  to  being  moved  to  their  known  dis¬ 
placements.  That  is  to  say  that  all  reference  gages  are  at  the  same  displace 
ment  when  they  are  read  for  their  zero  reference  readings. 

Phase  adjustment 

52.  As  described  previously,  results  of  the  least  squares  harmonic 
analysis  are  expressed  as  amplitude  and  phase.  A  data  series  can  be  recon¬ 
structed  from  this  analysis  by  using  the  following  equation: 

J 

f'(nAt)  =  a  +  2  Amp.  cos  [u).(nAt  -  PHS.)]  (3 

i=l  1 

where 

f'(nAt)  =  reconstructed  data  series 

a  =  mean  of  data  series 
o 

Amp^  =  amplitude  of  component  i 
uk  =  2 n  period  of  component  i 
PHS^  =  phase  angle  of  component  i 

53.  It  is  a  program  option  to  compare  this  reconstructed  data  series 
with  the  original  data  on  tape.  Another  option,  as  previously  mentioned,  is 
to  compare  the  original  data  with  a  reference  data  set  such  as  prototype  data 
data  from  another  test,  etc.  However,  it  is  not  possible  to  compare  the 
original  data  with  both  a  reference  data  set  and  the  analysis  in  the  same 
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execution  run  of  the  program.  The  program  must  be  executed  twice  to  exercise 
both  options. 

54.  In  comparing  the  original  data  with  reference  data,  components  from 

previous  analysis  are  input.  Phases,  amplitudes,  periods,  and  mean  water 
level  are  input  to  Equation  3.  Thus  a  data  series  is  generated  for  comparison. 
It  is  unlikely  that  values  of  phase  angles  for  the  comparison  data  will  agree 
with  that  of  the  original  data;  therefore  the  phase  for  the  two  data  sets 
should  be  adjusted.  Whenever  original  data  are  to  be  compared  with  reference 
data,  the  program  considers  the  first  gage  that  it  analyzes  to  be  a  deepwater 
reference  gage.  The  deepwater  reference  gage  should  most  closely  represent 
the  tide  that  is  being  applied  to  the  estuary.  Thus  the  deepwater  gage  should 
be  remote  enough  from  the  estuary  such  that  modification  or  changes  in  the 
estuary  would  not  affect  data  on  the  deepwater  reference  gage.  Therefore  the 
amplitudes  should  remain  constant  from  data  series  to  data  series  on  the  deep¬ 
water  gage  for  a  particular  tide.  The  phase  angles  should  also  remain  con¬ 
stant  relative  to  one  another.  Thus,  with  the  deepwater  gage  remaining  con¬ 
stant  from  data  series  to  data  series,  it  is  possible  to  adjust  phases  so  that 
they  will  closely  relate  to  this  deepwater  reference  gage.  This  is  done  by 
subtracting  the  phase  of  the  component  of  the  reference  data  from  the 

phase  of  the  component  of  the  original  data  and  adding  this  difference 

to  all  phases  of  all  components  of  the  reference  data  for  all  gages.  This 
effectively  adjusts  the  phases  for  the  entire  run  if  the  phases  for  all  gages 
on  the  reference  data  set  are  correlated.  It  is  only  necessary  to  take  the 
difference  of  the  M2  component  phase  since  for  the  deepwater  reference  gage 
the  phase  gage  angles  for  each  component  should  remain  constant  relative  to 
one  another.  It  also  should  be  noted  that  the  M2  component  is  the  most  in¬ 
fluential  component  (that  is,  the  component  with  the  largest  amplitude)  and 
that  component  should  have  the  least  error  in  phase. 

In-line  assembly  language 

55.  The  EAI  Pacer  100  FORTRAN  compiler  allows  the  programmer  to  use 
some  EAI  Pacer  100  assembly  language  instructions  in  line  with  FORTRAN  coding. 
The  only  three  assembly  language  instructions  used  in  this  program  are  the 
Jump,  Link,  and  Octal  instructions.  These  instructions  have  been  explained 

in  detail  in  paragraphs  16-20  of  PART  II  and  will  not  be  repeated  here. 


Program  Inputs 


Hi 


Teletype  inputs 

56.  Several  program  initializations  are  input  by  teletype.  The  func¬ 
tions  and  pertinent  related  information  of  each  of  these  inputs  are  defined 
in  this  section. 

a.  Force  Shift  -  If  reference  data  are  to  be  compared  with  the 
model  data,  a  shift  in  phase  is  performed  on  the  reference  data 
so  that  it  will  more  closely  relate  to  the  model  data.  This 
automatic  shift  may  or  may  not  be  desirable  so  an  overriding 
force  shift  may  be  performed  to  further  enhance  or  counteract 
the  automatic  shift.  This  shift  is  input  in  degrees.  Thus  a 
value  of  0  for  the  force  shift  causes  no  further  shift  action 
to  be  taken  or  a  value  of  180  would  cause  a  shift  of  180  deg 
for  the  reference  data. 

b.  Taper  Residual  -  If  it  is  desirable  to  apply  a  taper  cosine 
bell  function  to  the  residual  of  the  model  data  and  the  com¬ 
parison  data  before  power  spectrum  is  performed,  this  would  be 
input  as  "Y."  If  not,  an  "N"  would  be  input. 

c.  Calibration  Fit  Option  -  As  mentioned  previously,  a  least 
squares  linear  fit  is  performed  on  the  reference  gages  (usually 
three)  to  determine  scale  factors  for  the  tidal  gages.  This 
option  allows  the  fit  to  be  performed  on  any  two  of  the  refer¬ 
ence  gages  or  all  three  of  them.  A  value  of  1  would  cause  the 
least  squares  fit  to  be  performed  on  the  reference  gage  with  a 
positive  displacement  and  the  reference  gage  at  zero.  A  value 
of  2  would  cause  the  least  squares  fit  to  be  performed  on  the 
reference  gage  with  a  negative  displacement  and  the  reference 
gage  at  zero.  A  value  of  3  would  cause  the  least  squares  fit 
to  be  performed  on  the  reference  gage  with  a  positive  displace¬ 
ment  and  the  reference  gage  with  the  negative  displacement.  A 
value  of  4  would  cause  the  least  squares  fit  to  be  performed  on 
all  three  of  the  previously  mentioned  reference  gages. 

d.  Number  of  Components  for  the  Analysis  -  The  number  of  least 
squares  harmonic  analysis  components  is  input. 

e.  Component  Period  -  A  component  period  in  prototype  hours  is 
input  for  each  component. 

f.  List  Update  Calibrations  -  If  this  is  answered  by  "Y,"  the  up¬ 
date  calibration  coefficients  for  the  tidal  gages  are  listed  as 
they  are  computed.  If  "N"  is  input,  the  coefficients  are  not 
listed. 

g.  Reference  Data  to  be  Input  -  If  this  is  answered  by  "Y,"  refer¬ 
ence  components  will  be  read  from  cards  as  each  data  channel  is 
reduced.  These  reference  components  will  be  used  to  generate 
reference  data  to  compare  with  the  particular  gage  data. 

h.  Number  of  Reference  Components  to  be  Input  -  If  reference  data 
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are  to  be  input,  the  number  of  components  to  be  read  from  cards 
is  input  here. 

i.  How  Often  Update  Calibration  -  The  number  of  data  scans  to  skip 
between  computing  update  calibration  coefficients  is  input.  If 
the  value  is  1,  update  calibration  coefficients  will  be  com¬ 
puted  for  every  scan.  It  follows  that  a  value  of  3  would  cause 
new  coefficients  to  be  computed  for  every  third  data  scan. 

2-  Label  -  If  a  label  for  the  tabular  output  and  plots  is  desired, 
it  is  input.  A  maximum  of  20  characters  may  be  input. 

k.  Number  of  Records  to  Skip  -  If  it  is  desired  to  skip  some  data 
before  the  analysis  is  done,  the  number  of  data  records  (on 
tape)  to  be  skipped  is  input. 

l.  Tape  Unit  -  The  magnetic  tape  unit  on  which  the  data  tape  is 
located  is  input  here.  It  will  be  either  14  or  15  for  an  EAI 
Pacer  100  with  two  tape  drives. 

m.  Average  Model  Depth  -  The  average  model  depth  is  input  as 
model  feet.  This  is  used  in  scaling  a  distorted  model. 

n.  File  Name  -  The  data  file  name  is  input  as  ASCII  characters. 

o.  Number  of  Data  Cycles  to  Analyze  -  The  number  of  cycles  of 
data  to  analyze  for  this  program  execution  is  input. 

Sense  switch  settings 

57.  Sense  switches  are  single  throw  switches  located  on  the  control 
panel  of  the  Pacer  100.  They  are  manually  set  and  can  be  read  by  the  CPU. 

The  FORTRAN  statement,  IF  (SENSW  (1))  GO  TO  10,  causes  program  execution  to 
branch  to  statement  number  10  if  sense  switch  A  or  switch  (1)  is  set  or  is 
.TRUE.  It  follows  that  sense  switches  are  treated  as  logical  variables.  The 
advantage  to  using  sense  switches  is  that  they  may  be  set  or  reset  while  a 
program  is  executing  without  impeding  program  execution.  The  functions  of  the 
sense  switches  in  this  program  are  as  follows: 

a.  Sense  Switch  A  [SENSW  (1)]  -  If  sense  switch  A  is  reset,  a  plot 
of  the  original  data  is  generated  for  each  channel.  If  set, 
the  plot  is  not  generated. 

b.  Sense  Switch  B  [SENSW  (2)]  -  If  sense  switch  B  is  reset,  an 
overlay  plot  of  the  original  data  and  comparison  data  is  gen¬ 
erated.  The  mean  of  the  comparison  data  will  be  included  in 
the  comparison  data  and  the  mean  of  the  original  data  will  be 
included  in  the  original  data. 

c.  Sense  Switch  C  [SENSW  (3)]  -  If  sense  switch  C  is  reset,  an 
overlay  plot  of  the  original  data  with  its  mean  removed  and  the 
comparison  data  with  its  mean  removed  is  generated. 

d.  Sense  Switch  D  [SENSW  (4)]  -  If  sense  switch  D  is  reset,  a  plot 
of  the  residual  of  the  original  data  and  the  comparison  data  is 
generated . 
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e.  Sense  Switch  E  (SENSW  (5)]  -  If  sense  switch  E  is  reset,  the 
power  spectral  density  option  is  performed  on  the  residual  with 
a  plot  of  the  results  generated. 

f.  Sense  Switch  F  [SENSW  (6)]  -  The  program  will  skip  data  chan¬ 
nels  as  long  as  sense  switch  F  is  set.  Once  the  switch  is 
reset,  the  program  will  continue  to  process  data  channels  in 
sequential  order. 

g.  Sense  Switch  G  [SENSW  (7)]  -  If  sense  switch  G  is  reset,  the 
original  scaled  data  will  be  listed  in  tabular  form. 


Magnetic  tape 

58.  Testing  parameters,  gage  identification,  scaling  factors,  test 
identification,  and  test  data  are  all  recorded  or  either  derived  from  informa¬ 
tion  recorded  on  magnetic  tape.  A  general  tape  format  is  described  in  this 
section  with  more  specific  information  available  in  Appendix  F.  The  test  data 
are  recorded  in  files  on  magnetic  tape.  There  are  two  general  file  formats 
containing  a  variable  number  of  records. 

59.  The  two  file  formats  can  be  called  a  calibration  file  and  a  non¬ 
calibration  file.  In  a  calibration  file,  several  records  are  devoted  to  gage 
identification  and  data  scaling  coefficients.  Since  this  information  remains 
constant  for  a  series  of  data  tests,  there  is  no  need  to  record  it  in  every 
data  file.  A  calibration  file  and  a  noncalibration  file  are  similar  except 
that  some  records  included  in  the  calibration  file  are  excluded  from  the  non¬ 
calibration  file.  Thus  this  exclusive  information  recorded  in  a  calibration 
file  is  applicable  to  all  proceeding  data  files  until  updated  by  another  cali¬ 
bration  file. 

60.  Each  file  is  separated  by  an  end-of-file  (EOF)  mark  and  the  file 
formats  are  as  follows: 

a.  File  Identification  Record  -  This  record  is  five  elements  long. 
The  first  three  elements  contain  a  six-character  ASCII  code 
denoting  the  file  name.  The  fifth  element  denotes  the  length 
of  the  next  record. 

b.  First  File  Header  Record  -  The  fifth  element  of  the  file  iden¬ 
tification  record  contains  the  elemental  length  (usually  20) 
of  this  record.  This  record  contains  testing  and  model  iden¬ 
tification  information.  The  20th  element  of  this  record  con¬ 
tains  the  length  of  the  next  record. 

c.  Second  File  Header  Record  -  The  20th  element  of  the  first  file 
header  record  contains  the  elemental  length  of  this  record. 

This  record  contains  testing  parameters  and  model  control 
information.  The  model  control  information  is  not  used  in 
reducing  the  data. 


d.  Transducer  Flag  Record  -  This  record  contains  two  flags  for 
each  sensor  sampled  during  data  acquisition.  These  flags  con¬ 
tain  data  scaling  information  revealing  the  length  of  the 
calibration  coefficient  record  and  calibration  voltage  record. 
This  record  exists  only  in  a  calibration  file. 

e.  Calibration  Coefficient  Record  -  This  record  length  is  deter¬ 
mined  by  the  number  of  gages  having  calibration  coefficients 
associated  with  them.  A  flag  in  the  transducer  flag  record 
determines  whether  a  gage  has  coefficients  in  the  coefficient 
record.  There  are  11  coefficients  for  each  gage  that  has  co¬ 
efficients  recorded  in  this  record.  If  no  gages  have  entries 
in  this  record,  the  record  is  not  recorded  on  tape.  This 
record  exists  only  in  a  calibration  file. 

f.  Gage  Name  Record  -  This  record  contains  two  elements  (4  ASCII 
characters)  for  each  gage.  These  ASCII  characters  are  used 
for  gage  labeling  and  identification.  This  record  exists  only 
in  a  calibration  file. 

g.  Gage  Coordinate  Record  -  This  record  contains  three  floating 
point  coordinates  for  each  gage:  an  X  and  Y  coordinate  deter¬ 
mining  location  and  Z  coordinate  determining  water  depth  for 
each  gage  that  is  recorded.  This  record  exists  only  in  a 
calibration  file. 

h.  Reference  Potentiometer  Calibration  Coefficient  Record  -  This 
record  is  not  used  by  this  program  and  is  skipped  during  pro¬ 
gram  execution.  This  record  exists  only  in  a  calibration  file. 

i.  Calibration  Voltage  Record  -  This  record  is  not  used  by  this 
program  and  is  skipped  during  program  execution.  This  record 
exists  only  in  a  calibration  file. 

Zero  Reference  Reading  Record  -  This  record  contains  the  zero 
reference  readings  for  each  tidal  gage  and  each  tidal  reference 
gage.  This  record  exists  only  in  a  calibration  file. 

k.  Data  Records  -  The  remaining  records  are  data  records.  The 
length  of  these  records  is  equal  to  the  number  of  the  data 
gages  times  the  number  of  data  samples  per  tidal  cycle  divided 
by  the  number  of  data  records  pe-r  tidal  cycle.  The  number  of 
data  records  is  equal  to  the  number  of  tidal  cycles  sampled 
times  the  number  of  records  per  tidal  cycle. 


Cards 

61.  As  mentioned  previously,  if  original  data  are  to  be  compared  with 
reference  data,  the  components  of  the  reference  data  are  input  on  cards.  This 
is  the  only  information  that  is  input  on  cards  and  if  reference  data  are  not 
to  be  compared,  no  cards  are  needed  for  program  execution.  The  necessary  in¬ 
formation  to  be  input  from  cards  is  a  period,  phase,  and  amplitude  for  each 
component  of  the  reference  data  as  well  as  a  mean  of  the  data.  There  will  be 
one  card  for  each  component  in  the  following  FORTRAN  format,  3F10.5.  The 


the  first  "F"  field  will  be  for  the  period  in  prototype  he.  irs,  the  second  "F" 
field  will  be  for  phase  angle  in  degrees,  and  the  third  "F"  field  will  be  for 
amplitude  in  prototype  feet.  The  last  card  will  contain  the  average  depth  or 
"mean"  in  model  feet.  Thus,  if  the  reference  data  were  to  be  generated  for 
three  components  having  the  following  periods,  phase,  and  amplitudes: 


Component 

No. 


be  as  follows: 


Card  #1 
Column  1 


Period 

hr 

Phase 

deg 

Amplitude 

ft 

10 

45 

5 

5 

50 

2 

2.5 

90 

0.5 

3.52  ft, 

the  cards 

for  that  particular 

1 

11 

21 

10.0 

45.0 

5.0 

5.0 

50.0 

2.0 

2.5 

90.0 

0.5 

3.52 

Program 

Outputs 

Teletype 

62.  All  teletype  outputs  generated  by  this  program  are  to  aid  in  pro¬ 
gram  operations.  The  teletype  inputs  described  in  paragraph  56  are  input 
upon  request  from  teletype  outputs;  that  is,  each  teletype  input  is  requested 
by  a  teletype  output  in  a  question-answer  fashion. 

63.  The  only  other  teletype  outputs  generated  by  this  program  are  gage 
numbers.  As  each  gage  is  processed,  the  gages ' s  number  is  output  to  the  tele¬ 
type.  This  is  done  so  that  the  operator  can  monitor  program  progression  and 
alter  sense  switch  settings  for  different  gages  if  desired. 

64.  Examples  of  teletype  output  as  well  as  teletype  inputs  are  illus¬ 
trated  in  Plates  10  and  11.  Notice  that  all  teletype  outputs  except  the  gage 
number  listings  are  followed  by  question  marks.  All  other  listings  represent 
teletype  inputs.  Notice  also  that  Plate  10  depicts  the  analysis  of  three  com¬ 
ponents  and  Plate  11,  only  one  component.  The  reference  data  consisting  of 


two  components  are  input  for  Plate  11  whereas  reference  data  are  not  input 
for  Plate  10. 

Printer 


65.  The  first  page  of  printer  output  is  test  identification  information 
This  information  comes  from  the  first  two  header  records  on  magnetic  tape. 

66.  The  test  identification  information  is  followed  by  a  listing  of 
update  calibration  coefficients  if  such  a  listing  is  requested.  Plate  12 
illustrates  the  test  identification  output  including  the  list  option  for  cali¬ 
bration  coefficients. 

67.  For  each  tide  gage  analyzed,  a  page  of  output  is  generated  similar 
to  Plates  13-15.  The  difference  between  Plates  13  and  14  is  that  sense  switch 
G  was  .TRUE,  when  Plate  14  was  generated  causing  the  two  cycles  of  scaled 
prototype  data  to  be  analyzed  to  be  listed.  In  Plate  15,  a  taper  function  was 
applied  to  the  residual  data  prior  to  performing  the  power  spectral  density; 
thus  the  statistics  of  the  residual  data  are  listed.  Also,  gage  and  test 
identifications  are  printed  on  each  analysis  page  so  if  a  page  is  misplaced 

it  can  still  be  identified.  The  mean  and  variance  of  the  original  and  re¬ 
sidual  data  as  well  as  results  of  the  least  squares  harmonic  analysis  are 
listed.  Plate  13  illustrates  a  three-component  analysis,  Plate  14  two  compo¬ 
nents,  and  Plate  15  four  components. 

68.  Plate  15  is  an  example  of  an  output  page  when  reference  data  are 
compared  with  the  original  data.  The  difference  between  Plate  13  and  Plate  15 
is  that  the  shifts,  mentioned  in  paragraphs  50-54,  are  listed  in  Plate  15. 
Graphics 

69.  Plate  16  is  an  example  of  a  one  cycle  original  data  plot. 

Plates  17-19  all  are  plots  of  the  original  data  (solid  line)  and  comparison 
data  (dotted  line) .  Plates  17  and  18  have  the  means  of  the  original  data 
removed  from  both  the  original  data  and  the  comparison  data.  Plates  18  and 

19  have  comparison  data  input  from  cards  while  Plate  17  uses  the  least  squares 
harmonic  analysis  results  for  comparison  data.  Plate  19  has  the  mean  of  the 
original  data  included  and  the  comparison  data  (dotted  line)  is  shifted  by 
90  deg.  Plate  17  depicts  only  one  cycle  of  data  while  Plates  18  and  19  depict 
two  cycles  of  data. 

70.  Plate  20  is  an  example  of  a  residual  plot,  which  is  simply  the 
comparison  data  subtracted  from  the  original  data.  This  residual  plot  is  the 
residual  of  the  two  data  sets  in  Plate  17. 


71.  Plates  21  and  22  are  examples  of  the  power  spectral  density  result 
Both  plates  are  plots  of  the  data  sets  in  Plate  20  with  Plate  22  illustrating 
the  results  of  the  taper  function  performed  on  the  residual  prior  to  the  powe 
spectral  density. 

72.  Plate  23  is  an  example  of  velocity  data.  This  plate  is  a  plot  of 
the  scaled  original  velocity  data.  The  original  data  plot  is  the  only  output 
that  this  program  will  generate  for  velocity  data  or  temperature  data. 

Plate  24  is  an  example  of  a  temperature  plot. 
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***  WAVE  ANALYSIS  BASED  ON  SIPPLE  SEARCH  TECHNIQUES  *** 
DAVE  DATA 

CHANNEL  8  GAGE  ARR.  0 


XY  COORDINATES  .0000 

UATER  DEPTH  IS  60.0000 
GAGE  ♦  IS  RD06 

CRESTS 

VALUE  POSITION 


3.198 
3.131 
3.096 
3.107 
3.091 
3.068 
3.051 
3.062 
3.068 
3.102 
3. 107 
3.085 
3.080 
3.074 
3.080 
3.085 
3.068 
3.051 
3.096 
3.091 
3.091 
3.119 
3.107 
3.125 


UAVEHEIGHTS 

AVERAGE 

RMS 


.000 
5.876 
11.852 
17.728 
23.703 
29.679 
35.654 
41.630 
47.605 
53.481 
59.457 
65.433 
71.408 
77.384 
83.260 
89.235 
95.211 
101.186 
107.062 
113.038 
119.014 
124,889 
130.865 
136.841 


.0000 


TROUGHS 

VALUE  POSITION 

2.888 


-3.232 
-3.117 
-3.122 
-3. 192 
-3. 107 
-3.056 
-3. 102 
-3. 122 
-3. 122 
-3.142 
-3.112 
-3. 122 
-3.112 
-3. 107 
-3. 138 
-3.117 
-3.092 
-3.067 
-3.097 
-3.082 
-3.112 
-3.092 
-3.082 
-3.117 


(JAVEHEIGHT 
VALUE  POSITION 


8.864 
14.839 
20.715 
26.790 
32.766 
38.642 
44.618 
50.593 
56.569 
62.445 
68 . 420 
74.296 
80.272 
86.247 
92.223 
98.099 
104.075 
110.050 
115.926 
121.902 
127.877 
133.753 
139.729 


MODEL  -  INCHES 


33.33  PCT 


UAVEPER10D 


.993  +✓-  .010 

.993 

1.002  +✓-  .011 

MODEL  -  SECONDS 


.685  +/- 


.042 


6.430 
6.248 
6.218 
6.300 
6.  198 
6. 125 
6,153 
6.185 
6. 190 
6.243 
6.219 
6.208 
6.192 
6.181 
6.217 
6.203 
6.161 
6.118 
6. 192 
6.173 
6.203 
6.212 
6.190 
6.242 


1.494 
7.370 
13.345 
19.221 
25.297 
31.272 
37.148 
43. 124 
49.099 
55.075 
60.951 
66.926 
72.902 
78.878 
84.754 
90.729 
96.705 
102.680 
108.556 
114.532 
120.507 
126.383 
132.359 
138.335 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
2! 
22 

23 

24 


PROTOTYPE  -  FEET 

6.208  +4-  .060 

6.209 

6.265  +✓-  .067 

PROTOTYPE  -SECONDS. 


5.928  +✓- 


.362 


PLATE  3 


>1 


»< 

Si 


ivL-.ls'l-.'L-. 


***  WAVE  ANALYSIS  BASED  ON  SIMPLE  SEARCH  TECHNIQUES  *** 


WAVE  DATA 

CHANNEL  12 

GAGE  ARR. 

0 

XY  COORDINATES 

.0000 

.0000 

WATER  DEPTH  IS 

60.0000 

GAGE  #  IS  RDI2 

MODEL  - 

INCHES 

PROTOTYPE  -  FEE 

UAVEHEIGHTS 

AVERAGE 

1.077  +/- 

.011 

6.734  +✓- 

RMS 

1 

.077 

6.734 

33.33  PCT 

1.088  +✓- 

.002 

6.802  +/•- 

MODEL  - 

SECONDS 

PROTOTYPE  -SECONDS 

UAVEPERIOD 

. 685  +/- 

.042 

5.929  +✓- 
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PLATE  6 


URVt  MtI  ftNftLYSIS  8flSED  DN  SIMPLE  SEARCH  TECHNIQUES  *** 
CHANNEL  10  GAGE  ARR.  0 

XY  COORDINATES  .0000  ,0000 

WATER  DEPTH  IS  60.0000 

GAGE  •  IS  RDI0 

***  ONLY  23  PERIODS  OF  DATA  FOUND  *** 


UAVEHEIGHTS 

AVERAGE 

RMS 


MODEL  -  INCHES 

1.026  +✓-  .018 
1.029 


33.33  PCT  1.048  +/-  .007 

MODEL  -  SECONDS 


UAVEPERIOD 


.684  +/- 


PROTOTYPE  -  FEET 

6.428  +/-  .113 

6.429 

6.552  +✓-  .046 

PROTOTYPE  -SECONDS. 


5.927  +/'- 


f  '.v.  «*.y«y ^ v.* . .  i'.v!v.' . 


,*  V.V.  •  .  ■*. 


PLATE  7 


***  DAVE  ANALYSIS  BASED  ON  SIMPLE  SEARCH  TECHNIQUES  *** 
MODEL  WAVE  DATA 

CHANNEL  1  GAGE  ARR.  0 

XY  COORDINATES  .0000  .0000 

UATER  DEPTH  IS  60.0000 

GAGE  #  IS  RD01 


CRESTS  TROUGHS  UAVEHEIGHT 


VALUE 

POSITION 

VALUE 

POSITION 

VALUE 

POSITION 

2.624 

1.021 

2.627 

2.449 

-.003 

1.755 

1 

2.624 

2.694 

2.628 

4.776 

-.003 

3.756 

2 

2.624 

5.225 

2.627 

5.797 

-.002 

5.511 

3 

2.624 

6.083 

2.628 

8.981 

-.003 

7.552 

4 

SEARCH  ERROR- 

-PERIOD 

5 

2.626 

9.430 

2.626 

9.349 

.000 

9.389 

5 

2.624 

10.002 

2.627 

10.859 

-.002 

10.451 

6 

2.626 

11.471 

2.627 

11.635 

-.000 

11.553 

7 

2.624 

12.451 

2.626 

12.492 

-.002 

12.492 

8 

2.589 

15.595 

2.914 

17.595 

-.326 

16.615 

9 

2.570 

19.881 

2.929 

20.208 

-.359 

20.044 

10 

2.593 

22.371 

2.926 

22.657 

-.333 

22.535 

11 

2.573 

24.453 

2.873 

24.984 

-.300 

24.739 

12 

2.536 

27.025 

2.877 

27.515 

-.341 

27.270 

13 

2.578 

29.434 

2.916 

29.883 

-.339 

29.679 

14 

2.551 

32.047 

2.875 

32.373 

-.324 

32.210 

15 

2.555 

34.496 

2.883 

34.904 

-.327 

34.700 

16 

2.568 

36.823 

2.924 

37.272 

-.356 

37.068 

17 

2.615 

39.272 

2.953 

39.721 

-.338 

39.517 

18 

2.601 

41.722 

2.941 

42.252 

-.340 

42.008 

19 

2.596 

44.294 

2.939 

44.865 

-.343 

44.579 

20 

2.595 

46.947 

2.919 

47.396 

-.324 

47. 192 

21 

2.574 

49.274 

2.914 

50.009 

-.340 

49.642 

22 

2.558 

52.050 

2.913 

52.581 

-.356 

52.336 

23 

2.565 

54.622 

2.904 

54.989 

-.339 

54.826 

24 

2.572 

56.949 

2.904 

57.480 

-.332 

57.235 

25 

2.553 

59.358 

2.887 

59.970 

-.335 

59.684 

26 

2.570 

62.093 

2.860 

62.460 

-.289 

62.297 

27 

2.565 

64.583 

2.832 

64.991 

-.268 

64.787 

28 

2.535 

67.114 

2.885 

67.441 

-.349 

67.277 

29 

2.552 

69.523 

2.895 

69.849 

-.342 

69.686 

30 

MODEL 

-  INCHES 

PROTOTYPE 

-  FEET 

UAVEHEIGHTS 

AVERAGE 

-1.463 

+/- 

.882 

-.244  +✓- 

147 

RMS 

1.709 

.285 

33.33  PCT 

-.344 

+/- 

.664 

-.057  +/- 

111 

MODEL 

-  SECONDS 

PROTOTYPE  -SECONDS. 

UAVEPERIOD 

1.691  +/-  .252  2.392  +/-  .357 


PLATE  8 


MODEL  WOVE  D ATP 


a  50  00  60 .00  70.00 

TIME  PROTOTYPE  SECONDS 


FORCE  SHIFT::? 

90. 

TAPER  RESIDUAL  ? 

NO 

CALIBRATION  FIT  OPT1  0  N:  ? 

A 

#COttPS  TO  ANALIZE  FOP? 

3 

COMP  PER:? 

12.4226 

COMP  PEP:? 

6.2123 

COMP  PER:? 

A.  1*4  P'2 

LIST  UPDATE  CALS? 

YFS 

REFERENCE  DATA  TO  BE  INPUT? 
NO 

HON  OFTEN  UPDATE  CAL  ? 

1 

LAPEL- ? 

MODEL  TIDE  DATA 
*P EC  TO  SKIP:? 

n 

TAPE  UNIT:? 

1A 

AVEP4RE  MODEL  DEPTH:? 

.2 

FILE  N4ME7 
OPI  0P3 

»  OF  CYCLES  TO  ANALYSE:? 

2 

GAOF  *:  I 

RARE*:  2 

CASE*:  3 

RACE*:  A 

GAOE  *=  5 

RARE*:  6 

RARE*:  7 

RARE*:  p 


PLATE  10 


FORCE  SHI  rT:? 

9f. 


TAPE?  RESIDUAL  ? 

YES 

CALIBRATION  FIT  OPTION? 

A 

i»COM?S  TO  A NALIZ  F  FOR? 

1 

COMP  P F.R  -? 

12. APIS 

LIST  'JPDATF  CALS? 

MO 

REFERENCE  OAT  A  TO  RE  INPUT? 

YFS 

NUMBER  of  REFEPEUCE  COMPONENTS  :? 
2 

HOW  OFTEN  'JPOATF  CAL  ? 

2 

LABELS 

MODEL  TIOE  nATA 

/»PRC  TO  SKIP:? 

2 

TAPE  UNIT:? 

15 

AVERAGE  "OREL  DEPTH:? 

FILE  NA-'E? 

OP  I  0',.T> 

#  OF  CVCLES  TO  ANALYSE:? 


I 

GAGE#:  1 

GAGE*;  2 

GAGE*:  5 

GAGE#:  A 

GAGE*:  A 

GAGE*:  A 

G  A  r  E  *  :  7 

GAGE*;  c 

n 


GAGE*; 


***  TIDE  RECORD  ANALYSIS  PROGRAM  *#* 


TEST  RUN  3  282  TH  DAY  OF  1976  AT  1000  HOURS 

MODEL  PLAN  VER  GAGE  ARR.  1 

HOR.  SCALE  1-300  VER.  SCALE  1-  60  WATER  TEMP.  67.00 

MODEL  PROTOTYPE 

PERIOD  1154.5  SECONDS  12.4  HOURS 

15  CHANNELS  OF  60  SAMP /PER.  FOR  1  PERIODS 
•REC./PER.  2 


UPDATE 

CALIBRATIONS 

INTERCEPT 

SLOPE 

SCAN  * 

1 

COEF 

M 

.21 18377E-01 

. 5700070E-03 

SCAN  * 

2 

COEF 

U 

.2210692E-01 

•5698883E-03 

SCAN  * 

3 

COEF 

9 

.21 18377E-01 

.5700070E-03 

SCAN  * 

4 

COEF 

• 

.2210692E-01 

.5698883E-03 

SCAN  * 

5 

COEF 

* 

.21 18377E-01 

.5700070E-03 

SCAN  ♦ 

6 

COEF 

« 

.2210692E-01 

•5698883E-03 

SCAN  # 

7 

COEF 

m 

.21 18377E-01 

.5700070E-03 

SCAN  * 

8 

COEF 

m 

.2210692E-01 

.5698883E-03 

SCAN  * 

9 

COEF 

M 

.21 18377E-01 

.5700070E-03 

SCAN  ♦ 

10 

COEF 

m 

.2210692E-01 

.5698883E-03 

SCAN  * 

11 

COEF 

m 

.2210692E-01 

.5698883E-03 

SCAN  * 

12 

COEF 

m 

.21 19377E-01 

.5700070E-03 

SCAN  * 

13 

COEF 

M 

.2210692E-01 

.5698883E-03 

SCAN  * 

14 

COEF 

M 

.2210692E-01 

.569B8B3E-03 

SCAN  * 

15 

COEF 

m 

.2210692E-01 

•5698883E-03 

SCAN  * 

16 

COEF 

m 

.2210692E-01 

.5698883E-03 

SCAN  * 

17 

COEF 

M 

.2210692E-01 

.56988B3E-03 

SCAN  * 

18 

COEF 

m 

. 239670 IE-01 

.5700070E-03 

SCAN  * 

19 

COEF 

m 

. 2297306E-0 1 

. 568369  IE-03 

SCAN  * 

20 

COEF 

m 

. 239670 IE-01 

.5700070E-03 

SCAN  * 

21 

COEF 

at 

. 239670  IE-01 

. 570087BE-83 

SCAN  ♦ 

22 

COEF 

* 

.2210692E-01 

.5698883E-03 

SCAN  * 

23 

COEF 

m 

.23101 15E-01 

.5715380E-03 

SCAN  * 

24 

COEF 

m 

.23101 15E-01 

.5715380E-03 

SCAN  ♦ 

25 

COEF 

m 

. 239670  IE-01 

. 5700070E-03 

SCAN  ♦ 

26 

COEF 

M 

. 239670  IE-01 

.5700070E-03 

SCAN  # 

27 

COEF 

m 

. 239670  IE-01 

.5780070E-03 

SCAN  * 

28 

COEF 

■ 

.2488957E-01 

.5698883E-03 

SCAN  * 

29 

COEF 

m 

•2488957E-01 

.5698883E-03 

SCAN  ♦ 

30 

COEF 

n 

. 239670  IE-01 

.5700070E-03 

SCAN  • 

31 

COEF 

m 

. 2488957E-0 1 

.56S8883E-03 

SCAN  * 

32 

COEF 

■ 

.2488957E-01 

.5E98883E-03 

SCAN  ♦ 

33 

COEF 

* 

.2396701E-01 

.5700070E-03 

SCAN  * 

34 

COEF 

* 

.2574830E-01 

. 566369  IE-03 

SCAN  * 

35 

COEF 

m 

.2574830E-01 

. 568369  IE-03 

SCAN  * 

36 

COEF 

• 

.2488957E-01 

.5698883E-03 

SCAN  * 

37 

COEF 

■ 

.2589186E-01 

.5715380E-03 

SCAN  ♦ 

38 

COEF 

at 

.2574830E-01 

. 568359  IE-03 

SCAN  • 

39 

COEF 

■ 

.2574830E-01 

. 568369  IE-03 

PLATE  12 


***  TIDE  ANALYSIS  PROGRAM  w 
MODEL  TIDE  DATA 

CHANNEL  1  GAGE  ARRANGEMENT  1  RUN  NUMBER  3 

XY  COORDINATES  .0000  .0000 

(JATER  DEPTH  IS  -.1200 

GAGE  *  IS  1 


STATISTICS  OF  THE  ORIGINAL  DATA 

MEAN  -  1679494E+00  VARIANCE  -  .  .843I425E+00 


****  HARMONIC  ANALYSIS  **** 

(HRS)  (FT)  (DEG) 

PERIOD  AMPLITUDE  PHASE 

1.242E+01  1.287E+00  9.165E+81 

6.210E+00  1 .485E-02  -7.290E+01 

4. 140E+00  1.558E-02  -1.158E+02 


JTATISTICS  OF  THE  RESIDUAL  DATA 

MEAN  -  1285430E-04  VARIANCE  •  .  .2595286E-03 
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*#*  TIDE  ANALYSIS  PROGRAM  *** 

MODEL  TIDE  DATA 

CHANNEL  1  GAGE  ARRANGEMENT  1  RUN  NUMBER  3 

XY  COORDINATES  .0000  .0000 

WATER  DEPTH  IS  -.1200 

GAGE  •  IS  1 


STATISTICS  OF  THE  ORIGINAL  DATA 

MEAN  -  -. 1592577E+00  VARIANCE  -  ,  .8496584E+00 


****  HARMONIC  ANALYSIS  **** 

(HRS)  (FT)  (DEG) 

PERIOD  AMPLITUDE  PHASE 

1.242E+01  1.298E+00  9.131E+01 

6.21 0E+00  1 . 587E-02  -8 . 70 1 E+0 1 

4. 140E+00  1.476E-02  -1.067E+02 

3. 105E+00  1.451E-02  4.761E+01 

COMPUTED  SHIFT  IN  MODEL  HOURS  -  .4510917E-01 

FORCE  SHIFT  IN  DEGREES  -  .9000000E+02 


STATISTICS  OF  THE  RESIDUAL  DATA 

MEAN  -  . 6143107E-03  VARIANCE  -  .  .24646B5E+01 

STATISTICS  OF  THE  TAPERED  RESIDUAL  DATA 

MEAN  -  - . 5559496E+00  VARIANCE  -  .  .3516450E+01 


PLATE  15 


PROTOTYPE 


.00 

(HRS 


12.00  16.00  20.00  24.00 

TIME  PROTOTYPE  (HRS) 


16.00 


RUN  NUMBER  3  GRGE  NUMBER 


000 '03 


000'9l 

S  ** 


000 '51 

(01  ) 


000'Q 


A933N3 


000 'h 


000' 0 


GfiGE 


a  . 

1 


Constant  value  of  i  component 

a  Mean  of  water-surface  elevation  data  series 

o 

[Aj  Vector  of  coefficients 

AR  Actual  reading 

2  2  th 

Amp^  Amplitude  (a^  +  b^)  for  i  component 
Constant  value  for  i^  component 
( C ]  Real,  symmetric  matrix  of  sine  or  cosine  terms 
E  Error  that  least  squares  procedure  minimizes 

f  Water-surface  elevation  data  series 

f'  Water-surface  elevation  data  series  less  noise  associated  with  trans¬ 
fer  function 

fftj  Real  part  of  FFT 

fft^  Imaginary  part  of  FFT 

[F]  Vector  of  products  of  observed  signal  and  sine  or  cosine  terms 
FFT  Results  of  Fast  Fourier  Transform 

H  Average  wave  height  for  record 

H  Average  wave  height  of  highest  x  percent  of  wave  height,  e.g., 

is  significant  wave  height 

H  Root-mean-square  wave  height  for  record 

i  Integer  variable 

j  Integer  variable 

J  Total  number  of  components 

KE  Prototype  water  elevation  at  time  zero  reference  readings  made 

Semidiurnal  principal  lunar  tidal  harmonic  component 
MSF  Model  scaling  factor 

MWL  Reference  mean  water  level 

n  Integer  variable 

N  Total  number  of  data  samples 

PHSi  Phase  [ARCTAN  (b^a^)]  for  i*"*1  component 
PSD  Result  of  power  spectral  density 

SE  Scaled  water  elevation 

SFC  Scale  factor  constant 

SFS  Scale  factor  slope 

SF1  Scale  factor  constant 

SF2  Scale  factor  linear  term 

SF3  Scale  factor  quadratic  term 


A2 


Wave  period 

Wave  period  for  component 

Average  wave  period  for  record 
Unsealed  water  elevation 
Virtual  reading 

Specified  percent  of  highest  wave  heights,  normally  x 

Cartesian  coordinate 

Cartesian  coordinate 

Water  depth,  ft 

Zero  reference  reading 

Time  interval  for  consecutive  observation 
Noise  associated  with  transfer  function,  etc. 

Standard  deviation  of  average  wave  height 
Standard  deviation  of  average  wave  period 
2n/T.  ,  radian  frequency 


nnnnnnnnnnnonnnnnnnnnnnnnonnn 


PACE  1  C  UAVE  DATA  ANALYSIS  PROGRAM 

FORTRAN  COMPILER  REV.  LEV.  J00 


SENSE  SWITCHES  ARE  SWITCHES  LOCATED  ON  THE  CONTROL  PANEL 
OF  THE  PACER  100.  THEY  CAN  BE  EITHER  ON  OR  OFF  (.TRUE.  OR 
.FALSE.).  FOR  OTHER  MACHINES  SENSW  SHOULD  BE  DEFINED  AS 
BEING  AN  EIGHT  ELEMENT  LOGICAL  ARRAY#  AND  SHOULD  BE 
INITIALIZED  BY  THE  USER. 

SENSE  SWITCH  A  FOR  LIMITED  OUTPUT 

SENSE  SWITCH  B  TO  CHANGE  •  PERIODS  TO  ANALYZE 

SENSE  SWITCH  C  TO  PLOT  DATA 

SENSE  SWITCH  D  TO  USE  DISK  SCRATCH  FILE 

SENSE  SWITCH  E  ENDS  AUTOMATIC  MODE 

SENSE  SWITCH  F  TO  STORE  RESULTS  ON  MAG  TAPE 

SENSE  SWITCH  G  &  A  FOR  SUrtlARY  OUTPUT 

SENSE  SWITCH  H  TO  SKIP  GAGES 

MAXIMUM  LIMITS  FOR  THIS  PROGRAM 
100  WAVE  PERIODS 
6000  DATA  POINTS  PER  CHANNEL 
1920  DATA  POINTS  PER  RECORD  ON  MAG  TAPE 
64  GAGES  OF  DATA 

THE  FOLLOWING  DEFINE  STATEMENTS  ARE  PECULIAR  TO  THE  EA I  PACER 
100  COMPILER.  THE  DEFINE  STATEMENT  USED  IN  THE  FOLLOWING 
MANNER  FORCES  THE  VARIABLE  DEFINED  TO  BE  ASSIGNED  TO  THE 
OCTAL  MEMORY  CELL  CONTAINED  IN  PARENTHESIS.  MEMORY  AND 
RUN  TIME  ARE  SA\«D  BY  USING  THIS  STATEMENT  ON  THE  PACER  100. 
THE  STATEMENTS  MAY  BE  OMITTED  FOR  USE  ON  ANY  OTHER  MACHINE. 

DEFINE  BUF (*200) 

DEFINE  Il(* 254) 

DEFINE  I0DEV<*255) 

DEFINE  NPCT<*256) 

DEFINE  IU1000<*25?) 

DEFINE  1ERR<*260) 

DEFINE  IRL('26l) 

DEFINE  NRSKIP<*262) 

DEFINE  IH0RZ<*263) 

DEFINE  IVERT<*264) 

DEFINE  1RUNN0(*265) 

DEFINE  NCH<*266) 

DEFINE  NRECPD<*267) 

DEFINE  IGAGAR(*270) 

DEFINE  ISPUPD<*271) 

DEFINE  N0UPD<*2?2) 

DEFINE  NSPTST<*273) 

DEFINE  MM<*274) 

DEFINE  NR0DS<*275) 

DEFINE  NSAMP<*276) 

DEFINE  LCK*277) 

DEFINE  LC2<*300) 

DEFINE  NCHAN<*301) 

DEFINE  NNC302) 

DEFINE  ICC (*303) 


PAGE  2 


C 


UAVE  DATA  ANALYSIS  PROGRAM 


'1 

!*> 

i 


K 

'j 


DEFINE  JJJJ('304) 
DEFINE  1RLD<’305) 
DEFINE  II <'306) 
DEFINE  IERCD<*307) 
DEFINE  I2<'310> 
DEFINE  I3<'311) 
DEFINE  K< '312) 
DEFINE  N<'313) 
DEFINE  IR<'314) 
DEFINE  1P<'315) 
DEFINE  IIII<*316> 
DEFINE  IV<'317) 
DEFINE  19 <'320) 
DEFINE  NC0EF<'321) 
DEFINE  NCALV<'322) 
DEFINE  NC<*323) 
DEFINE  NPER < ' 324) 
DEFINE  IPOTPT < ' 325) 
DEFINE  I <'326) 
DEFINE  1ER<'32?) 
DEFINE  NPTPLT<'330) 
DEFINE  I8<'331) 
DEFINE  ICHSKP<*332) 
DEFINE  IST<'333) 
DEFINE  14<'334) 
DEFINE  J<'335) 
DEFINE  IF1<'336) 
DEFINE  IF <'337) 
DEFINE  IF3<*340) 
DEFINE  KK3<'341) 
DEFINE  KKl<’342) 
DEFINE  KK2<'343) 
DEFINE  1 1 1 < ' 344) 
DEFINE  J1<'34S) 
DEFINE  111 <'346) 
DEFINE  IBX<'347) 
DEFINE  NS <'350) 
DEFINE  NPS<'35l) 
DEFINE  KGAGE < ' 352) 
DEFINE  I GAGE <’353) 
DEFINE  IU1020<'354) 
DEFINE  IU1030<'355) 
DEFINE  IU1040<’356) 
DEFINE  1 1 1 <'357) 
DEFINE  L2<*360) 
DEFINE  L2P<'361) 
DEFINE  I2D<*362) 
DEFINE  L<’363) 
DEFINE  NPEAK<*364) 
DEFINE  NTROF<*365) 
DEFINE  IF2<'366) 
DEFINE  NP<'367) 
DEFINE  ICV(’370) 
DEFINE  IXX<’371) 
DEFINE  Jll<'372) 
DEFINE  K1 1<'373) 
DEFINE  IJK< '374) 


B4 


.'.V-j 


. „■ 


*1 


y:.i 

-*  i 
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C 


DAVE  DATA  ANALYSIS  PROGRAM 


DEFINE  NPK(’375) 

DEFINE  K9(*376) 

DEFINE  I FLAG ('377) 

DIMENSION  RESULT (64) 

DIMENSION  IPOITl(ll).  IPOIT2(10) 

DIMENSION  NAME (3) ,  ISAVEO) ,  I DATA (6000), IBUF <  1920). 2P<2) . 

&IX>«( 100) .TROFPS( 100).PEAKPS( 100) ,KT1TLE(20)» 1YYY( 100) , 
&C00R(3.64) .P0T(64) ,X( 1 1) ,2Y( 1 1) .S2( 1 1) . 

&>«X( 100) . YYY( 100) ,DATA(3000) ,2(3000) . IX( 100) , IY( 100) , IUC100.2) 
& IGAGEN (2, 64) , BUF (22) , I FLAG (2,64). COEF (11.64), ICALV( 42. 50) . 

&UH( 100) 

DIMENSION  DIST(50) ,DEPTH(50) . IKFTI 1 (4.2) . IKFTI2(6.2) 
EQUIVALENCE  (X( 1) ,BUF( 1) ) . (BUF( 12) ,2Y( 1 ) ) . ( IDATA( 1 ) , 

&DATA( 1) .2(1) ) , ( IX( 1) . 

21BUF ( 1 ) ) , ( IY( 1 ) , I BUF ( 101) ) . ( IU( 1 . 1) . IBUF (201 ) ) . 

3(UH(1). IBUF(401)), 

&OW(  1 ) ,  IBUF (60 1 ) ) ,  (YYY(  1 )  .  IBUF (80 1 ) ) , 

5(IYYY(1),IBUF( 1001) ) , ( 1XXX( 1) . IBUF( 1101)). 

6 ( TROFPS ( 1 ) . IBUF ( 1 20 1 ) ) . ( PEAKPS ( 1 ) . I BUF ( 1 40 1 ) ) 


EQUIVALENCE  (UVHGT, IDATA(66) ) . (tJMLL. 1DATA( 130)) 

LOGICAL  I OUT, ICAL.SENSU. ICHSKP 

LOGICAL  ISFRCT 

SCALED  FRACTION  IFRACT 

DATA  IUNIT2/' 14/. PCT/33. 33333333333/ 

DATA  NAME(1)/2HA1/,NAME(2)/2H11/,NAME(3)/2H11/ 

DATA  SAME/9999./ 

DATA  IRESP/1HY/ 

DATA  IP0IT1/15. 14. 13. 12. 11.0.1.2.3.4,5/ 

DATA  IPOIT2/0. 16. 17. 18. 19. 10.9.8.7.6/ 

DATA  PIE/3. 141592654/, HBP IE/6 . 283 1 85307/, FORP IE/12.56637061/ 
DATA  IKFTI 1/2H  ,2H  ,2H  ,2H  .2HAD.2HJU.2HST.2HED/ 

DATA  IKFTI2/2H  ,2H  ,2H  ,2H  ,2H  ,2H  .2HUA.2HVE.2H  H.2HEI. 

&2HGH.2HT  / 

DATA  KTITLE/2H  ,2H  ,2H  ,2H  ,2H  ,2H  ,2H  ,2H  ,2H  , 

&2H  ,2H  .2H  ,2H  ,2H  ,2H  .2H  ,2H  ,2H  ,2H  ,2H  / 

GO  TO  8 

1  FORMAT ( 15HAUT0MATIC  MODE?,/) 

FORMAT ( /, 7HDEPTH -? , /) 

FORMAT(F10.5) 

2  FORMAT(Al) 

1  FORMAT ( 20 A2) 

9  FORMATOIHANY  OUTPUT?./) 

3  FORMAT(10HFILE  NAME?./) 

6  FORMAT (1H?) 

FORMAT( 1H1. 17X.28H#**  UAVE  RECORD  ANALYSIS  ***./, 85X. 

&5HM0DEL, 1 IX. 9HPR0T0TYPE , /, 1X.20A2.29X.7HPERIOD  , 
ft2(4X.F7.3,0H  SECONDS) ./.9H  TEST  RUN.2I6.14H  TH  DAY  OF  197. 

& 1 1 . 3H  AT.I5.6H  HOURS, 20X, 10HUAVEHEIGHT.F8.3.7H  INCHES. 6X, 
&F6.3.5H  FEET./.12H  MODEL  PLAN  .4A2.3X.9HGAGE  ARR..I5.3X. 
&9HGEN.  POS.,16. 15X.12.12H  CHANNELS  0F.13.14H  SAMP/PER.  FOR. 
fcl4..8H  PERIODS./. 14H  HOR.  SCALE  1-.I3.15H  VER.  SCALE  1«. 

&I3. 13H  UATER  TEMP. .F6.2. 16X, 10H*REC./PER.. I5.3X. 

&7HSUL  ADJ.F6.3) 

F0RMAT(1H1.1X. 

X52H###  UAVE  ANALYSIS  BASED  ON  SIMPLE  SEARCH  TECHNIQUES  . 
«,3H)W*./,  1X.20A2./.8H  CHANNEL. 2X.  13, 5X.9HGAGE  ARR.  ,2X,  14./) 
FORMAT ( 15H  XY  COORDINATES. 2X,2(F9. 4.2X) ,/. 15H  UATER  DEPTH  IS. 
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8.2X.F9.4./.12H  GAGE  #  IS  .2A2) 

68  FORMAT*/. 10X.7H  CRESTS. 14X.SH  TROUGHS. 13X. 1 1H  UAVEHEIGHT./.3X 
13<3X.6H  WlLUE.2X.9H  POSITION, 2X) ) 

69  FORMAT* 1X.5H*****,  15.36H  TH  UAVEHEIGHT  >  BREAKING  UAVEHE1GHT) 
7000  FORMAT*/, 14HDISK  I/O  ERROR) 

5012  FORMAT* 1H  ,6<F10.3. IX) . IX. 14) 

5011  FORMAT* 1H  ,//,23H  SEARCH  ERROR— PERIOD  ,15) 

5013  FORMAT* 1H  /.21X. 14HM0DEL  -  INCHES. 15X, 16HPR0T0TYPE  -  FEET. 
8./.12H  UAVEHEIGHTS./.5X, 10HAVERAGE  ,2*F10.3,5H  +/-  ,F10.3, 
&5X) ./.5X.5H  RMS.13X,F10.3,20X.F10.3.//.5X,F5.2.5H  PCT  . 
&2(F10.3,5H  +/-  .F10.3.5X) , 

8./.21X. 15HM0DEL  -  SECONDS, 1 IX, 19HPR0T0TYPE  -SECONDS. 

&//, 11H  UAVEPERIOD,/. 10X,2*5X.F10.3.5H  +/-  .F10.3)) 

C  INPUT  TAPE  UNIT 
8  TYPE  6201 

6201  FORMAT*/, 11HTAPE  UNIT-?./) 

ACCEPT  6202,1 UN  IT 

6202  F0RMAT<08) 

IF* IUNIT.EQ. *  14) IUNIT2-' 15 
ICHSKP-.FALSE. 

C  LIST  CHANNELS? 

TYPE  8209 

6209  FORMAT*/, 14HLIST  CHANNELS?./) 

ACCEPT  2.1 

IF* I .EQ. IRESP) ICHSKP-.TRUE. 

I OUT-. FALSE. 

IODEV-16 
C  ANY  OUTPUT? 

7060  TYPE  9 

ACCEPT  2.1 

IF* I .EQ. IRESP) I OUT-. TRUE. 

C  INPUT  TITLE 
TYPE  8031 

8031  FORMAT*/, 7HTITLE-?./) 

ACCEPT  3311, <KT!TLEU). 1-1.20) 

T  CPPPT*  P^l  CP 

C  DATA  SCALED  FRACTIONS? 

TYPE  7381 

7381  FORMAT*/. 22HDATA  SCALED  FRACTIONS?./) 

ACCEPT  2. 1 

IF* I. EQ. IRESP)  1SFRCT-.TRUE. 

C  INPUT  MODEL  DEPTH 
TYPE  900 

ACCEPT  901, DEPTH 1 
CALL  SELMT* IUNIT) 

IRL-20 

C  INPUT  DATA  FILE  NAME 
TYPE  3 

ACCEPT  3311, *IDATA*I). 1-1.3) 

C  POSITION  TO  DATA  FILE 

CALL  POSMT* 1DATA<1).3. IRL. .FALSE. • IERR) 

IF* IERR.EQ.0)GO  TO  5 
7  TYPE  6 
J  .8 
5  IRL-20 
C  READ  HEADER 
C  READ  FIRST  HEADER 
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78  CALL  RDMT( IDATA<51), IRL. .TRUE. . IERR) 

L  .99 
ICAL-.TRUE. 

IFi IDATA<65) .GE.0) ICAL-. FALSE. 

C  READ  SECOND  HEADER 

CALL  RDMTC I DATA ( 101). IDATA<70) . .TRUE..IERR) 

IF< IERR.NE.0)GO  TO  7 
IF(.NOT.SENSU<6))GO  TO  7384 
C  IF  OUTPUTTING  TAPE  URITE  FILE  NAME 
CALL  SELMT( IUNIT2) 

CALL  Nl-f  MT< 1DATA< 1) .3. IRL. .FALSE. . IERR) 

L  .99 

C  OUTPUT  FIRST  HEADER 

CALL  URTMTUDATA<51).  IRL.. TRUE.,  IERR) 

L  .99 

C  OUTPUT  SECOND  HEADER 

CALL  URTMT< IDATA( 101) , IDATA(?0) . .TRUE., IERR) 

L  99 

CALL'SELMT(IUNIT) 

7384  NRSKIP-0 

IF ( .N0T.SENSLK2) )GO  TD  111 

C  INPUT  NUMBER  OF  PERIODS  TO  ANALYZE  AND  NUMBER  OF  RECORDS  TO 
C  SKIP 

TYPE  6299 

6299  FORMATS,  2 1H#R£C  &  #PER  TO  SKIP-?./) 

ACCEPT  6252. IDATA<105).NRSKIP 
6252  F0RMAT(2I7) 

C  RETRIEVE  USABLE  INFORMATION  FROM  FIRST  TTJO  HEADERS 
111  LC1-IDATA<54)/1000 

LC2-1DATA<54)-LC 1*1000 
CUL“FLOAT( IDATA<62) )/100. 

TEN>-FL0AT(IDATA(63) )/10. 

UML-UMLL 

KFACT-2 

IF< I DATA  ( 129) .GE.0)KFACT-1 
UAVPD“DATA<51) 

IHORZ-IDATA(60) 

IVERT-IDATA<61) 

C  COMPUTE  SCALING  FACTORS 
H0R2 -FLOAT ( IHORZ) 

VERTFT “FLOAT  < I VERT) 

OMEGA -VERTFT/HORZ 
SUL -5. 1 2*UAVPD*UAVPD 
C-SUL 

902  ARG-TUOP IE*DEPTH 1/SUL 
SUA “C*TANH ( ARG ) 

I F  <  ABS  <  SUA-SUL ) -0 . 00 1  *OJL )  903 . 903 , 904 
904  SUL-(SUA+SUL)/2. 

J  .902 

903  THARG -TANH ( ARG ) 

1RUNNO- I DATA (53) 

ARGU-ARG*OMEGA 
THARGU-TANH ( ARGU) 

TIME “SORT (HORZ*THARG/THARGU) 

VERT-VERTFT/12. 

SA-UAVPD*T I ME 
RR“UVHGT*VERT 
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NCH-IDATA(64) 

NRECPD - I DATA ( 1 04) 

1GAGAR«IDATA(116) 

ISPUPD»IDATA< 103) 

NOUPD-IDATAO05) 

NSPTST  ■ ISPUPD*NOUPD 
C  COMPUTE  RECORD  LENGTHS 
NRODS-IDATA( 110) 

NSAMP ■ I SPUPD/NRECPD 
1RLD- ISPWPD*NCH/NRrCPD 

C  DETERMINE  KENETIC  VISCOSITY  OF  FRESH  WATER  FOR  COfPUTING  KEULEGAN 
C  FACTORS 

AA -TEMP-60. 

BB-ABS(AA) 

IF (BB- 10 .) 5050,5050. 5005 
5005  DD-BB-10. 

IF (AA) 5014. 5050. 5021 
5014  VISC-(1.41+(.0254*DD>)*. 00001 
J  .5060 

5021  VISC- < 1 . 059- ( . 0 129*DD) )*. 0000 1 
J  .5060 

5050  VISC-( 1. 217- (.01?6*AA)+( .00 1?*8B))». 00001 

5060  IF( .NOT. IOUT)GO  TO  13 

C  LIST  INFORMATION  FROM  FIRST  TWO  HEADERS 

12  WRITE  (I0DEV.14)  (KTITLE(KKI) .KKl -1 ,20) .UAVPD.SA. IRUNNO. 

&LC2. LC 1 . I DATA (55) . UVHGT. RR . ( I DATA (KK 1 ) . KK 1 -56. 59) . IGAGAR, 

& I DATA (115). NRODS . ISPUPD.NOUPD. IHORZ, I VERT. TEMP. NRECPD. CUL 
IF(SENSW(7))UR1TE  (I0DEV.6992)  ( IKFTI 1 ( I.KFACT). 1-1,4) , 
&(IKFTI2(I.KFACT). 1-1.6) 

6992  FORMAT (//.  16X.27HUAVEHE1GHT  PROTOTYPE  -  FEET.20X. 

&20HUAVEPER I OD  -  SECONDS . 1 4X. 4A2 , /. 1 X. 7HCHANNEL , 4X, 

&14H AVER AGE  HEIGHT. 

&9X. 1 1HHIGHEST  1/3. 10X, 12HMODEL  PERI0D.7X. 16HPR0T0TVPE  PERIOD, 
&6X.6A2. 1 IX, 12HCH  CODE  PER) 

C  IF  NOT  CALIBRATION  FILE  SKIP  READING  UNNECESSARY  RECORDS 
13  IF(ICAL)GO  TO  1? 

J  .16 

C  READ  TRANSDUCER  FLAGS 
1?  IRL-NCH*2 

CALL  RDMT( I FLAG (1, 1), IRL. . TRUE . . IERR) 

L  .99 

C  COMPUTE  NUMBER  OF  GAGES  HAVING  CALIBRATION  COEFFICIENTS 
NCOEF-0 
NCALV-0 

DO  6291  l-l.NCH 
L  ,.26 

6291  CONTINUE 

C  READ  CALIBRATION  COEFFICIENTS 
IF (NCOEF .EQ . 0) GO  TO  6292 
K-NCH*22 

CALL  RDMT(COEF(l.l).K.. TRUE.. IERR) 

L  99 

C  READ  GAGE  NAMES 

6292  CALL  RDMT(IGAGEN( 1.1). IRL.. TRUE.. IERR) 

L  .99 

C  READ  WEIGHTED  MEAN  DEPTHS 
IRL-NRGDS*2 


PAGE  7  C  UAVE  DATA  ANALYSIS  PROGRAM 


IF<KFACT.NE.2)G0  TO  8987 
CALL  RDMT<D£PTH. IRL. .TRUE. . IERR) 

L  .99 

C  READ  DISTANCES  FROM  UAVE  GENERATOR 
CALL  RDMT(D 1ST. IRL.. TRUE.. IERR) 

L  .99 

C  READ  GAGE  COORDINATES 
8987  K«NCH*6 

CALL  RDMT<COOR(l.l),K.. TRUE. .IERR) 

L  .99 

C  READ  POTENTIOMETER  COEFFICIENTS 
IF (NRODS.EQ.B)GO  TO  20 
CALL  RDMT<POT(l). IRL.. TRUE.. IERR) 

L  .99 

C  READ  CALIBRATION  VOLTAGE  READINGS 
28  I F (NCALV . EO . 0) GO  TO  16 

IRL-NCALV*42 

CALL  PDMT< ICALV< 1 . 1 ) . IRL . . TRUE . . IERR) 

L  .99 

16  IF ( .NOT. IOUT)GO  TO  76 

IF  < NRSK IP . LE . 0) GO  TO  8391 
C  SKIP  RECORDS  IF  ANY  NOT  DESIRED  IN  ANALYSIS 
DO  8392  K-l.NRSKIP 
CALL  SRFMT(. TRUE.. IERR) 

L  99 
8392  CONTINUE 

C  INITIALIZE  COUNTERS  FOR  GAGE  PROCCESSING 
8391  NCOEF-0 
NCALV-0 
IPOTPT-0 

NPER-NOUPD*NRECPD 
IF ( . NOT. SENSU<4) ) GO  TO  7001 
C  COPY  DATA  TO  TEMPORARY  DISK  FILE 
CALL  SELMD<'21.0) 

CALL  NMrMD(NAME< 1) .3.88. IER) 

IF< < IER.EQ.0) .OR. ( IER.EQ.6) )GO  TO  7002 
7004  TYPE  7000 
J  .7 

7002  CALL  SAVEMD  < I SAVE < 1 ) ) 

DO  7003  Il-l.NPER 

CALL  RDMT<IBUF(1).IRLD,. TRUE. .IERR) 

L  SO 

CALL ’ URTMD ( IBUF ( 1 ) . IRLD. . TRUE . . IERR) 

L  .8999 

7003  CONTINUE 
IF(IER.EQ.6)G0  TO  7001 
CALL  UEFMD(. TRUE., IERR) 

L  .8999 

C  COMPUTE  NUMBER  OF  POINTS  TO  PLOT.  CAN  PLOT  NO  MORE  THAN  3000  PTS 
7001  NPTPLT -NSPTST 

IF (NSPTST. GT. 3000) NPTPLT-3000 
DELT-UAVPD/FLOAT ( ISPUPD)*TIME 
ZP<1)«0. 

ZP(2) -DELT*FLOAT (NPTPLT) 

DO  40  I-l.NCH 

C  BUMP  POINTERS  FOR  CALIBRATION  COEFICIENTS  fc  CALIBRATION  VOLTAGES 
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NTROF-0 

NPEAK-0 

IF (.NOT. ICHSKP)GO  TO  8110 
TYPE  8010.1 

8010  FORMAT (/'.3HCH  .13,/') 

C  INTERPRET  FLAGS  FOR  THIS  GAGE 
8110  IF1-1FLAG(1.1)/100 

IF2-1FLAG( 1 , I)-( IF1#100) 

IF(IF2.EQ.1)G0  TO  42 
GO  TO  40 

42  IPOTPT-IPOTPT+1 

IF2-(IFLAG(2. 1)/100)*100 
IF3-((1FLAG(2,1)/'10)*10)-IF2 
IF2-IFLAG(2.  D-IF2-IF3 
IF(IF2.EQ.2)G0  TO  40 
IF3-IF3/'10+1 
IF( 1F3.EQ. 10) IF3-4 
C 
C 

C  MOVE  CALIBRATION  COEFFICIENTS  TO  APPLICABLE  BUFFERS 
DO  41  1 1-1- 1 1 

41  S2  < 1 1 ) -COEF (II, NCOEF ) 

IF(IF3.LE.3)G0  TO  3303 

K-l 

L-l 

L  .2100 

K-12 

L-22 

L  .2100 
DO  48  11-1,11 

48  ZY( 1 1) -ZY( I l)#POT( IPOTPT) 

C  SKIP  THIS  GAGE 

3303  IF(SENSU(8) )GO  TO  40 
13-1 

1F(.N0T.SENSU(4))G0  TO  49 
C  DEMULTIPLEX  GAGE  DATA  FROM  DATA  RECORDS 
CALL  RESETD ( 1SAVE ( 1 ) . . TRUE . ) 

49  DO  43  Il-l.NPER 
IRL-IRLD 

IF(.N0T.SENSU(4))G0  TO  7007 
CALL  RDMD( IBUF( 1) , IRL. .TRUE. , IERR) 

L  .8999 
J  .7008 

7007  CALL  RDMT(IBUFd),  IRL,.  TRUE.,  IERR) 

L  .99 

7008  DO  44  14-1. IRLD.NCH 

IF ( .NOT. ISFRCT)GO  TO  8001 
J-IBUF( 14) 

LA  J 

STA  IFRACT 
SOUT-IFRACT 
SOUT-SOUT* 10000. 

IBUF( 14) -SOUT 
6001  IDATA( I3)-IBUF(I4) 

44  13-13+1 
43  CONTINUE 

IF(SENSU(4))GO  TO  957 


BIO 
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C  REWIND  TAPE  TO  BEGINNING  OF  DATA  RECORDS 
7086  DO  45  Il-l.NPER 

CALL  SRRMT < .TRUE . » I ERR) 

L  99 
45  CONTINUE 
957  III-I 

C  INITIALIZE  COUNTERS  FOR  WAVE  ANALYSIS  ROUTINE 
IERCD-8 
IBX-1 
NS -NO UP D 
NPS-1SPUPD 

6110  BREAK “COOR (3. 1 1 1 ) /l . 28 
IF(SENSU<?) )GO  TO  943 
S“C00R<3.  D+CUL 
C  LIST  GAGE  INFORMATION 

URITE  < I0DEV.66)  <KTITLE(KK1>.KK1-1.20). III.IGAGAR 
URITE  < IQDEV.67)  COORC 1 , 1) ,COOR<2. 1) ,S. IGAGEN< 1, I) , 
ftIGAGEN(2.I> 

943  IF<.NOT.SENSU(l))URI7E  (I0DEV.68) 

IXXX<1) “-32766 

C  WAVE  HEIGHT  ANALYSIS  ROUTINE  BEGINS  HERE 
C  FIND  FIRST  PEAK  UI7HIN  PERIOD  UIDTH 
DO  6100  Ill-l.NPS 

IF(IXXX(1)-IDATA< 1 11) >921.6100,6100 

921  IXXX(l)-IDATAdll) 

IX(1)-I11 

6100  CONTINUE 

C  COMPUTE  SEARCH  UINDOU  UIDTH,  30X  MORE  THAN  PERIOD 
LZ«0.3*FLOAT<NPS)40.5 
LZP-NPS+L2 

C  LOOP  TO  FIND  MAXIMUM  AND  MINIMUM  PAIRS  (LOOKING  FOR  NS  PAIRS) 
DO  6200  Ill-l.NS 
IX<Ill+l)-0 
IY< 1 1 1+1) -0 
IZD-0 

C  FIRST  PAIR? 

IFdl 1-1)922,922,6210 

C  BRACKET  SEARCH  UINDOU  6  ELEMENTS  BEFORE  FIRST  MAXIMUM  FOUND 

922  J-IX(Ill)-6 
IF(J.LT.1)J-1 

j  6220 

C  BRACKET'SEARCH  UINDOU  2  ELEMENTS  PAST  LAST  MINIMUM  FOUND 
6210  J-IY(Ill-l)+2 
6220  KK1-J+LZP 

K-MIN0(KK1,NSPTST) 

I F  <  J +NPS-NSPTST) 6230,6230, 6320 
6238  IXXX< Ill) “-32766 
IYYY< 1 1 1) -32766 

C  FIND  MAXIMUM  &  MINIMUM  UITHIN  UINDOU 
DO  6240  L-J.K 

IF ( IXXX( 1 1 l)-lDATA(L) ) 923,6250.6250 

923  IXXX( Ill) “IDATA(L) 

IX< 1 1 1 >  — L 

J  .6240 

6250  IFdYYYd  1  l)-IDATA(L))  6240,6240,924 

924  IYYYdll)“IDATA<L) 

IY(I11)-L 
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6240  CONTINUE 

C  CHECK  CREST  -  TROUGH  PAIR  FOR  ACCEPTABILITY 
C  DOES  MINIMUM  FOLLOU  MAXIMUM? 

IF( IYYYC 1 1 1 )-IXXX< Ill) >925,925,6500 
C  FIRST  PAIR? 

925  IF<I 11-1)926,926,6260 

C  REDUCE  SEARCH  BY  1  ELEMENT 

926  K-IX<I11>-1 
J  .6230 

C  DISTANCE  BETUEEN  THIS  MAXIMUM  &  LAST  MAXIMUM  ONE  INITIAL  SEARCH 
C  UINDOU  U1DTH? 

6260  IF<IXnil)-IXail-l)-L2P)92?. 927,6270 

C  REDEFINE  SEARCH  UINDOU  BEGINNING  AT  MAXIMUM  FOUND  &  GOING  1  PERIOD 
92?  lYYYd  11) -32766 
KK1-IX<111)+1 
KK2-IX<  11 D4ISPUPD 
KK2 -M I N0 ( KK2 » NSPTST) 

C  DOES  SEARCH  UINDOU  OVERLAP  END  OF  DATA  SERIES? 

IF  <KK  HNPS-NSPTST) 626 1,6261, 6320 
C  DETERMINE  MINIMUM  UITHIN  UINDOU 

6261  DO  6280  L-KK1.KK2 

IF ( I YYY( 1 1 1 ) - 1DATA <L) ) 6280 , 6280. 928 

928  IYYY<I11)-IDATA(L) 

IY< Ill) -L 

6280  CONTINUE 
J  .6500 

C  DISTANCE  BETUEEN  THIS  MINIMUM  &  LAST  MINIMUM  30X  OF  PERIOD? 

6270  IF< IY( 1 1 1)-IY< 1 1 1-1 )-L2) 6290. 6290. 929 

C  REDUCE  SEARCH  UINDOU  BY  SETTING  END  OF  SEARCH  UINDOU  TO  LOCATION 
C  THAT  MINIMUM  FOUND 

929  IXXX(Ill)— 32766 
KKl-IYdll) 

KK1-M1N0(KK1. NSPTST) 

IF  < J+NPS-NSPTST) 6262.6262. 6320 
C  FIND  MAXIMUM  UITHIN  SEARCH  UINDOU 

6262  DO  6300  L-J.KKl 

IF<IXXX< I 11)~IDATA<L))930. 6300.6300 

930  IXXX< Ill ) -IDATA<L) 

IXdll)-L 

6300  CONTINUE 
J  .6500 

C  EXPAND  SEARCH  UINDOU  BY  2  ELEMENTS  BEFORE  MINIMUM  FOUND  6  2 
C  ELEhENTS  PAST  MAXIMUM  FOUND 
6290  J-lY<Ill)+2 
K-IX(Ill)-2 

C  ADD  1  TO  ERROR  DETERMINING  COUNTER 
I2D-I2D+1 

C  IF  COUNTER  GREATER  THAN  10  OUTPUT  SEARCH  ERROR  MESSAGE 
IF ( I2D- 10)6230.6230. 6000 
6000  IERCD-2 

IF(S£NSU(?) )GO  TO  942 
URITE  (1ODEV.5011)  Ill 
942  IYYY<I11)-IDATA<L) 

C  ACCEPTABLE  PAIR  HAS  BEEN  FOUND 
6500  IU(lll.l)-(IX<lll)+lY(lll)+l)/2 
1U(I 11.2) -111 
XB-FLOAT ( I YYY( Ill)) 
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C  SCALE  NUMBERS  THAT  ARE  TO  BE  PRINTED 
GO  TO  <61. 61. 40,62), IF3 


61 

L 

.72 

L 

.3033 

L 

.72 

J 

.956 

62 

N- 

11 

L 

.74 

L 

.8033 

L 

.74 

956 

XXX < I 1 1) -FXX 

UH  <110  -XXX(1 11)-YYY(  Ill) 

C  STORE  MAXIMUM  &  MINIMUM  FOR  PLOTS 
63 10  PEAKPS (Ill)- (FLO AT ( IX< 1 1 1 )- 1 ) ) #DELT 
TROFPS (Ill) -< FLOAT < lY(Ill)-l)) *DELT 
I F  <  PEAKPS ( 1 1 1 ) . LE . ZP  <  2 ) ) NPE AK -NPEAK+ 1 
I F ( TROFPS ( 1 1 1 ) . LE . ZP ( 2) ) NTROF -NTROF+ 1 
IF(SENSU( 1) )GO  TO  6200 
TT3-<FL0AT<IU(I11. 1)-1))*DELT 
1F(UH(111).GE. BREAK)  URITE  (I0DEV.69)  Ill 
URITE  ( IODEV.5012)  XXX<  Ill)  ,PEAKPS<  Ill).YYY(Ill).TROFPSail) 
1UH<I11),TT3,IU<I11,2) 

6200  CONTINUE 
NP-NS 
J  .6330 

C  OUTPUT  ERROR  MESSAGE  THAT  FEUER  PAIRS  FOUND  THAN  SEARCHED  FOR 
6320  NP-lll-1 
IERCD-1 

IF(SENSU(?))GO  TO  6330 
URITE  (I0DEV.6263)  NP 

6263  FORMATS.  1 1H  ***  ONLY  , 14, 26H  PERIODS  OF  DATA  FOUND  ***) 

C  CALCULATE  STATISTICS  TO  BE  OUTPUT 
6330  IF < IBX) 6400, 932. 932 
932  S-0. 

SA-S 

ICV-0 

DO  1000  IXX-l.NP 
J1 1-NP-IXX 
DO  1200  Kll-l.Jll 
IJK-1XX+K1 1 

Q -FLOAT ( IX( IJK)-IX(K1 l))/FLOAT< IXX) 

RR “FLOAT < IY( IJK)-IY(K1 1) )/FLOAT< IXX) 

S-S+Q+RR 
SA -SA+Q*Q+RR#RR 
1200  ICV-ICV+2 
1000  CONTINUE 

PDEST-S/FLOAT<ICV) 

SDP-SQRT( ABS <SA-S#PDEST) /(FLOAT < ICV) - 1 . ) ) 

NPK-NP 
L  .3001 
SMEN-PMEN 
VAR-PVAR 

URMS  -SORT  <  VAR+SMEN*SftN ) 

DO  3100  IXX-l.NP 
DO  3200  K9-1.NP 

IF(UH< IXX)-UH(K9) ) 3900.3200. 3300 
3300  TT-UH(K9) 
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UH(K9)  -UH(IXX) 

UH( IXX) -TT 
KK1-IU<K9.1) 

KK2-IU(K9.2) 

1U(K9. 1 ) =]U( IXX. 1) 

11KK9.2) " 1U( IXX.2) 

HJ(IXX.1)«KK1 
ILK  IXX.2)  -KK2 
3200  CONTINUE 
3100  CONTINUE 

NPCT-IFIX<PCT*FLOAT<NPV100.0+.5> 

NPK-NPCT 
L  .3001 
SOUT -SMEN  A*ERT 
VAROU  -SORT ( VAR ) 

VARO-VAROU/VERT 
URMSO-URMS/VERT 
PMOUT-PfEN/VERT 
PVARO -SORT ( PVAR ) 

PVAO-PVARO/VERT 
PODE-PDEST*DELT/riME 
SDPO-SDP*DELT/TIME 
DELTOU»DELT*PDEST 
SDPOUT -SDP*DELT 
RESULT  <1 ) -PMEN 
IF(KFACT.NE.2)GO  TO  5061 
C  COMPUTE  KEULEGAN  COEFFICIENTS 

DEPTH  1  -DEPTH  <  IPOTPT)  +CULWERTFT 
AAL-5. 1 24789 1?*UAVPD#UAVPD 
BBB-TUOP 1E*DEPTH 1 
CCL-AAL 

5001  CCC -EXP  <  BBB/CCL ) 

SSH-l./CCC 

UL-AAL*( ( CCC-SSH )  ✓  ( CCC+SSH ) ) 

1 F ( ABS (UL-CCL) - .001) 5020. 5020*  5010 
5010  CCL-UL 

J  .5001 

5020  AAA -FORP 1 E*DEPTH 1 AL 
CCC  » ( ULAJAVPD )  *UML 
ALPHA -2 . 5*SGRT ( P I E*V I SC/UAVPD 1 
IFIAAA-10. >6020.6020.6010 
6010  ALPHA -ALPHA/CCC 
J  .6060 

6020  SINHY-EXP(AAA) 

IF(AAA-3. >6040.6040.6030 
6030  SINHY-SINHY/2. 

J  .6050 

6040  SINHY-<SINHY-I./'SINHY)/2. 

6050  ALPHA -ALPHA* (UL*SI NHY+TUOP I E*(JhL )  ✓  <  S I NHY*LL+fORP I E*DEPTH 1 > 
8./CCC 

6066  FACTOR-EXP (ALPHA*DIST< IPOTPT)) 

PhENK-PMEN*F ACTOR 
6061  IF<  .N0T.SENSIK7)  )GO  TO  3305 
IF<KFACT.EQ.2)G0  TO  8960 

C  OUTPUT  SUftlARY  OUTPUT  UITHOUT  KEULEGAN  COEFFICIENT  CONSIDERED 
URITE  ( IODEV.3320)  IGAGENt 1. I) , 1GAGENC2. 1) .SMEN, VAROU. 

&PMEN . PVARO . PODc . SDPO . l'ELTOU ,  SDPOIJT .  1.  IERCD.NP 
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3320  FORMAT ( 2X. 2A2. 4X. 2 ( 2 ( F7 . 3 . 5H  +✓-  ,F7.3.2X).1X).23X.I3.4X. 

8,12.  IX.  13) 

J  .3310 

C  OUTPUT  SUhWARY  OUTPUT  UITH  KEULEGAN  COEFFICIENT  CONSIDERED 
0960  WRITE  < I0DEV.8961 )  IGAGEN< 1. I) . IGAGEN(2. D.SMEN.VAROU.PttN.PVARO 
&PODE.SDPO.DELTOU.SDPOUT.PMENK. I, IERCD.NP 
8961  FORMAT (  2X, 2A2 » 4X. 2 ( 2 ( F? . 3 . 5H  +/-  .F7.3.2X) IX) . 1X.F7.3. 15X. 

8,13. 4X.  12.  IX.  13) 

J  .3310 

C  OUTPUT  STANDARD  OUTPUT 

3305  WRITE  ( IODEV.5013)  SOUT . VARO . SfEN . VAROU . URMSO . URMS . PCT, PMOUT. 
1PVAO.PMEN. 

1 PVARO , PODE . SDPO . DELTOU . SDPOUT 
IF(KFACT.EQ.2)URITE  (I0DEV.8972)  FACTOR, PMENK 
8972  F0RMAT(3X, 18HATTENUATI0N  FACTOR. 2X,F9.6,5X, 

8.20HADJUSTED  WAVE  HEIGHT.2X.F7.3) 

3310  IF ( .NQT.SENSU(3) )GO  TO  40 
C  SCALE  DATA  TO  BE  PLOTTED 
DO  8007  K9-1.NPTPLT 
L-NPTPLT-K9+1 
XB -FLOAT ( I DATA ( L ) ) 

GO  TO  (8003,8003,40.8005), IF3 
8003  L  .72 

J  .8006 
8805  N-ll 

L  .74 
6006  Z(L)«FXX 
8007  CONTINUE 
TT3-12. 

IF (NOUPD . GT. 24) TT3 -FLOAT (NOUPD)  Si . 

C  PLOT  DATA 

CALL  M0DE(7,TT3.8.5,SAf€) 

CALL  MODE < 2. TT3. SAME. SAME) 

CALL  M0DE(3,8.5,SAhE,SAME) 

CALL  SCAN(2P,DELT. -2.480) 

CALL  M0DE(-8,TT3.DX,XB) 

DX-DELT/DX 

CALL  SCAN ( DELT , Z . -NPTPLT , 840) 

CALL  DRAW ( DX, Z . NPTPLT ,841) 

CALL  MODE (5, SAME. SAME. 440.) 

CALL  NOTE (PEAKPS . XXX, 42 . -NPEAK ) 

CALL  NOTE ( TROFPS , YYY , 42 . -NTROF ) 

CALL  AXES(22.2,22HTIME  PROTOTYPE  SECONDS. 

125.3, 25HELE VAT I ONS  PROTOTYPE  FEET) 

C  PLOT  TITLE 

CALL  MODE (3, 9. 4, SAME. SAME) 

CALL  N0TE<2. .8.7.4HRUN  .4) 

CALL  NOTE (SAME .SAME, IRUNNO.0) 

CALL  NOTE ( SAME , SAME . 7H  GAGE  .7) 

CALL  NOTE (SAME, SAME. IGAGEN( 1 . 1) ,4) 

CALL  MODE(4. .2. . 167,8. ) 

CALL  NOTE(2. .9..KTI7LE.40) 

CALL  MODE(4,.l.. 067.0.) 

CALL  DRAU(0..0., 1.8000) 

J  .40 

6400  CONTINUE 
40  CONTINUE 
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IF( .N0T.SENSLK6) )GO  TO  7382 
C  OUTPUT  RECORDS  TO  OUTPUT  TAPE  FILE 
CALL  SELMT ( IUN1T2) 

N=NCH#2 

C  OUTPUT  TRANSDUCER  FLAG  RECORD 

CALL  URTMT< I FLAG <1 , 1 ) , N, . TRUE . , I ERR) 

L  99 
N«NR0DS*2 

C  OUTPUT  GAGE  NAME  RECORD 

CALL  URTMT( IGAGENC 1. 1) ,N, .TRUE. . I ERR) 

L  .99 

IF(KFACT.NE.2)G0  TO  8983 
C  OUTPUT  WEIGHTED  MEAN  DEPTH  RECORD 
CALL  URTMT(DEPTH.N. .TRUE. . IERR) 

L  .99 

C  OUTPUT  DISTANCES  FROM  WAVE  GENERATOR 
CALL  URTMT<DIST,N.. TRUE.. IERR) 

L  .99 
8983  L-=NCH*6 

C  OUTPUT  GAGE  COORDINATE  RECORD 

CALL  URTMT(COOR( 1. 1)  ,L, .TRUE.. IERR) 

L  .99 

C  OUTPUT  ANALYSIS  RESULT  RECORD 

CALL  WRTMT < RESULT ( 1 ) , N , .TRUE., IERR) 

L  .99 

CALL  UEFMTC. TRUE., IERR) 

CALL  UEFMT ( . TRUE . , IERR) 

7382  I F < SENSU ( 3 ) ) CALL  DRAU<0. 0,0, 9999) 
IF(SENSU(5) )G0  TO  8 
C  POSITION  TO  NEXT  INPUT  TAPE  FILE 
76  CALL  SELMT ( I UN IT) 

CALL  SFFMT  < . TRUE . . I ERR ) 

IRL*5 

CALL  RDMT ( I DATA  < 1 ) . I RL . . TRUE . . IERR ) 

IF < IERR. EG.' 100002) GO  TO  77 
IF< IERR.NE.0)GO  TO  7 
J  .S 
C 

C  LINKS 
C 

C  AVERAGE  CALIBRATION  VOLTAGE  READINGS 
2100  OCT  0 

DO  2901  LC1-1.11 
LC2eK+LCl-l 
i<K  1  =0 

J=L+IP0IT1(LC1) 

L  .2902 
XP-1. 

1F( (LCl .EQ. 1 ) .OR. <LCl .EQ. 1 1) )GO  TO  2901 

j=l+ipgit2<:lci) 

L  .2902 
XP-2. 

IF(LC1 .NE.6)G0  TO  2901 
J-L+20 
L  .2902 
XP-3. 

2901  BUF  <LC2)  «=FLOAT  <K1<  1 )  /'XP 
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J  .2100.2 
2902  OCT  0 

KK 1 -KK 1 + I CALV <  J . NC ALV) 

J  .2902,2 

C  INTERPRET  GAGE  FLAGS 
2S  OCT  0 

I 1-1FLAG(2.  D/10000 

I F ( (1 1 . EQ . 1 ) . OR . ( 1 1 . EQ . 3) ) NCOEF -NCOEF+1 
I F  < 1 1 . GT . 2 1 NCAL V-NCALV+ 1 
J  .26,2 

C  DETERMINE  IF  DISK  I/O  ERROR  OCCURED 
8999  OCT  0 

IF( IERR.NE.01GG  TO  7004 
J  3999  p 2 

C  DETERMINE  IF 'MAGNETIC  TAPE  I/O  ERROR  OCCURED 

99  OCT  0 

IF< IERR.NE.01GG  TO  7 
J  99,2 

C  LINEAR' 8.  QUADRATIC  FITS 
72  OCT  0 

FXX-  ( ( S2  <  3 )  *XB+S2  <  2 ) 1  *XB+S2  <  1)  1  xMERT-tCUL 
J  .72,2 

C  SPLINT  ROUTINE  FOR  SPLINE  FITTING  PROCEDURE 
74  OCT  0 
MM-0 
K=-l 
XP=XB 

1F0^-X(  11)100. 170, 110 

100  MM— 1 

xp-xai 

J  .170 

110  IF(XP-X(N) 1130, 150, 140 
120  IFO-P-XCIO  1160. 170.130 
130  K-K+l 

J  .120 
140  MM— 1 
XP-XCN) 

150  K-N 

160  K-K-l 

170  HTl=XP-X(K) 

HT2-XP-X(K+1) 

PR0D=HTl*HT2 

DX»X(I<+1)-X(K) 

DELY»(2Y<K+1)-2Y<K))/DX 
S3°(S2(K+1)-S2(K) l/DX 
FPPXX-S2 <  K 1 +HT 1 *S3 
DELSQS=(S2(K)  +S2  (K+l  )-tf  PPXX)  /6 . 

FXX-2Y  ( IO  +HT  l*DELY+PROD*DELSQS 
FPXX=DELY+(HTl+HT2)*DELSQS+PR0D*S3/6. 
IF(MM.EQ.0)GO  TO  1S0 
FXX-FXX+FPXX*  <  XB-XP 1 
180  CONTINUE 

FXX-FXX*VERT+C(JL 
J  .74.2 

C  COMPUTE  MEAN  AND  VARIANCE 
3001  OCT  0 
sum*©. 


B17 
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SUM1-0. 

DO  3032  NN-l.NPK 
SUMM=SUm+UH(NN) 

3002  SUM1«SUM1+UH<NN)*UH(NN) 
PMEN=SUMM/FLOAT<NPK) 

PVAR - ( SUM1  -S(JMM*PMEN  > /FLOAT (NPK > 
J  .3001,2 
C 

6033  OCT  0 

YYY<I11)«FXX 

XB-FL0AT<IXXX(I11)? 

J  .8933.2 
??  CALL  MONOUT 
END 


PROGRAM  SIZE  *  ' 44226 


■  vvvvvv 


1.  Wave  data  are  recorded  m  data  files  on  magnetic  tape  by  ADACS  data 
acquisition  programs.  These  data  files  on  magnetic  tape  are  used  as  input  to 
the  wave  analysis  program. 

2.  Two  general  file  formats  are  used  by  the  wave  analysis  program--a 
calibration  file  and  a  noncalibration  file.  Several  records  containing  gage 
identification  and  data  scaling  factors  are  included  in  a  calibration  file 
that  are  excluded  from  a  noncalibration  file.  Since  gage  identification  and 
data  scaling  factors  remain  constant  for  a  series  of  data  tests,  there  is  no 
need  to  record  this  information  in  every  data  file.  Thus  this  information, 
which  is  recorded  exclusively  in  a  calibration  file,  is  applicable  to  all  sub' 
sequent  data  files  until  updated  by  another  calibration  file.  That  is  to  say, 
before  analyzing  a  particular  noncalibration  file,  the  wave  analysis  program 
must  obtain,  for  the  file,  the  correct  calibration  coefficients  which  are  con¬ 
tained  in  the  last  preceding  calibration  file.  It  should  be  noted  that  float¬ 
ing  point  variables  require  two  words  for  storage  and  two  alpha  characters  are 
stored  in  one  word. 

3.  All  data  files  are  separated  by  an  end-of-file  (EOF)  mark  with  the 
individual  file  and  records  having  the  following  format: 

a.  File  Identification  Record.  This  record,  which  is  included  in 
both  calibration  and  noncalibration  files,  is  used  simply  for 
file  identification.  It  contains  five  words  with  the  following 
definitions : 


Word  Data  Type 


1  Alpha 

2  Alpha 

3  Alpha 

4  Integer 

5  Integer 


Definition 


File  name  ASCII  Characters  (3A2) 

File  name  ASCII  Characters 

File  name  ASCII  Characters 

File  Type  Designator  Always  =  3 

Record  Length  of  Header  Record  #1 


b.  Header  Record  #1 .  This  record,  which  is  included  in  both  cali¬ 
bration  and  noncalibration  files,  contains  testing  identifica¬ 
tion  and  testing  parameters.  It  contains  20  words  with  the 
following  definitions: 


Data  Type 


Definition 


Word 

Data  Type 

Definition 

6 

Alpha 

Model  identification  (4A2) 

7 

Alpha 

Model  identification 

8 

Alpha 

Model  identification 

9 

Alpha 

Model  identification 

10 

Integer 

Horizontal  model  scale 

11 

Integer 

Vertical  model  scale 

12 

Integer 

Reference  mean  water  level  in  prototype 
feet  x  10 

13 

Integer 

Model  water  temperature  x  10 

14 

Integer 

Number,  N+M,  of  transducers  sampled  for 
this  test 

15 

Logical 

.TRUE,  if  this  is  a  calibration  file 

16 

Floating 

Model-generated  wave  height  in  prototype 

Point 

feet 

17 

Not  used  by  this  program 

19 

Not  used  by  this  program 

20 

Integer 

Record  length  of  Header  Record  #2  =  50 

Header 

Record  #2. 

This  record  contains  further  testing  param- 

eters . 

It  consists  of  50  words  with  the  following  definitions: 

Word 

Data  Type 

Definition 

1 

Floating 

Point 

Generated  wave  period  in  model  seconds 

2 

Not  used  by  this  program 

3 

Integer 

Number,  J,  of  data  samplings  per  gage  per 
wave  period 

4 

Integer 

Number,  K,  of  magnetic  tape  data  records  per 
wave  period 

5 

Integer 

Number,  L,  of  wave  periods  of  data  recorded 

6 

• 

• 

Not  used  by  this  program 

9 

• 

Not  used  by  this  program 

10 

Integer 

Number,  N,  of  wave  gages  sampled  for  this 
test 

11 

• 

• 

• 

• 

Not  used  by  this  program 

15  Not  used  by  this  program 

16  Integer  Transducer  arrangement  number 

17  Integer  Generated  wave  direction  in  degrees  from 

north 

18  Not  used  by  this  program 


28 


Not  used  by  this  program 
(Continued) 
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Definition 


t 


Word 

Data  Type 

29 

Logical 

30 

Floating 

31 

Point 

32 

Floating 

33 

Point 

. 

• 

50 

Transducer  Flag  ] 

.TRUE,  if  weighted  mean  depth  record  and 
gage  distance  record  recorded 
Wave  generator  length  in  model  feet 


Wave  machine  length  in  model  feet 
Not  used  by  this  program 


Not  used  by  this  program 


calibration  file.  Its  record  length  is  two  times  the  fourteenth 
element  (N+M)  in  Header  Record  #1.  This  record  contains  two 
integer  flags  for  each  computer  channel  read  during  a  data  scan 
including  sequential  wave  gages  that  are  scanned  but  not  used. 
The  flags  are  set  up  as  follows: 


(1)  FLAG1  =  Transducer  type  +  (computer  channel  number)  x  100  . 

Transducer  types  vary  from  0  to  99  with  1  identify¬ 
ing  wave  gages.  For  any  other  value  besides  1,  this 
analysis  program  does  not  consider  the  gage  data. 


(2)  FLAG2 


2  Gage  not  used 

0  Linear  calibration  fitting  procedure 
1  Quadratic  calibration  fitting  procedure 
9  Spline  calibration  fitting  procedure 


+  ((number  of  calibration  voltage  readings  per 
channel)  x  100 J 


+ 


=  0  No  entry  for  gage  in  calibration  voltage 
record  or  calibration  coefficient  record 

=  10,000  entry  for  gage  in  calibration  coeffi¬ 
cient  record  only 

=  20,000  entry  for  gage  in  calibration  voltage 
record  only 

=  30,000  entry  for  gage  in  both  calibration 
voltage  record  and  calibration  coefficient 
record 


The  two  (N+M)  words  in  this  record  have  the  following 
definitions : 


Word  Data  Type 


Definition 


1 

2 

3 


Integer  Flag  1  for  first  wave  gage 

Integer  Flag  2  for  first  wave  gage 

Integer  Flag  1  for  second  wave  gage 


(Continued) 


Integer  Flag  2  for  second  wave  gage 


Integer 

Integer 

Integer 

Integer 


2(N+M)-1  Integer 
2(N+M)  Integer 


Flag  1  for  last  (N  )  wave  gage 
Flag  2  for  last  (N  )  wave  gage 
Flag  1  for  first  of  other  type  gages 
Flag  2  for  first  of  other  type  gages 


Flag  1  for  last  (M  of  other  type  gages 
Flag  2  for  last  (M  )  of  other  type  gages 


Calibration  Coefficient  Record.  The  length  of  this  record  is 
obtained  by  multiplying  22  times  the  number  of  gages  indicated 
as  having  entries  in  this  record  by  Flag  2  in  the  Transducer 
Flag  Record.  This  record  is  included  only  in  a  calibration 
file.  Eleven  floating  point  calibration  coefficients  exist  in 
this  record  for  each  gage  having  calibration  coefficients.  Th 
reader  is  reminded  that  a  floating  point  number  requires  two 
words  in  memory.  The  words  in  this  record  have  the  following 
definitions : 


Data  Type 


Definition 


1 

Floating 

Mathematical 

fit 

coefficients 

for 

first 

Point 

wave  gage 

2 

• 

• 

3 

Floating 

Point 

4 

• 

• 

• 

• 

• 

• 

• 

23 

• 

Floating 

Mathematical 

fit 

coefficients 

for 

second 

Point 

wave  gage 

24 

• 

• 

25 

Floating 

Point 

26 

• 

• 

: 

• 

45 

Floating 

Mathematical 

fit 

coefficients 

for 

third 

Point 

wave  gage 

46 

• 

• 

22N-21 

Floating 

Mathematical 

fit 

coefficients 

for 

last 

Point 

(N1"*1)  wave 

gage 

(Continued) 

22N-20 


Data  Type 


Definition 


22N+1 


22N+2 


Floating  Mathematical  fit  coefficients  for  first 
Point  of  other  type  gages 


22(N+M)-21  Floating  Mathematical  fit  coefficients  for  last 


22(N+M) 


Point 


(M  )  of  other  type  gages 


f.  Gage  Name  Record.  Record  length  is  twice  the  fourteenth  ele¬ 
ment  (N+M)  of  Header  Record  #1.  This  record  is  included  only  in 
a  calibration  file  and  contains  a  four  character  name  or  number 
for  each  gage  scanned  during  data  acquisition.  The  2(N+M)  words 
of  this  record  have  the  following  definitions: 


Data  Type 


Definition 


$1 


Alpha  First  wave  gage  name  (2A2) 


Second  wave  gage  name 


2(N+M)-1 
2 (N+M) 


Last  (N  )  wave  gage  name 
First  of  other  type  gage  name 


Last  (M  )  of  other  type  gage  name 


Weighted  Mean  Depth  Record.  Record  length  is  two  times  the 
tenth  (N)  element  in  Header  Record  it 2.  This  record  is  included 
only  in  a  calibration  file  if  the  twenty-ninth  element  of  Header 
Record  #2  is  .TRUE.  It  contains  an  entry  for  each  wave  gage 
representing  the  average  model  depth  in  feet  between  the  gage 
and  the  wave  generator.  The  2N  words  in  this  record  have  the 
following  definitions: 


Data  Type 
Floating  Point 


Floating  Point 


_ Description _ 

Weighted  mean  depth  for  first 
wave  gage 

Weighted  mean  depth  for  second 
wave  gage 


(Continued) 
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Word 


2N-1  Floating  Point  Weighted  mean  depth  for  last 

(N^)  wave  gage 
2N 

h.  Gage  Distance  Record.  Record  length  is  two  times  the  tenth  (N) 
element  in  Header  Record  #2.  This  record  is  included  only  in  a 
calibration  file  if  the  twenty-ninth  element  of  Header  Record 
#2  is  .TRUE.  It  contains  an  entry  for  each  wave  gage  represent¬ 
ing  the  distance  between  each  gage  and  the  wave  generator  in 
prototype  feet.  The  2N  words  in  this  record  have  the  following 
definitions : 


Data  Type 


Description 


Floating  Point 


Distance  for  first  wave  gage 


Floating  Point 


Distance  for  second  wave  gage 


2N-1  Floating  Point  Distance  for  last  (Nttx)  wave  gage 


i.  Gage  Coordinate  Record.  Record  length  is  six  times  the  four¬ 
teenth  element  (N+M)  of  Header  Record  // 1.  This  record  is  in¬ 
cluded  only  in  a  calibration  file  and  contains  three  floating 
point  coordinates  (X,  Y,  and  prototype  water  depth)  for  each 
gage  scanned  during  data  acquisition.  The  6 (N+M)  words  in  this 
record  have  the  following  definitions: 


6N-5 

Floating  Point 

X  coordinate  for  last 

6N-4 

• 

• 

6N-3 

• 

Y  coordinate  for  last 

6N-2 

• 

• 

6N-1 

• 

Water  depth  for  last 

6N 

• 

• 

(Continued) 
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(Nth) 

(Nth) 

(Nth) 


wave  gage 
wave  gage 
wave  gage 


Data  Type 


Description 


Floating  Point  X  coordinate  for  first  of  other  type 
gages 


6(N+M)-1 


6(N+M) 


Water  depth  for  last  (M  )  of  other 
type  gages 


Reference  Potentiometer  Calibration  Coefficient  Record.  Record 
length  is  twice  the  tenth  element  (N)  of  Header  Record  #2.  This 
record  is  included  only  in  a  calibration  file  and  contains  the 
slope  coefficient  for  the  reference  potentiometer  of  each  wave 
gage  calibration  stand.  The  2N  words  in  this  record  have  the 
following  definitions: 


Data  Type 


Definition 


Floating  Point  Potentiometer  coefficient  for  first 


wave  gage 


Potentiometer  coefficient  for  second 


wave  gage 


Potentiometer  coefficient  for  last  (N  ) 
wave  gage 


Calibration  Voltage  Record.  The  length  of  this  record  is  deter¬ 
mined  by  multiplying  U2  times  the  number  of  gages  indicated  as 
having  entri  s  in  this  record  by  Flag  2  in  the  Transducer  Flag 
Record.  This  record  is  included  only  in  a  calibration  file  and 
contains  calibration  voltage  readings  for  wave  gages  calibrated 
under  computer  control  prior  to  data  acquisition.  The  words  in 
this  record  have  the  following  definitions: 


Calibration 


Data  T\ 


Steps 


Integer 


Zero 

+A 

+B 

+C 

+D 

+E 

+D 

+C 

+B 

(Continued) 


First  wave  gage  sensor 


Word 


Calibration 
Data  Type  Steps 


Gage 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 


Integer 


Integer 


+A 

Zero 

-A 

-B 

-C 

-D 

-E 

-D 

-C 

-B 

-A 

Zero 

Zero 

+A 

+B 

+C 

+D 

+E 

+D 

+C 

+B 

+A 

Zero 

-A 

-B 

-C 

-D 

-E 

-D 

-D 

-E 

-D 

-C 

-B 

-A 

Zero 

Zero 

+A 


First  wave  gage 
potentiometer 


Second  wave  gage  sensor 


45N-44  Integer  Zero  Last  (N  )  wave  gage  sensor 


T  t  h  S 


45N-23  Integer  Zero 


Last  (N^)  wave  gage 
potentiometer 


(Continued) 


Yi 


1 

P 

f 


Word 

Data  Type 

Calibration 

Steps 

Gage 

45N+1 

Integer 

Zero 

First  other  type 

gage 

sensor 

45(N+M)-23 

Integer 

Zero 

Last  (M*"*1)  other 
potentiometer 

• 

type 

gage 

45 (N+M) 

Integer 

Zero 

* 

Last  (M1"*1)  other 
potentiometer 

type 

gage 

1.  Data  Records.  The  length  of  this  record  can  be  determined  by 
multiplying  the  third  element,  J,  of  Header  Record  #2  by  the 
14th  element,  N+M,  of  Header  Record  #1  and  dividing  this  product 
by  the  fourth  element,  K,  of  Header  Record  #2.  These  data 
records  are  included  in  both  calibration  and  noncalibration 
files.  Since  the  number  of  data  records  in  these  files  varies 
and  depends  upon  the  number  of  wave  periods,  L,  acquired  during 
each  wave  test,  the  number  of  data  records  can  be  determined  by 
multiplying  the  fourth,  K,  and  fifth,  L,  elements  of  Header 
Record  # 2.  The  (J/K)(N+M)  words  in  this  record  have  the 
following  definitions: 

Word  Data  Type 

1  Scaled  fraction 

2 


_ _ Description  _ _ 

Sample  from  first  scan  of  first 
wave  gage 

Sample  from  first  scan  of  second 
wave  gage 


N 

N+l 

N+2 


Scaled  fraction 


Scaled  fraction 
Scaled  fraction 


Sample  from  first  scan  of  last 
(N*'*1)  wave  gage 

Sample  from  first  scan  of  first 
of  other  type  gage 
Sample  from  first  scan  of  second 
of  other  type  gage 


N+M 

(N+M)+l 

(H+M)+2 


Scaled  fraction  Sample  from  first  scan  of  last  of 

other  type  gage 

Scaled  fraction  Sample  from  second  scan  of  first 

wave  gage 

Scaled  fraction  Sample  from  second  scan  of  second 

wave  gage 
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1.  An  output  tape  may  be  generated  by  the  wave  analysis  program  con¬ 
taining  wave  analysis  results,  model  and  wave  test  identification,  wave  test 
parameters,  and  gage  identification.  This  output  tape  is  generated  to  provide 
permanent  files  of  reduced  wave  data  for  each  model  wave  test.  These  reduced 
data  are  used  to  produce  computer-generated  tables  or  plots  for  comparing 
results  of  various  wave  tests  conducted  in  a  particular  wave  model.  The  out¬ 
put  tape  is  simply  a  copy  of  pertinent  records  contained  on  the  input  tape 
with  the  addition  of  one  record  containing  the  results  of  the  wave  analysis 
program.  There  is  no  differentiation  between  a  calibration  file  and  a  non¬ 
calibration  file  on  the  output  tape.  That  is  to  say,  those  files  which  were 
excluded  on  noncalibration  input  tape  files  are  included  on  each  output  tape 
file. 

2.  Header  Record  #1,  Header  Record  it 2,  Transducer  Flag  Record,  Gage 
Name  Record,  Weighted  Mean  Depth  Record,  Gage  Distance  Record,  and  Gage  Co¬ 
ordinate  Record  are  included  in  every  output  tape  file.  These  records  are 
simply  copied  directly  from  input  tape  to  the  output  tape  using  the  same  for¬ 
mats  that  are  defined  in  Appendix  C  for  these  records.  The  last  record  in¬ 
cluded  in  the  output  tape  file  is  the  Analysis  Result  Record.  This  record 
length  is  twice  the  tenth  element,  N,  of  Header  Record  #2.  This  record  con¬ 
tains  the  average  highest  one-third  of  the  wave  heights  for  each  wave  gage. 
There  is  one  floating  point  entry  in  this  record  for  each  wave  gage.  The  2N 
words  in  this  record  have  the  following  definitions: 

Word  Date  Type  _ _  Definition _ 

1  Floating  Point  Wave  height  for  first  wave  gage 

2  • 

3  Floating  Point  Wave  height  for  second  wave  gage 

4 

•  •  • 

•  •  • 

•  •  • 

2N-1  Floating  Point  Wave  height  for  last  (N  )  wave  gage 

2N 
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C  TIDE  DATA  ANALYSIS 


FORTRAN  COMPILER  REV.  LEV.  J00 


C 

C  SENSE  SUITCH  A  DELET  ORIGINAL  DATA  PLOT 

C  SENSE  SWITCH  B  DELET  ORIGINAL  &  COMPARISON  DATA  PLOT  WITH  MEAN 
C  ADDED  TO  BOTH  DATA  SERIES 

C  SENSE  SUITCH  C  DELET  ORIGINAL  8.  COMPARISON  DATA  PLOT  WITH  MEAN 
C  REMOVED  FROM  BOTH  DATA  SERIES 

C  SENSE  SUITCH  D  DELET  RESIDUAL  PLOT 
C  SENSE  SUITCH  E  DELET  POUER  SPECTRAL  DENSITY  8.  PLOT 
C  SENSE  SUITCH  F  SKIP  CHANNELS 
C  SENSE  SUITCH  G  LIST  DATA 
C 
C 

C  PROGRAM  LIMITS 

C  1680  DATA  POINTS  PER  GAGE  MAXIMUM 

C  20  LEAST  SQUARES  HARMONIC  ANALYSIS  COMPONENTS  MAXIMUM 

C  98  DATA  CHANNELS  MAXIMUM  48  TIDE  GAGES  MAXIMJM 

C  DATA  RECORD  LENGTH  3362  WORDS  MAXIMUM 

P 

C 

DIMENSION  PERHLDC20) 

D I MENS  ION  ZP ( 3 ) . ZR ( 3 ) . IDATA ( 3362 ) . IBUF ( 3362 ) , 

8,IZAP<2)  ,PERIOD(20)  ,OMEGA<20) ,  AMP(20) ,PHS(20) . 

8.INDEXC 10) ,ZAP(2) .PSDC( 1681 ) ,S(5) .AAA(3,4) . IFLAG(2.98) . 

SiXCOR  ( 98 ) ,  YCOR  ( 98 ) . ZCQR ( 98 ) , 

8.IGAGEN(2.98)  -COEF (11,50)  ,NSUL<48) 

DIMENSION  DATA<1680),LABELC10), 

8.Z(1681).CH0LD(41.41) 

EQU I VALENCE  ( I BUF ( 1 ) . PSDC  < 1 ) . CHOLD (1,1)) 

EQU I VALENCE  ( TX. I ZAP  < 1 ) ) . < I DATA (1 ) , DATA ( 1 ) . Z ( 1 ) ) 

SCALED  FRACTION  I VOLT 
LOGICAL  1ADJCH, 1REFD 
LOGICAL  ICAL.SENSU, ITAPER. IGTYPE 
DATA  TPP/0./ 

DATA  PI/3. 1415927/. MODULO/"  100/.SAM£/9999./,Tl/0./,T2/0./ 
DATA  LMASK3/' 74 17/ 

DATA  DX/0./.LMASK/' 10421/.LMASK2/' 177777/ 

DATA  IRESP/1HY/ 

DATA  I ODE V/ 1 6/ , RADDEG/57 . 2957795/ . TUOPI/6. 283 1 853/ 

0  TYPE  6501 

6501  FORMAT(/. 13HF0RCE  SHIFT-?./) 

ACCEPT  901.FSHIFT 
1 TAPER -.FALSE. 

TYPE  6092 

6092  FORMAT (/. 16HTAPER  RESIDUAL  ?./) 

ACCEPT  2. I 
2  FORMAT(Al) 

I F  < I . EQ . I RESP ) I TAPER - . TRUE . 

TYPE  6202 

6202  FORMAT (/.24HCAL IBRATION  FIT  OPTION-?./) 

ACCEPT  6252.N0PTN 
TYPE  8111 

8111  FORMAT </,22H*C0MPS  TO  ANALIZE  FOR?./) 

ACCEPT  6252, INCANL 
6252  F0RMAT<217) 
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C  TIDE  LATA  ANALYSIS 


DO  S391  I-l. INCANL 
TYPE  6254 

6254  FORMAT (/.10HCOMP  PER-?./) 

6391  ACCEPT  901 .PERHLD( 1 ) 

IGTYPE  -.FALSE. 

TYPE  6210 

6210  FORMAT!/. 17HLIST  UPDATE  CALS?./) 

ACCEPT  2# 1 

IF ( I .EQ. IRESP) IGTYPE-. TRUE. 

IREFD-.TRUE. 

TYPE  6306 

6306  FORMAT (/.27HREFERENCE  DATA  TO  BE  INPUT?./) 
ACCEPT  2.1 

IF ( I .ME . IRESP)GO  TO  6307 
1REFD-. FALSE. 

TYPE  6308 

6308  FORMAT (/.32HNUM8ER  OF  REFERENCE  COMPONENTS-?./) 
ACCEPT  6252, NC ANAL 

6307  TYPE  101 

101  FORMAT (/.22HHOU  OFTEN  UPDATE  CAL  ?./) 

ACCEPT  6252.KALUPD 
TYPE  6392 

6392  FORMAT (/.7HLABEL-?./) 

ACCEPT  6393. (LABEL <I). 1-1.10) 

6393  FORMAT  <  10A2) 

TYPE  2092 

2092  FORMAT (/, 14H*REC  TO  SKIP-?./) 

ACCEPT  6252.NRSK IP 
TYPE  6297 

6297  FORMAT (/. 11 HTAPE  UNIT-?./) 

ACCEPT  6298,11 

CALL  SELMT(Il) 

6298  FORMAT (08) 

TYPE  900 

900  FORMAT (/, 2 1 HAVERAGE  MODEL  DEPTH-?./) 

ACCEPT  901. DEPTH 

901  FORMAT(F10.5) 

IRL-0 

TYPE  3 

3  FORMAT (/. 10HF I LE  NAME?,/) 

ACCEPT  3311. (IDATA(I). 1-1.3) 

3311  FORMAT (3A2) 

C  POSITION  TAPE  TO  FILE 

CALL  POSMT (I DATA ( 1 ) . 3 . I RL . . F ALSE . , I ERR ) 

IF ( IERR.EQ. ' 100004)GO  TO  5 
IF ( IERR.EQ.0)GO  TO  5 
7  TYPE  6 

6  FORMAT ( 1H?) 

J  .8 
5  IRL-20 

C  READ  FIRST  HEADER 

78  CALL  RDMTC IDATA<51). IRL. .TRUE., IERR) 

1CAL-.TRUE. 

IF ( IDATA(65) .GE.0) ICAL-. FALSE. 
lc ,IERR.NE.0)GO  TO  7 
C  READ  SECOND  HEADER 

CAi  I.  RDMT(IDATA<  101 ) .  IDATAC70) . . TRUE..  IERR) 
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IF(1ERR.NE.0)GO  TO  7 
TYPE  6299 

6299  FORMAT C/.24H*  OF  CYCLES  TO  ANALYSE-?./') 

ACCEPT  6252. I DATA ( 105) 

C  BREAK  OUT  DATE  AND  TIME 

I DATA C 1 ) ■ I DATA ( 54 ) / 1 000 
IDATAC2) -IDATA (54)- ( IDATAC 1)*1000) 

CUL -FLO AT  < I DATA  C  62 ) ) ✓ 1 00 . 

TEMPA -FLOAT (I DATA ( 63 ) ) ✓ 1 0 . 

DATAC3) -DATAC51) 

IHORZ- IDATA (60) 

IVERT- IDATA (61) 

UAVPD  -DATA  ( 3 ) /3600 . 

C  COMPUTE  MODEL  SCALES  TO  CONVERT  TO  PROTOTYPE 
H0R2 -FLOAT ( IH0R2) 

VERT-FLOAT ( I VERT) 

Tl-VERT/'HORZ 

SUL -5. 1 2*DATA ( 3 ) *DATA  <  3 ) 

C-SUL 

902  ARG-6 . 283#DEPTH/SUL 
SUA-C*TANH(ARG) 

I F ( ABS ( SUA-SUL ) -0 . 00 1 *SUL ) 903 , 903 . 904 
904  SUL -0 . 5* ( SUA+SUL ) 

J  .902 

903  THARG-TANH(ARG) 

ARGU-ARG*T1 
THARGU-TANH ( ARGU) 

T I ME -SORT ( HORZ*THARG/THARGU) 

VERT-VERT/12. 

DATAC4) -DATA  ( 3 ) *T I ME/3600 . 

C  MOVE  (REAL)  TIDE  HEIGHT 
1DATAC9) -1DATAC66) 

IDATAC 10) -IDATA(67) 

DATA(6)-DATA(5)*VERT 
C  BREAK  OUT  TESTING  PARAMETERS 
NCH-IDATA<64) 

NRECPD- IDATAC 104) 

IGAGAR-1DATAC 1 16) 

ISPUPD- IDATAC 103) 

NOUPD- IDATAC 105) 

NSPTST- I DATA ( 1 03 ) * I DATA ( 1 05 ) 

1RLD-1DATAC 1 03  )*NCH/' IDATAC  104) 

I RUNO - I DATA ( 53 ) 

C  OUTPUT  TEST  HEADER  INFORMATION 

URITE  (I0DEV.15)  I RUNO. IDATAC2). IDATAC 1) . IDATAC55) . IDATAC56) . 
&IDATAC57) . IDATAC5B) . IDATAC59) , IGAGAR. IHORZ. IVERT, TEMPA. 
tDATA C  3 ) . DATA  C  4) . NCH . I SPUPD . NOUPD . NRECPD 
15  FORMATC . 13X.36H***  TIDE  RECORD  ANALYSIS  PROGRAM  **». 
&/V/.9H  TEST  RUN.1X. I5.3X. I3.14H  TH  DAY  OF  197.I1.4H  AT  . 

&I4.6H  HOURS./'. 1 1H  MODEL  PLAN. 1X.4A2.3X.9HGAGE  ARR..1X.I4. 
IS.VM  HOR.  SCALE  1-.I3.15H  VER.  SCALE  1-. 

4I3.2X. 11HUATER  TEMP. . 1X.F5.2.//'. 16X.5HM0DEL. 1 3X. 9HPR0T0TYPE . 
K/.7H  PERI0D.7X.F7. 1.8H  SECONDS. 5X. F7. 1.6H  HOURS./. 

II3.13H  CHANNELS  OF  . 13. 15H  SAMP/PER.  FOR  .I3.8H  PERIODS. 

1/.11H  *REC./PER., IX. 14) 

IF ( IGTYPE)URITE  CIODEV.106) 

106  FORMATC /'/'/’/'. 20H  UPDATE  CALIBRATIONS. 5X.9H INTERCEPT. 6X.5HSL0PE) 
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C  TIDE  DATA  ANALYSIS 


1RL-NCH*2 

C  READ  TRANSDUCER  FLAGS 

CALL  RDIiT(  IFLAG<  1.1).  IRL.  .TRUE. .  I  ERR) 

IF( IERR.NE.0)GO  TO  7 

IF-5 

Jl-0 

NRE-0 

NCOEF-0 

C  COMPUTE  *  OF  CAL  COEFFICIENTS  &  •  OF  REFERENCE  ELEVATIONS 
DO  6291  I-l.NCH 

I2-IFLAG<1. 1)-<<IFLAG(1.I)/'100)*100) 

I1-IFLAG(2,I)/'10000 
1F<12.EQ.2)NRE-NRE+1 
IFC ( 1 1 .EQ. 1) .OR. ( 1 1 .EQ.3) )NCOEF-NCOEF+l 
IF(I2.NE.3)G0  TO  6291 
NRE-NRE+1 
Jl-Jl+1 
1F-IF+1 
INDEX(J1)-I 
1NDEX<1F)-NRE 
6291  CONTINUE 
C  READ  CAL  COEF 

IRL-NC0EF*22 

CALL  RDMT(COEF( 1, 1) . IRL. .TRUE. . I ERR) 

IF ( IERR.NE.0)GO  TO  7 
C  READ  GAGE  NAMES 
IRL«NCH*2 

CALL  RDMT(IGAGEN< 1. 1) . IRL. .TRUE. » IERR) 

IF< IERR.NE.0)GO  TO  7 
C  READ  GAGE  COOR. 

IRL«NCH*6 

CALL  RDMTUDATA(l).  IRL.  .TRUE..  IERR) 

IF< IERR.NE.0)GO  TO  7 
1-1 

DO  25  11-1.288.3 
XCOR ( I ) -DATA (II) 

YCOR< I) -DATA < 1 1+1) 

2C0R( I) -DATAC 1 1+2) 

25  1-1+1 

C  SKIP  2  RECORDS  POT  COEF  &  CAL  VOLTAGES  NOT  USED  BY  THIS  PROGRAM 
CALL  SRFMTC. TRUE.. IERR) 

CALL  SRFMT  ( . TRUE ..IERR) 

C  READ  REFERENCE  ELEVATIONS 

CALL  RDMT(NSULCl) >NRE. .TRUE. . IERR) 

IF(IERR.NE.0)GO  TO  7 
16  IADJCH-.FALSE. 

C  SORT  REFERENCE  TIDE  GAGES  IN  ORDER  OF  DISPLACEMENT 
NPER  »NOUPD#NRECPD 
DO  41  1-1.3 
DO  42  Kl-1.3 
NN-INDEXU) 

NP-INDEX(Kl) 

IF(ZCOR(NN) .GE.ZCOR(NP))GO  TO  42 
KS-lNDEX(Kl) 

INDEX(Kl) -INDEX( I) 

INDEXC I) -K5 
K5-INDEX(Kl+5) 
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lNDEX(Kl+5)“INDEX(I-f5) 

1NDEX(1+5)“K5 
42  CONTINUE 
41  CONTINUE 

ICl-lNDEX(l) 

IC2-INDEX(2) 

IC3-INDEX(3) 

PRANGE-ZCOR(ICl) 

Z£R0-ZC0R(IC2) 

ORANGE -ZCOR( IC3) 

NCCRC-0 

NSULRC-0 

NN-INDEX(6) 

ASUL “FLOAT (NSUL (NN) ) 

NN-INDEX(?) 

BSUL -FLOAT ( NSUL (NN)) 

NN-INDEX(8) 

CSUL -FLOAT (NSUL (NN) ) 

DELT -UAVPD/FLOAT ( lSPUPD)*TIf1E 
C  SKIP  DATA  RECORDS  IF  DESIRED 
IF (NRSKIP.EQ.0)GO  TO  2093 
DO  2091  IF-l.NRSKIP 
2091  CALL  SRFMT ( . TRUE . . IERR) 

2093  DO  40  I-l.NCH 

C  TYPE  PRESENT  GAGE  NUfBER  TO  TELETYPE 
TYPE  6281.1 

6201  FORMATS.  7HGAGE*-  .15) 

C  BREAK  PRESENT  CHANNEL  FLAGS 
IF-IFLAG(l.I) 

IF1-IF/100 
1F3-IF-IF1*100 
IF-IFLAG(2, I ) 

IF2-IF/10000 

IF( ( IF2.EQ. 1) .OR. (1F2.EQ.3) )NCCRC«NCCRC+1 
I F ( ( I F3 . EQ . 2 ) . OR . ( I F3 . EQ . 3 ) ) NSULRC -NSULRC+ I 
IF( ( IF3.EQ.2) .OR. ( IF3.EQ.4) .OR. ( IF3.EQ.5))G0  TO  46 
GO  TO  40 

46  IF2-1F-((IF/'10)*10) 

IF ( 1F2.EQ.2)G0  TO  40 
IF2“((IF/10)-<(IF/'100)*10))+2 
IF (SENSU(6) )G0  TO  40 
C  READ  &  DEMULTIPLEX  DATA  FOR  1  GAGE 
13-0 

ICALNM-1 
KALKT -KALUPD 
DO  43  Il-l.NPER 

6200  CALL  RDMT(I8UF(1).IRLD.. TRUE.. IERR) 

IF(IERR.NE.0)GO  TO  7 
K8-0 

DO  44  14-1. IRLD.NCH 
I3-I3+1 

C  COMPUTE  UPDATE  CALS 

GO  TO  (6001.6002.6241.6003) .NOPTN 
6001  ZP(l)-ZERO 

ZP(2) -PRANGE 
K5-IC2 
L  .6111 
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ZR< 1) -XB-BSUL 
K5-IC1 
L  .6111 
ZRC2)  ->®-ASUL 
K5-2 

L  .6261 
J  .6006 

6002  ZP Cl) -ZERO 
ZPC2) -ORANGE 
K5-IC2 

L  .6111 
ZR ( 1 ) -XB-BSUL 
K5-IC3 
L  .6111 
ZRC2) -XB-CSUL 
K5-2 

L  .6261 
J  .6006 

6241  ZPCD-PRANGE 
ZP (2) -ORANGE 
K5-IC1 
L  .6111 
ZRC 1) -XB-ASUL 
K5-IC3 
L  .6111 
ZR<2) -XB-CSUL 
K5-2 

L  .6261 
J  .6006 

6003  ZP( 1) -PRANGE 
ZPC2) -ZERO 
ZP(3) -ORANGE 
K5-IC1 

L  .6111 
ZRC 1) -XB-ASUL 
K5-IC2 
L  .6111 
ZRC2) -XB-BSUL 
K5-1C3 
L  .6111 
ZRC3) -XB-CSUL 
K5-3 

L  .6261 
6006  K5-I4-K8 


L  .6111 

1F<IF3.NE.2)G0  TO  6002 
REMOVE  ZERO  REFERENCE 

XB-XB-FLOATCNSUL CNSULRC) ) 

SCALE  DATA 

DATAC13) -<AAA<  i,S)4AAA(2.3>^*Sfeft+4dJL 

92  DATAC 13) -XB*C0EF(2»NCCRC)4C0EF<1»NCCRC) 
K8-K8+NCH 
CONTINUE 
IGTYPE-.FALSE. 
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C  SKIP  RECORDS  REVERSE  ON  MAG  TAPE 
DO  45  K2-1.NPER 
CALL  SRRMTC .TRUE. * IERR) 

IFC IERR.NE.0)GO  TO  7 
45  CONTINUE 
C 

C  SET  UP  FOR  PLOTS 
C 

KBARY-11 

T2-5.5 

KBARS-17 

Tl-8.5 

IFCNOUPD.LT. 10)GO  TO  6005 

T2-8.5 

KBARY-17 

Tl-25.5 

KBARS-51 

6005  CALL  M0DEC7.T1.T2.SAME) 

CALL  M0DEC2.T1, SAME. SAME) 

CALL  MODE ( 3. T2, SAME, SAME) 

T99-T2 

K2-2 

ZAP  C 1 ) “0 . 

ZAP (2) -DELT*FLOAT(NSPTST) 

CALL  SCAN ( ZAP , DELT , -K2 , 480 ) 

CALL  M0DE(-8.T1,DX,T2) 

DX-DELT/DX 

C 

C  SKIP  PLOT 

C  COMPUTE  ORIGINAL  DATA  PLOT  IF  DESIRED 
IF (SENSUC 1 ) )G0  TO  6059 
C 

CALL  SCAN (DELT. Z.-NSPTST. 840) 

CALL  DRAU ( DX, Z , NSPTST . 842 ) 

C  TEMP  GAGE 

I F  < I F3 . NE . 4) GO  TO  6242 

CALL  A>ES(20.2.20HTIME  PROTOTYPE  (HRS). 16.2. 

1 16HTEMP.  IN  DEGREES) 

J  .6243 

6242  IFC IF3.NE.5)G0  TO  7242 

CALL  AXESC20.2.20HTIME  PROTOTYPE  CHRS). 
18.2.8HVEL0CITY) 

J  .6243 

7242  CALL  AXESC20.2.20HTIME  PROTOTYPE  CHRS). 
8.26.2.26HTIDE  HEIGHT  PROTOTYPE  (FT)) 

6243  CALL  MODEC 10. FLOAT (LMASK) .SAME. SAME) 

CALL  FORMCKBARS. .5.KBARY. .5) 

L  .6251 

CALL  DRAUC0..0.. 1.8000) 

C  DO  NOT  DO  ANALYSIS  UNLESS  TIDAL  GAGE 
IF(IF3.NE.2)G0  TO  40 
C 

C  SET  UP  COMPONENTS  FOR  LEAST  SQUARES  HARMONIC  ANALYSIS 
6059  DO  6051  K2-1.INCANL 
6051  PERI0DCK2) •PERHLDCK2) 

INC-INCANL 
6091  J-I*2 
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C  PRINT  HEADER  FOR  PRESENT  GAGE 

URITE  (I0DEV.66)  (LABELCK2) .K2-1. 10) . IF1. IGAGAR. IRUNO 

66  F0RMATC31H1  ***  TIDE  ANALYSIS  PROGRAM  ***./.  IX.  10A2.//. 

&8H  CHANNEL. 2X. 13. 5X. 16HGAGE  ARRANGEMENT. 2X. I 4, 2X. 

&10HRUN  NUMBER. 14,/) 

UEE-ZCOR ( I ) +CUL 

URITE  ( 10DEV.67)  XCOR( I) ,YCOR( 1 ) .UEE, 1GAGENC 1. 1) . 

1IGAGENC2. I) 

67  FORMATC 15H  XY  COORDINATES, 2X.2CF9. 4.2X) ./. I5H  UATER  DEPTH  IS. 
&2X.F9.4./. 12H  GAGE  *  IS  ,2A2) 

ND1M-21 

C 

C  LIST  DATA  ? 

C 

1F(.N0T.SENSU(7))G0  TO  6247 

URITE  ( I0DEV.7881)  (Z(K2) ,K2-1 .NSPTST) 

7881  F0RMAT(8E15.7) 

C 

C  COMPUTE  8.  REMOVE  MEAN 
C 

6247  XBAR  -RMEAN ( Z , NSPTST) 

XVAR-RVAR ( Z . NSPTST, 1 . XBAR ) 

URITE  ( 10DEV.6295) 

6295  FORMATC//. 32H  STATISTICS  OF  THE  ORIGINAL  DATA) 

URITE  ( IODEV.6044)  XBAR.XVAR 
6044  FORMAT (8H  MEAN  -  .E14.7, 10X, 12HVAR1ANCE  -  .E14.7) 

DO  6009  J-l. NSPTST 

6009  Z(J)-Z(J)-XBAR 
C 

C  BEGIN  LEAST  SQUARES  HARMONIC  FIT 
C 

C  CALCULATE  OfEGA  -  TUOP I /PERIOD 
C 

6269  DO  6004  K2-1.INC 

6004  0MEGACK2) “TUOP 1 /PER 10DCK2) 

C  SET  UP  COUNTERS 
NC1-INC+1 
NC2-2*INC 
NCP-NC2+1 
C 

C  ZERO  MATRIX 
C 

DO  6010  K2-1.NDIM 
DO  6010  J-l.NDIM 

6010  CH0LDCK2. J) “0.0 
C  COMPUTE  MATRIX 

DO  6011  L-l. NSPTST 
TT -FLOAT (L- 1 ) *DELT 
DO  6012  J-l. INC 
A 1 -COS ( OMEGA ( J ) *TT) 

DO  6013  K2-J.INC 
A2 -COS ( OMEGA (K2)*TT) 

6013  CH0LD(J.K2)-CH0LD<J.K2>-HU*«2 
DO  6014  K2-NC1.NC2 
K0N1-K2-1NC 

A2 -SIN (OMEGA (KONl)*TT) 

6014  CH0LD(J,K2)-CH0LD(J,K2)+A1*A2 
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6012  CH0LD(J.NCP)-CH0LD(J.NCP)+A1*Z(L) 

DO  6015  J-NC1.NC2 
KONW-1NC 

Al-SIN(OMEGA(KONl)*TT) 

DO  6016  K2-J.NC2 

K0N1-K2-1NC 

A2 -S 1 N ( OMEGA  <  KON 1 ) *TT) 

6016  CHOLD (J.K2) -CHOLD (J.K2)+A1*A2 
6015  CHOLD( J.NCP) -CHOLD( J.NCP)+A1#Z(L) 

6011  CONTINUE 
C 

C  FILL  OUT  REMAINDER  OF  MATRIX  (  REAL  SYMMETRIC  ) 

C 

DO  6018  J-2.NC2 
Kl-J-1 

DO  6018  K2-1.K1 
6018  CH0LD(J.K2)-CH0LD(K2,J) 

C 

C  SOLVE  SYSTEM  OF  EQUATIONS  BY  GAUSS  JORDAN  METHOD 
C 

DO  6021  J-1.NC2 
Al-CHOLD(J.J) 

DO  6029  M-l.NCP 
6029  CHOLD(J.M)-CHOLD(J,M)/ttl 
DO  6023  NN-1.NC2 
IF(J.EQ.NN)GO  TO  6023 
A2-CHOLD(NN,J) 

DO  6024  Kl-l.NCP 

6024  CHOLD ( NN , K 1 ) -CHOLD (NN.K1) -A2*CH0LD  <  J , K 1 ) 

6023  CONTINUE 
6021  CONTINUE 
C 

C  END  INVERSE 
C 

DO  6031  Kl-l.INC 
KON1-K1+INC 

AIf>  (K 1 ) -SQRT(CHOLD (K 1 . NCP) **2+CH0LD (KON 1 . NCP) **2> 

6031  PHS (K 1 ) -ATAN2 ( CHOLD (KON l.NCP). CHOLD (K 1 .NCP)) *RADDEG 
C 

C  PRINT  FIT  RESULTS 
C 

WRITE  (IODEV.6072) 

6072  F0RMAT(/'//28H  ****  HARMONIC  ANALYSIS  ****/V6X. 

128H(HRS)  (FT)  (DEG)) 

WRITE  (IODEV.6032)  (PERIOD(KK) ,AMP(KK).PHS<KK) ,KK-1. INC) 

6032  FORMAT (5X. 29H  PERIOD  AMPLITUDE  PHASER IX. 1P3E12.4)) 
C 

C  END  HARMONIC  ANALYSIS 
C 

IF(IF3.EQ.4)G0  TO  6301 
C  IF  COMPAIRING  TO  ANALYSIS  GO  TO  6301 
1F( 1REFD)G0  TO  6301 
1NCHLD-INC 
INC-NCANAL 
PPHASE-PHS(l) 

C  READ  COMPARISON  DATA  FROM  CARDS 
DO  29  Kl-l.INC 
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29  READ ( 6 . 24) PER 1 OD ( K 1 ) . PHS ( K 1 ) . AMP  C K 1) 

READC6.241DEPH 
24  FORMAT<3F10.5) 

IF C IADJCH)GO  TO  6302 

C  COMPUTE  PHASE  SHIFT  SO  THAT  INPUT  COMPARISON  DATA  CONFORMS  TO  PHASE 
C  OF  ORIGINAL  DATA 

TPP - (PPHASE-PHS ( 1 ) ) *(PER 10D( 1 )  /"360 . ) 

IADJCH-.TRUE. 

6302  DO  6305  K1-1.1NC 

6305  OMEGACKl)«TUOPI/PERIOD(Kl) 

WRITE  (1ODEV.6402)  TPP.FSHIFT 

6402  FORMAT (/»33H  COMPUTED  SHIFT  IN  MODEL  HOURS  -  .E14.7.A 
&26H  FORCE  SHIFT  IN  DEGREES  -  .E14.7./0 
C  COMPUTE  COMPARISON  DATA 
6301  DO  6071  Kl-l.INC 
6071  PHS ( K 1 ) "PHS  C  K 1) /RADDEG 
DO  6036  Kl-1. NSPTST 

T2  .0 

TT “FLOAT (Kl-1) *DELT 
DO  6037  K2-1.INC 

6037  T2 -T2+AMP ( K2 ) *C0S  COMEGA ( K2 ) * (TT-TPP ) -PHS ( K2  > -FSH I FT) 

6036  PSDC(K1)-T2 

C 

C 

C  PLOT  ORIGINAL  &  GENERATED  DATA 

C 

C 

C  REMOVE  ft  AN  PLOT 
C 

IF(.N0T.SENSU(3))G0  TO  6321 
5010  DO  5000  Kl-1, NSPTST 
Z(K1)-Z(K1)+XBAR 

IF(.N0T.1R£FD)PSDC(K1)-PSDC(K1)+DEPH 

IF(IREFD)PSDC(K1)-PSDC(K1)+XBAR 

5000  CONTINUE 
Kl-0 

C  DELET  PLOT 

5001  IF(SENSU(2))GO  TO  6061 
C 

6321  CALL  SCAN (CELT. Z.NSPT5T. 840) 

CALL  SCAN (DELT.PSDC. -NSPTST. 840) 

IF(K1.EQ.0)GO  TO  7001 

IF(SENSU<7) )GO  TO  7001 

CALL  AXESC20.2.20HTIME  PROTOTYPE  (HRS), 

138.2.38HTIDE  HEIGHT  PROTOTYPE (FT)  MEAN  REMOVED) 

J  .7002 

7001  CALL  AXES(20.2.20HTIFE  PROTOTYPE  (HRS), 

C26.2.26MTIDE  HEIGHT  PROTOTYPE  (FT)) 

7002  CALL  MODE ( 10. FLOAT (LMASK) .SAME. SAFE) 

CALL  FORM(KBARS, .5.KBARY. .5) 

CALL  MODE ( 10. FL0AT(LMASK3) . SAFE. SAME) 

CALL  DRAU( DX»  PSDC » NSPTST . 842 ) 

CALL  MODE( 1 0. FLOAT (LMASK2) .SATE, SAME) 

CALL  DRAU(DX.Z. NSPTST, 842) 

L  .6251 

CALL  DRAU(0..0., 1.8000) 

C 
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IF(K1.EQ.0)GO  TO  6061 
J  .5010 

6061  IF( .NOT. IREFD) INC-INCHLD 
C 

DO  6038  Kl-1. NSPTST 
6038  Z<K1)-ZCK1)-PSDC(K1) 

C  COMPUTE  MEAN  &  VARIANCE  OF  RESIDUAL 
XBAR -RMEAN ( Z . NSPTST) 

XVAR -RVAR ( Z . NSP  TST , 1 , XBAR ) 

URITE  ( IODEV.6043) 

6043  FORMAT (///. 32H  STATISTICS  OF  THE  RESIDUAL  DATA) 
URITE  (IODEV.6044)  XBAR. XVAR 
C 

C  PLOT  RESIDUAL 

C 

C 

C  DELET  PLOT 

IF(SENSU(4) )GO  TO  7003 
C 

CALL  SCAN (DELT.Z. -NSPTST. 840) 

CALL  DRAU(DX,Z. NSPTST, 842) 

CALL  AXES(20.2,20HTIME  PROTOTYPE  (HRS). 16.3. 

1 16HRESIDUAL  OF  DATA) 

CALL  MODE ( 1 0. FLO AT (LMASK).S Aft. SAME) 

CALL  FORMCKBARS, .5.KBARY. .5) 

L  .6251 

CALL  DRAU(0.,0.. 1.8000) 

C 

C  DO  PSD  ON  DIFFERENCE 

C 

C 

C  DELET  PSDC 

7003  IF(SENSU(5))G0  TO  40 

C  SET  UP  FOR  POUER  SPECRAL  DENSITY 

C 

NSPHLD -NSPTST 

I F ( NSPTST . GT . 1 024) NSPTST- 1 024 
T1 -FLOAT (NSPTST) 

Ml-IFIX(AL0G(Tl)/AL0G(2.)+l.E-5) 

NPTS 1 -2**M1 

I F ( NSPTST . EQ . NPTS 1 ) GO  TO  6041 
Ml-Ml+1 

NPTS 1 -NPTS 1+NPTS1 

6041  MPTS -NPTS 1 -NSPTST 
MPTS2-MPTS/2 
NPT2 -NPTS 1/2 

ZBAR -RMEAN ( Z . NSPTST) 

VAR -RVAR(Z. NSPTST. 1 .ZBAR) 

DO  6042  I FT- 1. NSPTST 

6042  Z ( I FT) -Z ( I FT) -ZBAR 
C  TAPER  DATA  ? 

IF( .NOT. ITAPER)GO  TO  6093 
Kl-NSPTST+1 
K2-NSPTST/2 
DO  6049  IFT-1.K2 

T 1 -2 . 6667* ( S I N ( FLOAT ( I FT- 1 ) * ( P I /FLOAT ( NSPTST) ))) **2 
Z(IFT) *Z( IFT)*T1 
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K3-K1-1 

6049  Z(K3)-ZCK3)*T1 

XBAR -RrEAN ( Z . NSPTST) 

XVAR-RVAR(Z.NSPTST. l.XBAR) 

WRITE  ( IODEV.6056) 

6056  FORMATS. 40H  STATISTICS  OF  THE  TAPERED  RESIDUAL  DATA) 
WRITE  (IODEV.6044)  XBAR.XVAR 
6093  DO  6045  Kl-1. NSPTST 
IFT-NSPTST-K1+1 
K2-IFT+MPTS2 

6045  Z(K2)-Z(1FT> 

K1-MPTS2+NSPTST+1 
DO  6046  K2-K1.NPTS1 

6046  Z(K2) *0. 

DO  6048  K2-1.MPTS2 
6048  Z(K2)-0. 

CALL  RPSDR(Z.PSDC.Ml) 

C 

C  PERFORM  PSD  AND  PLOT  RESULTS 
C 

K3«2**(M1-1) 

DO  6040  K2-1.K3 
6040  PSDC(K2)-PSDC(K2)*10.E5 

CALL  SCAN CDELT.PSDC.-K3. 840) 

K2-2 

ZAP ( 1 )  -0 . 

ZAP(2)-1./(2.*DELT) 

CALL  SCAN ( ZAP, DELT . -K2. 480 ) 

CALL  M0DE(-8,T1,DX.T2) 

DX- ( 1 . / ( DELT*FLOAT ( NPTS 1 ) ) )/DX 
CALL  DRAW(DX,PSDC.K3.842) 

CALL  AXESC9.2.9HFREQUENCY. 18.3. 18HENERGY  *  (10)  **  5) 
L  .6251 

CALL  DRAW(0.,0., 1,8000) 

NSPTST "NSPHLD 
40  CONTINUE 
8020  CALL  DRAW(0. 0.0. 9999) 

J  .77 
C 

C  LINKS 
C 

C  UNSCALE  SCALED  FRACTIONS 
6111  OCT  0 

K5-K54K8 

K3-IBUFCK5) 

LA  K3 
STA  I VOLT 
XB-IVOLT 
XB-XB* 10000.0 
J  .6111,2 

C  LEAST  SQUARES  POLYNOMIAL  FITTING  ROUTINE 
6261  OCT  0 

IF (KALKT.NE.KALUPD)GO  TO  103 

KALKT-0 

H-l. 


S(1)-FL0AT(K5) 
DO  6216  Kl-1, 2 
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S(Kl+l)-0. 

DO  6217  K2-1.K5 

6217  S(K1+1)-S(K1+1)+ZR(K2)**K1 
6216  CONTINUE 

DO  6218  K 1-1.2 
K2-K1-1 
AAACK1 ,3) “0. 

DO  6219  K3-1.K5 
IF(K2.EQ.0)GO  TO  6219 
H-ZR(K3)**K2 

6219  AAA(K1.3)-AAA(K1.3)+ZP(K3)*H 

6218  CONTINUE 

DO  6220  Kl-1.2 
K2-K1 

DO  6221  K3-1.2 
AAACK1 ,K3) -S(K2) 

6221  K2-K2+1 

6220  CONTINUE 

DO  6222  Kl-1.2 
H-AAA(Kl.Kl) 

DO  6236  K2-1.3 

6236  AAA(K1,K2)-AAA(K1.K2)//H 

DO  6224  K2-1.2 
IF(K1.EQ.K2)G0  TO  6224 
T1-AAACK2.K1) 

DO  6225  K3-1.3 

6225  A AA ( K2 , K3 ) - A AA ( K2 . K3 ) -T 1 *AAA (K1.K3) 

6224  CONTINUE 

6222  CONTINUE 

IF ( .NOT. IGTYPE)GO  TO  103 

URITE  ( IODEV. 105)  ICALNM. AAA( 1 ,3) , AAA<2,3) 

105  F0RMATC8H  SCAN  *  . I5.3X.7HC0EF  -  .2CE14.7. IX) ) 
103  KALKT-KALKT+1 

ICALNM- ICALNFf+1 
J  .6261,2 

C  PUTS  LABELS  ON  PLOTS 
6251  OCT  0 

T1-T99+1. 

CALL  F10DE(3.T1. SAME. SAFE) 

Tl-5.75 

IFCKBARY.GT.  1  DT1-8.75 
CALL  NOTEC  1 .5.T1 .LABELC 1 ) .20) 

CALL  NOTE (SAME. SAME, 12H  RUN  NUMBER  ,12) 

CALL  NOTE ( SAME. S AME. 1RUNO.0) 

CALL  NOTE (SAME. SAFE, 14H  GAGE  NUMBER  .14) 

CALL  NOTE (SAME. SAME. IGAGENC 1, I) .4) 

CALL  MODE ( 3, T99. SAFE. SAME) 

J  .6251,2 
77  END 


PROGRAM  SIZE  -  *36577 
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1.  Tide  data  are  recorded  in  data  files  on  magnetic  tape  by  ABACS  data 
acquisition  programs.  These  data  files  on  magnetic  tape  are  used  as  input  to 
the  tide  analysis  program. 

2.  Two  general  file  formats  are  used  by  the  tide  analysis  program--a 
calibration  file  and  a  noncalibration  file.  Several  records  containing  gage 
identification  and  data  scaling  factors  are  included  in  a  calibration  file 
that  are  excluded  from  a  noncalibration  file.  Since  gage  identification  and 
data  scaling  factors  remain  constant  for  a  series  of  data  tests,  there  is  no 
need  to  record  this  information  in  every  data  file.  Thus  this  information, 
which  is  recorded  exclusively  in  a  calibration  file,  is  applicable  to  all 
subsequent  data  files  until  updated  by  another  calibration  file.  That  is  to 
say,  before  analyzing  a  particular  noncalibration  file,  the  tide  analysis 
program  must  obtain,  for  the  file,  the  correct  calibration  coefficients  which 
are  contained  in  the  last  preceding  calibration  file.  It  should  be  noted  that 
floating  point  variables  require  two  words  for  storage  and  two  alpha  charac¬ 
ters  are  stored  in  one  word. 

3.  All  data  files  are  separated  by  an  end-of-file  (EOF)  mark  with  the 
individual  file  and  records  having  the  following  format. 

a.  File  Identification  Record.  This  record,  which  is  included  in 
both  calibration  and  noncalibration  files,  is  used  simply  for 
file  identification.  It  contains  five  words  with  the  follow¬ 
ing  definitions: 


Data  Type 

Alpha 

Alpha 

Alpha 

Integer 

Integer 


Definition 


File  name  ASCII  Characters  (3A2) 

File  name  ASCII  Characters 

File  name  ASCII  Characters 

File  Type  Designator  Always  =  3 

Record  Length  of  Header  Record 

n  =  20 


Header  Record  #1.  This  record,  which  is  included  in  both  cali¬ 
bration  and  noncalibration  files,  contains  testing  identifica¬ 
tion  and  testing  parameters.  It  contains  20  words  with  the 
following  definitions: 


Data  Type 


Definition 


Not  used  by  this  program 
Not  used  by  this  program 
Integer  I  st  run  number 

Integer  Date  of  test  =  Julian  date  +  last 

digit  of  year  *  100 

(Continued) 
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Data  Type 

Integer 

Alpha 

Alpha 

Alpha 

Alpha 

Integer 

Integer 

Integer 


Integer 

Integer 

Logical 

Floating  Point 
Floating  Point 

Integer 


_ Def ini t ion 

Time  of  test 

Model  identification  (4A2) 

Model  identification 
Model  identification 
Model  identification 
Horizontal  model  scale 
Vertical  model  scale 
Reference  mean  water  level  in  pro¬ 
totype  feet  x  10  at  time  zero 
reference  readings  are  made 
Model  water  temperature  x  10 
Number  N+M  of  transducers  sampled 
for  this  test 

.TRUE,  if  this  is  a  calibration 
file 

Model-generated  tide  height  in 
prototype  feet 
Model-generated  tide  height 
Not  used  by  this  program 
Not  used  by  this  program 
Record  length  of  Header  Record 

n  =  so 


Header  Record  #2.  This  record  contains  further  testing  param¬ 
eters.  It  consists  of  50  words,  with  the  following  definitions 


Data  Tvpe 


Definition 


Floating  Point 

Generated  tide  period  in  model 
seconds 

Floating  Point 

Generated  tide  period 

Integer 

Number,  J,  of  data  samplings  per 
gage  per  tide  period 

Integer 

Number,  K,  of  magnetic  tape  data 
records  per  tide  period 

Integer 

Number,  L,  of  tide  periods  of  data 
recorded 

• 

Not  used  by  this  program 

Integer 

Number,  N,  of  tide  gages  sampled 
for  this  test 

• 

Not  used  by  this  program 

Integer 

Transducer  arrangement  number 

• 

Not  used  by  this  program 

• 

Not  used  by  this  program 

d.  Transducer  Flag  Record.  This  record  is  included  only  in  a 

calibration  file.  Its  record  length  is  two  times  the  fourteenth 
element  (N+M)  in  Header  Record  #1.  This  record  contains  two 
integer  flags  for  each  gage  sampled  during  a  data  scan.  The 
flags  are  set  up  as  follows: 


(1)  FLAG1  =  Transducer  type  +  (computer  channel  number)  x  100  . 
Transducer  types  vary  from  0  to  99  with  type  2  identifying 
tide  gages,  3  references  tide  gages,  4  temperature  gages, 
and  5  velocity  gages.  For  any  valve  besides  2,  3,  4,  or  5 
this  analysis  program  does  not  consider  the  gage  data.  In 
all  following  records,  entries  for  the  various  gages  are 
ordered  exactly  as  they  are  ordered  in  this  record.  That  is, 
the  first  flag  pair  represents  the  first  gage  in  all  follow¬ 
ing  records;  the  second  flag  pair,  the  second  gage,  etc. 
Therefore  these  flags  are  important  in  determining  what  in¬ 
formation  in  following  records  is  necessary  for  the  execu¬ 
tion  of  this  program. 


(2)  FLAG2  = 


=  2 
=  0 
=  1 
=  9 


Gage  not  used 

Linear  calibration  fitting  procedure 
Quadratic  calibration  fitting  procedure 
Spline  calibration  fitting  procedure 


+  {(number  of  calibration  voltage  readings  per 
channel)  x  100} 


0  No  entry  for  gage  in  calibration  voltage 
record  or  calibration  coefficient  record 
10,000  entry  for  gage  in  calibration  coeffi¬ 
cient  record  only 

20,000  entry  for  gage  in  calibration  voltage 
record  only 

30,000  entry  for  gage  in  both  calibration 
voltage  record  and  calibration  coefficient 
record 


The  2  (N+M)  words  in  this  record  have  the  following  definitions: 


Word  Data  Type 


Definition 


1 

2 

3 

4 


Integer 

Integer 

Integer 

Integer 


Flag  1  for 
Flag  2  for 
Flag  1  for 
Flag  2  for 


first  gage 
first  gage 
second  gage 
second  gage 


2(N+M)-1  Integer 

2(N+M)  Integer 


Flag  1  for  the  last  (N+M)^  gage 
Flag  2  for  last  (N+M)*^  gage 


Calibration  Coefficient  Record.  The  length  of  this  record  is 
obtained  by  multiplying  22  times  the  number  of  gages  indicated 
as  having  entries  in  this  record  by  Flag  2  in  the  Transducer 
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Flag  Record.  This  record  is  included  only  in  a  calibration 
file.  Eleven  floating  point  calibration  coefficients  are  in 
this  record  for  each  gage  having  calibration  coefficients.  If 
tide  gages  have  entries  in  this  record,  they  will  not  be  uti¬ 
lized  by  this  program.  The  reader  is  reminded  that  a  floating 
point  number  requires  two  words  in  memory.  The  words  in  this 
record  have  the  following  definitions: 

Word  Data  Type  _ Definition _ 

1  Floating  Point  Mathematical  fit  coefficients  for 

first  gage  with  coefficients 

2  •  • 

3  Floating  Point 

4 


23  Floating  Point  Mathematical  fit  coefficients  for 

second  gage  with  coefficients 

24  •  • 

25  Floating  Point 

26 


22(N+M)-21 
22 (N+M) 


Floating  Point  Mathematical  fit  coefficients  for 
last  (N+M)th 


f.  Gage  Name  Record.  Record  length  is  twice  the  fourteenth  element 
(N+M)  of  Header  Record  #1 .  This  record  is  included  only  in  a 
calibration  file  and  contains  a  four  character  name  or  number 
for  each  gage  scanned  during  data  acquisition.  The  2(N+M)  words 
of  this  record  have  the  following  definitions: 


Word  Data  Type  _ Definition 

1  ,  Alpha  First  gage  name  (2A2) 

2 

3  Alpha  Second  gage  name  (2A2) 

4 


2(N+M)-1  Last  (N+M)th  gage  name 

2(N+M) 

Gage  Coordinate  Record.  Record  length  is  six  times  the  four¬ 
teenth  element  (N+M)  of  Header  Record  #1.  This  record  is  in¬ 
cluded  in  a  calibration  file  and  contains  three  floating  point 
coordinates  (X,  Y,  and  prototype  water  depth)  for  each  gage 
scanned  during  data  acquisition  sequenced  according  to  order 
appearing  in  the  Transducer  Flag  Record.  The  6(N+M)  words 
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PWR i 


in  this  record  have  the  following  definitions 


Word 

1 

2 

3 

4 

5 

6 

7 

8 


Data  Type 
Floating  Point 

Floating  Point 

Floating  Point 

Floating  Point 


_ Definition _ 

X  coordinate  for  first  gage 

Y  coordinate  for  first  gage 

Water  depth  for  first  gage 

X  coordinate  for  second  gage 


6(N+M)-1  Floating  Point 


Water  depth  for  last 


(N+M)th  gage 


h.  Reference  Potentiometer  Calibration  Coefficient  Record.  Record 
length  is  twice  the  tenth  element  (N)  of  Header  Record  //2.  This 
record  is  included  only  in  a  calibration  file  and  contains  the 
slope  coefficient  for  the  reference  potentiometer  of  each  wave 
gage  calibration  stand.  This  record  is  not  used  by  the  tide 
analysis  program.  The  words  in  this  record  have  the  following 
definitions : 


Word  Data  Type 


Definition 


1  Floating  Point  Potentiometer  coefficient  for 

first  gage  having  potentiom¬ 
eter  coefficient 

2 

3  Floating  Point  Potentiometer  coefficient  for 

second  gage  having  potentiom¬ 

eter  coefficient 

4 


2N-1  Floating  Point  Potentiometer  coefficient  for 

last  gage  having  potentiom¬ 
eter  coefficient 
2N 

i.  Calibration  Voltage  Record.  The  length  of  this  record  is  deter¬ 
mined  by  multiplying  42  times  the  number  of  gages  indicated  as 
having  entries  in  this  record  by  Flag  2  in  the  Transducer  Flag 
Record.  This  record  is  included  only  in  a  calibration  file  and 
contains  calibration  voltage  readings  for  wave  gages  calibrated 
under  computer  control  prior  to  data  acquisition.  This  record 
is  not  used  by  the  tide  analysis  program.  The  words  in  this 
record  have  the  following  definitions: 


F6 


Tvv^/  —  "  1,11  11  v*  ’r 

■/. 

•. 

Calibration 

*,  ) 

\'4 

r\  - 

Word 

Data  Type 

Steps 

Gaj>e 

■  j 

r/- 

'  1 

Integer 

Zero 

First  gage  sensor 

•’i 

b 

3 

+A 

+B 

4 

+C 

■ 

:  5 

+D 

■' 

6 

+E 

7 

+D 

IB 

1  8 

9 

+C 

+B 

• 

10 

+A 

11 

Zero 

X-! 

12 

-A 

L*.\ 

13 

-B 

14 

15 

-C 

-D 

-1 

li 

16 

-E 

7 

‘-V 

17 

-D 

18 

-C 

>o 

19 

-B 

20 

-A 

U 

21 

Zero 

22 

Integer 

Zero 

First  gage  potentiometer 

23 

24 

+A 

+B 

J 

25 

+C 

.V 

,  n 

4 

26 

+D 

ii 

27 

+E 

♦ 

28 

+D 

• 

9“ 

29 

+C 

• 

S'! 

30 

+B 

v\* 

31 

+A 

1 

32 

Zero 

• 

i* 

33 

-A 

• 

34 

-B 

• 

35 

-C 

- 

36 

-D 

*- 

:>■ 

37 

-E 

“ 

TH 

38 

-D 

39 

-C 

H 

►:•:■ 

40 

-B 

L  .•• 

41 

-A 

v.-1 

K.*/, 

42 

Zero 

. 

£ 

43 

Integer 

Zero 

Second  gage  sensor 

99 

44 

+A 

£ 

• 

• 

* 

• 

^7 

V. 

,  „  th 

1 

45(N+M)-44 

Integer 

Zero 

Last  (N+M)  gage  sensor 

t* 

«T*. 

* 

* 

\ 

£ 

(Continued! 

•; 

C*\ 

‘j 

-\ 
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• 

fc' 

14 

• 

rr 

** 

1 

• 

^v/.y-y^sv-. 

Word 


Data  Type 
45(N+M)-23  Integer 


Calibration 

Steps  _ Gage 

Zero  Last  (N+M)*'^  gage 

potentiometer 


45(N+M)  Integer 


Zero 


Last  (N+M)^  gage 
potentiometer 


Zero  Reference  Reading  Record.  Record  length  is  two  times  the 
eleventh  element  (N)  of  Header  Record  #2.  This  record  contains 
a  zero  reference  reading  for  each  tide  gage  and  each  reference 
tide  gage  sequenced  according  to  order  appearing  in  the  Trans¬ 
ducer  Flag  Record.  The  2N  words  in  this  record  have  the  follow¬ 
ing  definitions: 


Word  Data  Type 


Definition 


1 

2 

3 

4 


Floating  Point  Zero  reference  reading  for  first 

gage 

•  • 

Floating  Point  Zero  reference  reading  for  second 

gage 


N-l  Floating  Point  Zero  reference  reading  for  last 

(Nth)  gage 
N 

k.  Data  Records.  The  length  of  this  record  can  be  determined  by 
multiplying  the  third  element,  J,  of  Header  Record  #2  by  the 
fourteenth  element,  N+M,  of  Header  Record  #1  and  dividing  this 
product  by  the  fourth  element,  K,  of  Header  Record  #2.  These 
data  records  are  included  in  both  calibration  and  noncalibration 
files.  Since  the  number  of  data  records  in  these  files  varies 
and  depends  upon  the  number  of  wave  periods,  L,  acquired  during 
each  wave  test,  the  number  of  data  records  can  be  determined  by 
multiplying  the  fourth,  K,  and  fifth,  L,  elements  of  Header 
Record  #2.  Each  gage  is  scanned  in  the  order  in  which  their 
flags  appear  in  the  Transducer  Flag  Record.  The  (J/K)(N+M) 
words  in  this  record  have  the  following  definitions: 


Word 

Data  Type 

Definition 

1 

Scaled  fraction 

Sample  from  first  scan  of 

first  gage 

2 

• 

Sample  from  first  scan  of 

• 

second  gage 

(Continued) 

F8 


Word 

Data  Type 

Definition 

N+M 

Scaled  fraction 

Sample  from  first  scan  of 
last  (N+M)th  gage 

(N+M)+l 

Scaled  fraction 

Sample  from  second  scan  of 
first  gage 

(N+M)+2 

Scaled  fraction 

Sample  from  second  scan  of 

• 

• 

second  gage 

2(N+M) 

Scaled  fraction 

Sample  from  second  scan  of 

# 

last  (N+M)th  gage 

(J/K) (N+M) 

Scaled  fraction 

Sample  from  last  (J^)  scan 
of  last  (N+M)t^1  gage 
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APPENDIX  G:  EXTERNAL  FORTRAN  SUBROUTINES 


’v.  -  .V  • 


1.  Several  FORTRAN  library  subroutines  are  referenced  in  this  report. 
Listings,  which  are  voluminous,  are  not  provided  in  this  documentation.  There¬ 
fore  the  general  function  of  each  of  these  subroutines  is  described  concisely 
as  follows: 


Magnetic  Tape  Input  and  Output  Routines: 

(1)  Cali  SELMT.  This  routine  selects  and  initializes  a  speci¬ 
fied  magnetic  tape  unit. 

(2)  Cali  POSMT.  This  routine  positions  a  previously  selected 
magnetic  tape  to  a  specified  magnetic  tape  file. 

(3)  Call  RDMT.  This  routine  reads  one  record  from  a  previously 
selected  magnetic  tape. 

(4)  Call  SRFMT.  This  routine  skips  one  record  in  the  forward 
direction  on  a  previously  selected  magnetic  tape. 

(5)  Call  SRRMT.  This  routine  skips  one  record  in  the  reverse 
direction  on  a  previously  selected  magnetic  tape. 

(6)  Call  SFFMT.  This  routine  skips  one  file  in  the  forward 
direction  on  a  previously  selected  magnetic  tape. 

(7)  Call  NMFMT.  This  routine  names  a  specified  file  on  a  pre¬ 
viously  selected  magnetic  tape. 

(8)  Call  WRTMT.  This  routine  writes  one  record  on  a  previously 
selected  magnetic  tape. 

(9)  Call  WEFMT.  This  routine  writes  an  end-of-file  mark  on  a 
previously  selected  magnetic  tape. 

Disc  Input  and  Output  Routines : 

(1)  Call  SELMD.  This  routine  selects  and  initializes  a  speci¬ 
fied  disc  unit. 

(2)  Call  RDMD .  This  routine  reads  one  record  from  a  previously 
selected  disc. 


(3)  Call  SAVEMD.  This  routine  saves  the  status  and  position  of 
a  previously  selected  disc. 

(4)  Call  RESETD.  This  routine  is  used  to  restore  the  discs  to 
the  identical  status  and  position  as  they  were  when  previ¬ 
ously  called  SAVEMD  routine  was  executed. 

(5)  Call  NMFMD.  This  routine  names  a  specified  file  on  a  pre¬ 
viously  selected  disc. 

(6)  Call  WRTMD.  This  routine  writes  one  record  to  a  previously 
selected  disc. 

(7)  Call  WEFMD.  This  routine  closes  an  open  file  on  a  previ¬ 
ously  selected  disc. 
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Versatec  Plot  Routines: 


(1)  Call  DRAW.  This  routine  plots  a  specified  buffer  against 
another  specified  buffer. 

(2)  Call  NOTE.  This  routine  plots  Hollerinth  characters. 

(3)  Call  MODE.  This  routine  sets  specific  parameters  for  plot 

specifications . 

(4)  Call  AXES.  This  routine  plots  one  or  both  axes. 

(5)  Call  SCAN.  This  routine  computes  scaling  factors  for  data 

to  be  plotted. 

(6)  Call  FORM.  This  routine  plots  background  grids. 

Basic  Statistical  and  Data  Analysis  Package  Routines: 

(1)  RMEAN .  This  routine  computes  the  mean  of  a  specified  buffer 

(2)  RVAR .  This  routine  computes  the  variance  of  a  specified 
buffer. 

(3)  RPSDR .  This  routine  computes  the  power  spectral  density  of 
a  specified  buffer. 


